БОТ ЭК A AKADEMIA NAUK 
ТОТ FIZYKI 


ACTA PHYSICA 
` POLONICA 


DWUMIESIECZNIK 


Orders and inquires concerning 
Acta Physica Polonica 


— complete sets, volumes and single fascicules — 
as well as other | 


Polish scientific periodicals y 


published 
before and after the war, 
regularly and irregularly, 
are to be sent to: 


Export and Import Enterprise «ЕССН” | E 


"Warszawa 1, P.O. Box 154, Poland 
Ask for cataloques, folders and sample copies. | | 


"ET 
Амь 24 


- 


ВО КАКО РЕМТА ND AUK 
рет Үт ыт ЕТЕ AY X.I 


ACTA PHYSICA 
` POLONICA 


DWUMIESIECZNIK 


|... Volumen XVI 


Жа ыла. ud ۴ Midi ы 
ВАА 2-4 
" А ^ жы » 
а b Г 
Mic. 


А” 
+. Te 


i и 


Volumen XVI 


a nn 


CONTENTS 


FASCICULUS 1—2 


Gajewski R., On Transient Radiation of a Dipole Inside a Wave 
cuide MMC M D ПБ ооа ERR RAS 
Meclewski R. and Woj da E., The Role of the Emitter Shape 
in Cold Emission Research de 
Auleytner J., Point Sources for X-Ray HL Fass ҮС ж. 
Sujak B., On the Induced Reserve Response of the Oxidized 
Aluminium G—M Counter to Visible Light, the Joshi Effect and 
the Geiger Counter Sensitivity to the- Air-Elux 2 во 9 Сс 
Frackowiak M., Decay of Phosphorescence of Rigid Solutions 
Witkowski A., Theoretical Considerations on the Possibility 
of a New Method for Separating Mixtures and Isotopes 
Sawicki J., Neutron Polarization and Angular Distributions in 
(p,n) Reactions and Nuclear Optical Model AE s 
Babecki Ј., Jurkiewicz L, Massalski J. M. and Mie- 
sowicz M., The Transition Curve of the Electron Photon Com- 
ponent of Extensive Air Showers in Lead Absorber in Thick- 
hesses between O and 25 Сш epe o е 
Budzanowski A. and Grotows ki K., Multi-plate Ionization 
Chamber for the Detection of Slow Neutrons’ . дама Bu 
Filipkowski A. Gierula J. and Zielinski By Survey 
of the Hyperfragment Experimental Data . 


Laboratory Equipment and Techniques 


Małecki J. Surma M. and Gibalewicz J. Measurements 
of the Intensity of transient Magnetic Fields by the Faraday 
Effect ЭШ АБСО» SUR OE RNC о CE илк n oC 


Letters to the Editor 


Suffezynski M. On the 3D Electron Band in Body Centered 
melo RU deu DTE Lue A SES LIB n реа RO MEN. 
SuffezyäskiM. On the Width of the 3D Electron Band in Iron 
Wotruba K., Eine neue Gesetzmássigkeit auf dem Gebiete der 
plastischen Deformation 1 Dna АЕ РОУ: 
Jaśkiewicz A. Domain Formation in Ferroelectrics 


+ 


FASCICULUS 3 


Infeld L, The Equations of Motion in General Relativity Theory 
and the Action Principle > ана, «Аы За 


F 


`135 


189 


151 


157 
161 


165 
169 


177 


IV 


Olszewski S., On the Free-Electron Theory of Absorption 
Spectra of some Linear Conjugated Systems Ны. 


Letters to the Editor 


JaSkiewicz A. Domain Formation in Ferroelectrics . . . 
Dabrowski J. and Tulczyjew B., Note on the !?C(d.p)*C 
Reaction Near the 4 MeV Resonance WO; di я ДОШ 


FASCICULUS 4 


Rzewuski H., Inwestigation of the Sulphide Layers by the 
Electric Field Pulse Method ety Ба” t 

Kolos W., Electron Correlation in Molecules I. The Ethylene 
Molecule блв c^ о E ELM 

Ко1ов W., Electron Correlation in Molecules. ПІ. The Benzene 
Molecule cU Numen US s 

Rayski J., On a Screw Model of Particles in Iso-Space 


Laboratory Equipment and Techniques 


O'Connor D. and Blinowski K., An Automating Timing and 
Recording Curcuit for Use with a Newton Crystal Spectrometer 


Letters to the Editor 


Kotos W., On the One-Centre Interaction Integrals in the Mole- 
cular Orbital Method 


FASCICULUS 5 


Czerwonko J. Asymptotic Behaviour of Higher Factorial Mo- 
ments of Electron Photon Cascades at Large Dephts : 
Królikowski W., А General Equation of Scattering Amplitudes 
of Particles in the Field Theory CENE тын сығыр Сы 
Budzanowski А, Grotowski K., Janik J., Kolos W, 
Maniawski Е., Rżany H. Szkatuła А. and Wanic А, 


Estimation of the Height of the Potential Barrier of Hindered~ 


Rotation in the CH3SH Molecule by Means of Thermal Neutron 
Scattering foc Je cl py CE Quoc Hn tuere ЕРКЕ 
Miasek M., The Application of the Tight Binding Method to the 
Investigation of Energy Bands in Hexagonal Close-Packed 
Structure. I. SET enu (n "ӨЛ ТА РА. 
Olszewski S., Interaction in the Free-Electron Model. I. Theory 
Hrynkiewicz A. Z. and Waluga T. On Nuclear Magnetic 
Resonance in Flowing Liquid . 2 E cut Бл. +. 
Krajewski Т. Pietrzak J. and Suwalski J., Investigation 
of Operating Conditions of a Dielectric Resonant Amplifier 


FASCICULUS 6 


Sredniawa B., Sur une modification de Lorentz-Invariant for- 
malisme dans la théorie de fusion 


Sredniawa B., Théorie de fusion de deux particules de Majorana 


211 


227 


231 


237 


257 


267 
279 


293 


299 


305 


327 


335 


343 


369 


381 


389 


399 
407 


М дивани О., Упорядочение тройных-сплавов, как коопера- 
тивное явление 

Bazanski S, Lagrange Е Ai NO Motion of бінен 
Particles іп General Relativity Theory. II. 

Rzewuski J., Оп Inversions of Space and Time 

Miasek М, The Application of the Tight Binding Method i5 iuc 
Trivestigation of Energy Bands in Hexagonal Close-Packed 
Structure. II. { 

Jabłoński A. Decay of ТЫЫ of lune 


Letters to the Editor 


Cofta H., Abhàngigkeit der Spinwellendispersion von der magne- 
tischen Struktur des antiferromagnetischen Raumgitters 


AUTHOR INDEX 


Auleytner J., Point Sources for X-Rays 

Babecki J, Jurkiewicz L, Massalski J. M. An Mie- 
sowicz M., The Transition Curse of the Electron Photon Com- 
ponent of Extensive Air Showers "in Lead Absorber in Thick- 
nesses Between 0 and 25 cm 

Bazanski S., Lagrange Function für 3 Motion of Chereca 
Particles in General Relativity Theory. II. D 

Budzanowski A. and Grotowski K., Multi-plate онот 
Chamber for the Detection of Slow Neutrons 

Budzanowski A, Grotowski K, Janik 7, К бов W., 
Maniawski F, Rzany H., Szkatuta A. and Wanic A. 
Estimation of бе Height of T Potential Barrier of Hindered 
Rotation in the CH3SH Molecule by Means of Thermal Neutron 
Scattering 

Go f tar Н, Abhängigkeit да Spinwellendispersion v von de ere 
tischen Struktur des antiferromagnetischen Raumgitters 

Czerwonko J., Asymptotic Behaviour of Higher Factorial Mo- 
ments of Electron Photon Cascades at Large Dephts 5 

Dabrowski Ј. and Tulezyjew В., Note on the эс(арус 
Reaction Near the 4 MeV Resonance . . 

Filipkowski A. Gierula J. and Zieliński Р, Survey 
of the Hyperfragment Experimental Data . 

Frackowiak M., Decay of Phosphorescence of Rigid E 

Gajewski R., On Transient Radiation of a Dipole Inside a Wave 
Guide СЕРИ 

ak ea hee gei ША UN Waluga T. On Nodes Magnetic 
Resonanse in Flowing Liquid 4 - 

-Infeld L., The Equations of Motion in Gener "Relativity Theory 

and the Action Principle: . 

Jabtonski A, Decay of i ына of Solutions 

Jaśkiewicz A. Domain Formation in Ferroelectrics 

Kolos W. Electron Correlation in Molecules Т The Ethylene 


Molecule . . £C 


481 


35 


119 
423 


135 


335 


481 


УТ 


Kotos W., Electron Correlation in Molecules. II. The Benzene 
Molecule „жо зз М 5951 9 90 с о. o a 
Kotos W., Оп the One-Centre Interaction Integrals in the Mole- 
cular Orbital Method аза MESE а ERO 
Krajewski T, Pietrzak J. апа Suwalski J., Investigation 
of Operating Conditions of a Dielectric Resonant Amplifier 
Królikowski W., A General Equation of Scattering Amplitudes 
of Particles in the Field Theory ug Ja АА 
Matecki J, Surma М. and Gibalewicz J., Measurements 
of the Intensity of Transient Magnetic Fields by the Faraday 
Effect Te Or Ta En а 
Meclewski В. and Wojda L., The Role of the Emitter Shape 
in Cold Emission Research Е da ooa PR 
М дивани О., Упорядочение тройных сплавов, как коопера- 
тивное явление . р dq ие 
Miasek M., The Application of the Tight Binding Method to the 
Investigation of Energy Bands. in Hexagonal Close-Packed 
Structure: LC PE P МАТ O LEN арылған 
Miasek M., The Application of the Tight Binding Method to the 
Investigation of Energy Bands іп Hexagonal Close-Packed 
Structure. II. ЖЕСЕ ee Cte ere E 
O’Connor D. and Blinowski K., An Automating Timing and 
Recording Curcuit for Use with a Newton Crystal Spectrometer 
Olszewski S., On the Free-Electron Theory of Absorption 
Spectra of some Linear Conjugated Systems к Те 
Olszewski 5. Interaction in the Free-Electron Model. I. Theory 
Rayski J., On a Screw Model of Particles in Iso-Space У 
Rzewuski Н., Inwestigation of the Sulphide Layers by the 
Electric Field Pulse Method ae Ale: 
Rzewuski J., On Inversions of Space and Time "Rn 
Sawicki J., Neutron Polarization and Angular Distributions in 
(p,n) Reactions and Nuclear Optical Model Poe ses 
Suffczynski М. On the 3D Electron Band in Body Centered 
` Structure Se ақ ным т 
Suffczynski M., On the Width of the 3D Electron Band in Iron 
Sujak В. On the Induced Reverse Response of the Oxidized 
Aluminium G—M Counter to Visible Light, the Joshi Effect and 
the Geiger Counter Sensitivity to the Air-Flux AN 
$redniawa B., Sur une modification de Lorentz-Invariant for- 
malisme dans la théorie de fusion . . . ; о, IR 
Sredniawa B., Theorie de fusion de deux particules de Majorana 
Witkowski A. Theoretical Considerations on the Possibility 
of a New Method for Separating Mixtures and Isotopes 
Wotruba K, Eine neue Gesetzmüssigkeit auf dem Gebiete der 
plastischen Deformation | 


261 


299 


389 


327 


151 


25 


415 


343 


447 
293 
211 
369 
279 


237 
435 


93 
157 
161 


49 


399 


407 
79 


165 


~~ C 


Vol. ХУ (1957) ACTA PHYSICA POLONICA Fasc. 1—2 


ON TRANSIENT RADIATION OF A DIPOLE 
INSIDE A WAVE GUIDE (II) 


By Ryszarp GAJEWSKI 


Institute od Physics, Polish Academy of Sciences, Warszawa 


(Received July 12, 1956) 


A dipole situated inside a wave guide of an arbitrary cross section starts oscillating, 
its moment being М 1 (1) sin wt. Asymptotic formulae valid far from the front of the 
forerunner, are derived for the electric and magnetic fields of different modes. An analysis 
of the pulse propagation leads to the conclusion that, on account of the pulse beco- 
ming diffused, a theoretical limit exists for shortening the pulse width. 


1. Introduction 


In the first part of this paper (Gajewski, 1956) we analyzed the field radiated by 
a dipole placed inside a doubly infinite wave guide with perfectly conducting walls. 
The moment of the dipole was assumed to be 


M = M(t) (l — e^?) sin oy t. (1.1) 


The factor (1 — е:%) had to be introduced in order to ensure conditions physically 
possible at the front of the forerunner. The resulting field intensities Е,, and Н,, 


were then found to be 


CU pp det. | 1:2 
Eres = E Бл seis А (1.2) 
1 2 uU 
Hres == о Н 4E — H (1.3) 
Wy . 
1- 1 2 33 , | 
where E, Н, Е, Н, Е, Н, are field intensities corresponding to the following current 
densities 
jı = Газо 6 (r) 1(t) cos wot, (1.4a) 
ja = I ds, д (г) (4) е“ cos cxt, (1.4b) 
‚ із- Газо д (г) I(t) e7% sin wot. (1.4c) 


(3) 
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А general solution for the fields Е (оо, о, 0) and H (оь, 0, &) radiated by a dipole 
with current density / = Ids, (ғ) 1(t) exp — i (®t — а), (where @ = wp — ig) was 
obtained by expanding these fields in terms of the eigenfunctions of the wave guide. 

In the present paper we shall discuss asymptotic expressions valid at distances 
sufficiently far from the front of the forerunner. From the technical point of view, 
we are particularly interested in the manner in which the field carried by the main 
wave becomes steady. An analysis given by Rubinowicz (1954) as well as numerical 
computations (Cohn 1953, Gajewski 1955) prove that the arrival of the main wave 
(i.e. that travelling with the group velocity) does not involve an immediato building 
up of the field of a given mode; this building-up process takes some time. We 
want to know how does the front of the main wave change its shapealong the wave 
guide. 

To obtain a clear picture of what is going on we assume that the factor o in (1.1) 
is so large, that for all practical purposes 


M = М1 (ї) sin oy t, (1.5) 


which corresponds to current density j, (1,4a). 
The E, H field radiated by such a dipole is then identical to that previously deno- 


iz 
ted by E, Н. According to Г, 83 these fields are 


E = У) (Е) + ПЕ”), (1.6) 
н = >) (Н + n н”), (1.7) 


where n is a unit vector directed along the z-axis (which we assume parallel to the walls 
of the guide). | 
Next, we have 


Ef? = (Р, i — E TR, u) Pas (1.80) 
m 1 т | 
HP =—— Om poo (1.8b) 
m 1 т l 
EP ЕТІ (Р„ ио? — km” Rm ut?) grad; Pm — 
(m) 
Ser Qm Un n X grad, Фи, (1.8c) 
umo s d UM grad, Фи — 
Lar Ks." xe ee Dre 


тер (Ри шш — km? В, ш) nx grad, Pm- (1.84) 
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In these formulae Е, and К, are eigenvalues, while Ф, and Ф, are normalized eigen- 
functions of the eigenvalue problem f,, + Л» + xf = 0 at the cross section of the 
wave guide; k, and @,, correspond to the boundary condition Де = 0, K, and Ф,, 
correspond to the condition Эд» lc = 0. P, , Qm» and А, depend on the position 
of the dipole and are given by! 


Р, = — 1 2л М: grad, On (ху, Yo); 
Qn = — i 2л [M x grad, Ф, (xo, Yo)] п, 
„= — ı2nM: n Фи (a: Уо), 


where xy, yo, 20 = 0 аге the coordinates of the dipole. 
The functions u” and their derivatives are now obtained from I (3.13) by putting 


there wy = w, and о = 0. The result is 


| 1 1 Pap elei ашық 
vem Spe а (4-2уа-е) do, (192) 
L 


1 о 1 2 Massa Ж 

= =y Бей 2 — we 1.9b 

и; оета exp — 1 С 2 Vo: 23 do, (1.9%) 
7 4 


| 1 w? 1 . B cep 
ue Lic —————— ——ехр—? wt — — Vo? — a, do, (1.90). 
0 2 


и") - 44 үс exp — i (o а ү? -а4) ао, (1.94) 
2ni) о — 00 с 
5 t 


4 2 0 z 
um => 1 : PEE Ly IDE (= — e aa ) dw, (1.9е) 
с . 
. L 


d emm E = exp —i (o — = ү c? — o? | do, (1.9f) 


The path of integration L runs from +00 to —æ in the upper half -plane of the complex 
variable о. The functions U and their derivatives are obtained from (1.9) by replacing 
Op = ck, Бу Qp = cK,. Physical significance is ascribed to the real parts of the 


functions (1.9). 


In the first part of this paper the functions (1.9) have been evaluated in the form 


of convergent series. Their convergence, however, was slow in the vicinity of the front 


of the main wave. In the next paragraph we derive asymptotic expressions for (1.9a — 
1.9f), valid at distances sufficiently far from the front of the first disturbance. 
| 


1 in formulae I (3.9) the phase factar — і was omitted. 
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2. Asymptotic expressions for u™ and its derivatives 


In order to obtain asymptotic expressions for (1.9a — 1.9f) we apply the method 
of the steepest descent. Consider the simplest of the functions (1.9), the function (1.94). 
We write it in the form 


1 do 
On) јаса gt 
Е zi 0 — Wy’ (op 
L 
where 
g- so) —i(o— کے‎ er). (2.2) 


An asymptotic expression for (2.1), valid for sufficiently large values of t, was 
obtained by Pearson (1953). The integrand is uniquely defined on a Riemann's surface 
consisting of two sheets connected along the branch cut (-о,, --о,). We define 


the first sheet as that оп which at infinity Vo? — 02 ~ о; on the second sheet we 
have at infinity Vc? — 02 = —w. The function (2.2) has two saddle points 


(2.3) 


Fig. 1 


2 с 
where В = ей At these points, the path of steepest descent (Fig. 1) crosses the real 


: 3 ; 
axis at angles + 27: The coordinates of the points P, and P, at which the line crosses 


the branch cut are 


+ = Vi SR | (2.4) 
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We replace the integral (2.1) by a sum of two integrals taken along АВ and CD. Deno- 
ting 


I. = iE ет do 
ri 0 — 0, ' оро 
АВ 
1 do 
І, = — 
от fe 0 — 0 ' сш 
бр 
we obtain 
u ДЕБ RO, (2.6) 


where RÍ? is the residue of the integrand at the pole cx, thus 


КЇ? = 1 (t—t) exp- i (wt - = Yo? — 92), (2.7) 
where 
Wo >05 (2.8а) 
foe one (2.8b) 
ЙЛ ОР 225 / S WIPE 
Voter 


Asymptotic expressions for the integrals (2.5) can be obtained by simple calcu- 
lation, typical for the method of steepest descent. These expressions are 


1 1 p AOS ) 
jbena E gt dw == — 0) (2.9) 
i RE к ү227 pesce үгү 
Om 


оли 


where 


N = o V1 — В? 3 (2.11) 
Thus, according to (2.6) we have 
Ê е beatae 4) م‎ AD + RP sv 
НР ен оу n 
`, т 


(т) ш 


# 
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This result, valid for sufficiently large values of t (or n), was obtained by Pearson (1953). 
It can also be obtained by substituting into I (4.26) — (4.28) asymptotic expressions 
for the Bessel functions: 


2 2 1 
a 


(Rubinowicz, 1950). 
The remaining five of the six functions (1.9) can be evaluated in a similar manner. 


As an example we evaluate the function u”: 


do 


um = — ic st == 


: е == 
Va w O Oo 


ЕЕ + Кот?) — 
ee am. Vo 2— o? vaca J^ x 
AB 


l 
c|——— 7, + of E, Ко», (2.13 
Ben dic ot 1 Бі 


where 
| 1 : 2 
Rm — —ic pm a 1 (t —to) exp- i С Ve - a2), (2.14) 


Bearing in mind that according to (2.3) 


saat = al eis 
Vat at ава 
2 m үс? В? 
we obtain 
. ү1— 8 
um) — —i nere (Т, — L) + Ко, k (2.16) 


All the functions (1.9) can be expressed i ina a similar way, in terms of J, — I, or I, + Ij. 
The final result is 


XD ; 
et — 4) + ВА —o =) |. 


^ (2.172) 


= (1, E: 4; EE -2)) (2.17) 


um — 


ортаа E dea уа == |. 


(2.17с) 


-— "YN 


On Transient Radiation of a Dipole 9 


ug = ll Y (t — t) exp =i (or = үз), (2.174) 
о) i| e. VERSETS s. 
Ul. = » Га. (I, — I) + Of оз, 1 (t — to) exp - i (oer — = Va2 =) A 
(2176) 


w 


ЖЕ =—ı Г.) js Wo 1 (t — to) exp- i C - E Vo2 — 3)| 5 (2.17f) 


It can be easily seen that with э approaching infinity, both J, and J, tend to zero. 
Thus the residues represent the steady state solution of the problem. Keeping in mind 


the convention concerning the sign of the root Vw — w? (cf. footnote in Part I) 
we find that in the case of wy < w,, the corresponding mode becomes exponentially 
damped in the steady state. 

It may be noticed that in the case of а) > о, the denominator in (2.9) is equal 
to zero for t = tọ. Thus, at the front of the main wave I, becomes infinite. This means 
that, in the vicinity of the wavefront, asymptotic formula (2.9) is not valid. However, 
following a method developed by Pearson (1953), an asymptotic expression for J, can 
- be found which remains finite at the wavefront. The idea of Pearson’s method is that 

we do not substitute simply w = ©; into the denominator in (2.5а) but consider 
€) as variable along AB. After performing the integration (cf. Appendix), the following 
result is obtained \ ; 
= O OD e (2.18) 

ІШ V2 | 


where 
[ 7] Wo : 
= — [1 — — y1— р? 2.19 
| d Z m | ut В ) (2.19) 
so that d>0 if t>tp andd<0 if t< tọ; C denotes the Fresnel integral 


со 

€ (d) = f 

5 Addi | 

Substituting (2.19) into (2.17) we obtain expressions which are continuous at the front 

of the main wave. In this way, the process in which different modes become steady can 
be followed in detail. 


ez” dy = |€ (Al eG + 37). (2.20) 


3. Deformation of the pulse front 


Inside the wave guide let us choose a cross section z — const, and let us suppose 
that only one mode; say that corresponding to cut-off frequency w,,, can propagate 
. undamped. The first disturbance arriving at our cross section will be the forerunner 
travelling with the velocity c. Its structure is described by formulae I (4.20), (4.27), 
(4.28). It can be seen’ (cf. Rubinowicz 1956) that the amount of energy carried 
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by the forerunner rapidly decreases with the distance from the dipole. After the fo- 
rerunner has passed, we observe, in our cross section relatively weak oscillations until 
the moment to = 2/v,,, i.e., until the arrival of the main wave. This wave does not, 
however, arrive instantaneously but is, in a certain way, diffused and deformed when 
compared with the original rectangular shape of the pulse front envelope. Thus, the 
front of the main wave can be looked upon as a deformed pulse front. In this seetion, 
we shall examine this deformation. 

We take the simplest of the functions (2.17) namely, the function u,. Physically и, 
represents the axial component EC? of a ТМ wave (this follows from 1.8a with №, = 0). 
It can be shown (Kulinski, private communication) that the remaining functions 
behave, if sufficiently far from the forerunner, in a way similar to U. 

According to (2.6), (2.10) and (2.18) the function u, can be described in terms of 
an asymptotic expression valid sufficiently far from the front of the forerunner 


EE 


IE (19| л 
z= 4 —— — y (|4 — ال‎ 
үз cos (m х ( е е 
+ l(t— to) cos (wot — Ф), (3.1) 
where i 
© 
Po = y үш — Om» (3.2) 
% 
b or (3.3) 


The upper sign in (3.1) refers to the instants of time t < tọ, the lower to t > ty. 
At the moment of the arrival of the main wave we have 


2 5 
Сы 
Алты ЫК V5?—1 
Han b , 
00 КРЕС 1 
7] No On to 1 65 = c b b? = 1 (3.4) 


Noting that ==, corresponds to d = 0, and that |€(0)| = = we find that, - 


in the vicinity of 4, 


IE (14) КЕРЕ; 1 
үз enun TS gu 
provided that 
1 (52—17? 
z> 1л? en №, (3.6) 


where A, is the free space wave length corresponding to the frequency шу. 
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Condition (3.6) is practically always satisfied, and thus we infer from (3.5) that 
the second transient term in (3.1) can be neglected when compared with the first one, 
at least for £ = % (4 = 0). It can be seen however, that this holds true also for 4 far 
from %, since the first term in (3.1) can then be replaced according to (2.9) by 


4 


Comparison of (3.7) with the second term in (3.1) leads to the conclusion that, except 
for the case В = 1, which is beyond our intérest, this second term can always be neglec- 
ted. Thus, (3.1) can be written as 


d 
Uz = c Би соз ( — y (1]) + я --1(--%) cos (wg t — фо). (3.8) 


In order to examine the behaviour of u, in the vicinity of tj, we develop 7 into 
a Taylor series about the point £ = &. Introducing 


8 cos (n >= я. 1 А | 4 | (3.7) 


V2 an 1—bV1— BB — VEmd чел 


O=t—t, (3.9) 
we have 
a 97 1 [927 4 
7 ғғ No + (22), 9 + 9 (2 3) д ’ (3.10) 
so that 
Wa t 1 z (ЬЪ®— 13 
= olo ر‎ at کے رھ ج‎ I ра É (3.11) 


Here, we confined ourselves to the first terms of the expansion thus assuming 


1 [937 DE 
= EJ 0 3E (2: 92|, (3.12) 
which leads to 
joi tor ае quas (3.13) 
c b(b? — 1) pus 
We introduce a time dependent "frequency" 
Qul 906 و‎ 18 (3.14) 


гае 


and note that (3.13) implies 2 < wọ. 
л л 
` The function x (|d|) grows monotonously from x (0) = {| to х(20) == and can 


be considered as a phase shift weakly depending on time. According to (2.20) this 
phase shift is 


x (19) = arc tg т —— —5 4 (3.15) 
2 
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where S and C are the tabulated Fresnel integrals 


x 


5 = | sin tay 
0 


х 


C (x) = | lr dy. 
0 


The behaviour of the function и, can now be pictured as follows: 
When t < tọ, d tends to zero with t — tọ, and thus |€(|d|)| increases from almost 


1 
zero (for t < ty) to үз (for t = %). The function и, is, for t < ty 


€ (4 Wy t 7 
Ue = = cos |“ + 00+ pE — y (|d]) + A (3.16) 
We see that the “frequency“ of the oscillations is higher then wp . It decreases with 
time and becomes exactly w, at the moment t = ty. The amplitude increases with time 
and reaches, at t = tọ, the value 1/2. 
When t > fo, d increases with time from zero to infinity, hence, |С(|4))| decreases 


1 k 
from үз to zero. We can now write the function и, as 


—. аа, 


и, үз 


% % 


08 |в-Ф 93 ЫА (141) + Я + cog (mt + o È — фу). 
(3.17) 


Thus, we have a superposition of two oscillations with slightly different frequencies, 
Wy and wg — £2. This results in a certain type of beat characterised by decreasing ampli- 
tude and increasing beat frequency. Keeping in mind that 


Wo Lg 


52 + Po == Wo by, Ер (3.18) 


(3.17) can be written as 


“= | / 1+ у? — 2 у cos [ t ш — 7] cos (о; 8 — д), (3.19) 


y = ва 1. 
Vig үг 


and 6 is a phase shift slowly changing with time. 


where 


(3.20) 


PE 


зоа 


2 


с0а) | 


бекке 


On Transient Radiation of а Dipole 3 


In the vicinity of t = % we can expand d into a Taylor series. Noting that 


Sd mu T 
SY E Wer ла b% 


aed) 
2/2; ағасы 
әза er 3 e Wy c (b? — 1)*№ (b? — 2) 
E 4 Sy s ЖА , (3.21) 
D оу c (b? — 1) 
4-1/% E 9 (3.22) 
1 (ad да 
= er 3 et 
31 Ca (24) ә 


which is satisfied if the condition (3.13) is fulfilled. From (3.13) and (3.22), we 
find that our approximation (3.22) holds for d such that 


? 


we have 


provided that 


<< ; (3.23) 


/2z (b? — 1) 
ES E T (3.24) 
In this range, we have, according to (3.14) and (3.22), 
Q8—— d, (3.25) 


so that (3.19) becomes 


и; = / l + у? — 2 у cos EPA cos (а, 0 — 6). (3.26) 


According to (3.16) and (3.26) the envelope of the oscillations is 


л 


л 
1 + y? — 2 у cos E d? + 100) | ии. а 
j ical 


|ші| = 


For x >> l we have asymptotic formulae 


1 1 
50) р с 3 
Jean. | Spot, 
zac Re.“ 3.28 
Cc ера (3.28) 


LTE 
Their accuracy is better then 1% for x > 3. Thus, for |d| <3 we have x = = and 
| í 
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H : А 
у ^ —=——. Since y is much smaller than unity, we can expand the quadratic 
V2 z |d] 


root in (3.27). The result is l 
Ez Visa Не 
(3.29) 
t < te. 
Fear 2л |4] ‚ 


We see that under the assumption d = 3 the transient term in (3.8) becomes identical 

with (3.7) which means that the asymptotic expression (2.9) for J, fails in a region 

where |d| <3 and must be replaced there by the more accurate formula (2.18). 
We have so far discussed the expression (3.8) for и, treating it as а superpo- 


[и] Fr 


sition oftwo terms oscillating at slightly different frequencies. Another way of treating 
(3.8) is to make use of (2.20). It is now convenient to rewrite (3.8) in complex form 


ps ME СЕИ L(t бетон). (3.30) 


We then have for t < t, 


= T. € (Id])e? (ara). (3.31) 


Taking into account that 


l+i 
буе T : (3.32) 
where 
F( = f é3* dy = CQ) +15), _ (3.83) 
0 
and noting that F(—x) = —F(x) we obtain for t < tọ 
Tue Es a г Еа) etse). (3.34) 
ү? | 
Substituting for 7 the о кый (3.11) and using (3.25) апа (3.18), we get 
| 1 + i (ә = Ф + =) | 
и„ = — Е е i Я 3.35 
„л; (3.35) 


so that 


ш- "e Lg Fa, = x] / E & ca т Б a sa| . (8.36) 


It can easily be proved that (3.35) and (3.36) hold also for t > tg. These formulae have 


Nd A 2. L um 


قوفف قق قق وف 
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been derived in a different way by Ginzburg who inrestigated pulse propagation 
in the ionosphere (cf. Alpert, Ginzburg, Fajnberg, 1953). They follow from (3.1) as 
an approximation, provided conditions (3.6) and (3.13) (or 3.24) are satisfied. The 
function (3.36) is plotted against d in Fig. 2, its plot being, of course, identical with 
that of function (3.27). 

The reason that we decided to cite the *beats method* is that it allows for an 
interesting interpretation of the process as a superposition of a steady main wave with 
a decreasing transient wave. Another reason is that with this method a simple expression 
can be obtained for the frequency change in the range t < t, (cf. 3.14). 


4. Some practical conclusions 


Our results allow some practical conclusions to be drawn. 

a) “Surge voltage“ 

It follows from the curve of Fig 2. that the first field maximum which appears 
for 4 = 1.4 is about 1.17. This means that the highest voltage inside the wave guide 
does not exceed 117% of its steady state value. 

b) Transient time 


Fig. 2 


We note that (3.22) can be written as 


Я =т 4, (4.1) 


where 


Az b^ 


Ty ape 6—1). 


(4.2) 
may be called the standard transient time. p 
The actual transient time can be defined іп a number of ways. If we define it as 
the time required for the envelope to grow from, say, 10%, of its steady state value 
4 
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(de —2) to the-steady state value 1 (d~1) then in the formula At = т: Ad we must 
put Ad = 3 so that with such a definition - 


At = At’ = 3r. (4.3) 


If we define the transient time as that required for the envelope to grow from 1% 
of its steady state value to such a value that the oscillations of the envelope do not 


exceed п%, then, according to (3.29) 


; | 
m: 7° 4.4 
V2 x\d| % (up 
and 
Ad = 2 1 ЗЕ (4.5) =. 
ST n | UN 
so that 
Area И (4.6) 
лп, n 


Thus, e.g. with n = 5% one has Zt), = 9. 
с) *Frequency* of envelope oscillations. 
This frequency can be determined frotn (3.29), (3.25) and di 1) as 


d) Numerical example 
Let us take 


№ = 3 cm, de DEO d a y 

In this case inequality (3.6) is ус fulfilled. - 

| м, _ We check the range of d w 
За m can Lbe иннаа eee 


Ww 
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апа 
Дів, = От = 9 х 10-8 в, 


respectively. 


We now calculate the instantaneous “frequency“ of the envelope oscillations at 
the moment corresponding to, say, d — 4. We have 


0-54 а тақ 


Thus, the frequency (2 at that moment corresponds to a wave length 3 m. 


5. Pulses of finite duration 


All the considerations presented so far correspond to the case for which the 
dipole moment depends on time through the factor 1 (t) sin оё. In practice, rectan- 
gular pulses of finite length are used. In that case our results hold in the vicinity of the 
pulse front, provided the pulse duration T is sufficiently long. 

It is, however, possible to obtain results valid even for short pulses by treating 
the dipole moment as depending on time through 


l(t) зто — 1(t — T) sin w(t — T). (5.1) 
Here Т must satisfy the condition 
ФТ = 2лп (5.2) 


where n is an integral (this is to ensure the vanishing of (5.1) for t7 T). 
Hence, the solution for a cut-off pulse is obtained as a superposition of two so- 
lutions for pulses of infinite length, e.g. 


и (t,z) = u, (t, z) — u, (t—T,z). (5.3) 


All the remaining functions occurring in (1.8) can be obtained in exactly the same 
way. | 
If asymptotic formulae are used for и, and (3.6) is satisfied, then, in the range 


_ determined by (3.13) or (3.24) we have, according to (3.35), (4.1), (5.2) апа (5.3) 


ие ае еа). (5.4) 


сак ү T 
Mir ТШ : қ 
This expression approaches (3.35) if DEM 2 > 1 which can be written as, е. р. 


UE UNT 
T 


This means that the influence of the rear of the pulse can be neglected for à such that 


9 « T— 10r (5.5) 
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0 ا‎ Е м 
The pulse envelope is described by 


1 д TA 70 
«enn 


ae: 
As we remember, ® denotes the time counted from the moment & = z/v,,. 
When working with pulses of a finite length it is more convenient to count the time 


(5.6) 


from the moment tọ + m 


Introducing a new variable 


ТЕ ; 
ES ty арлы (5.7) 
and denoting 
3 = 2: = р (5.8) 
т D 
we have from (5.6) 
1 4 4 


Thus the pulse envelope is symmetric with respect to the middle of the pulse (D = 0) 
and its shape can readily be found by using of Cornu's spiral. It is possible to estimate 
the pulse duration Т for which the influence of the rear of the pulse on the shape of 


L 
the pulse front can be neglected. Substituting 9 = — into (5.5) we obtain 


2 
T >> 20т, or A — 20; (5.10) 
We shall now pass to the case of very short pulses. Formula (5.9) сап be rewritten 
as 
э | 
- 5 re 4 
jus = Teer) c rr DE ау {бый 
v 2 
pee D 


Suppose now 
A 
D> =. (5:12) 4 


In that range of D the first term іп the exponent can be neglected, thus 


y2 


4 
c 38:515 іл Dy sat .00390 ; 
CREME dy = Tp sng DA. (5.18) _ 
5 
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It can be seen that in the case of very short? pulses (4 <1) this formula provides 
an approximation valid in the full range of D (except for excessively high values of D 
where (4.8) would not be satisfied; however, this limitation is of no practical im- 
portance). 

Curves computed from (5.9) with A = 2, 1, 0.8 are plottedo n Fig. 3a, b, c. We see 
that for A < 1 the shape of |u| is approaching that of (5.13) 


Im 


= 5 
Е къ“ т. T e 


Fig. 3b 


? An example od such short pulses are millimicrosecond pulses described by Beck (1956). 


90 В. Gajewski 


The maximum value of the short pulse envelope сап be found from (5.9). Putting 
there D — 0 we get 


Іші шах = V2 | F (4) (5.14) 


A 


For A <2, | E N ES the accuracy of such an approximation being about 


10% for 4 = 2 and about 1% for 4 = 1. Thus, the maximum value of the pulse 
envelope is, for very short pulses 


A 
Ш = = (5.15). 


It can be found (e.g. from the Cornu spiral) that the maximum of the function ІЕ(«)| 
occurs for x = 1/2, and is equal to about 0.95. Thus, in order to obtain in the middle 


Fig. 3c 


of the pulse the highest possible maximum, the relative width of the applied pulse 
should be A == 2.4 which corresponds to |u| max = 1.34. 


The relative width о of the transmitted pulse can be defined as the distance be- 
tween the first two symmetrical minima of the envelope (cf. Fig. 3c). Thus, from (5.13) 


к=п (5.16) 
so that 
a-A=4 (5.17) 


Passing from relative widths A and & to actual widths (measured in, units of time) 
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T = tA and. A = та we obtain 
Ld? (5.18) 
From (5.15) and (5.18) it can be inferred that there is a theoretical limit for the 


shortening of pulses transmitted through wave guides. If the pulse duration T becomes 
too short then, after passing the wave guide, the height of the pulse may be quite small 


and its duration A may become much larger then T. 


From (5.18) it follows that 
=); (5.19) 


so that with 4<1 we have A> 47. 
The above conclusions may be of practical importance in short pulse microwave 
technique. 


€ 
Numerical example 


Let us take the values used in $ 4 во that т = 1 х 10 8s. We shall find the pulse 
duration required in order that the rear of the pulse not influence the shape of 
the pulse front. From (5.10) we get 


Т.->2 1071 == 0.2 us 


Formulae (5.13) and (5.15) are valid for 4 < 1 which corresponds to.pulses with 
Т< 1х 10-8 s. If we take T = 0.5 x 10-8 s (A = 0.5), then the width of the transmitted 
pulse will be, from (5.19), А = 8x 108 s. 

In order to obtain the highest possible maximum in the middle of the pulse, 
envelope, the pulse duration should be T = 2.4 т = 2.4 x 10 * s. 


2 


6. Concluding remarks 


In this paper we were interested in the behaviour of the field at instants of time 
close to £ and removed from the front of the first disturbance. Thus asymptotic 
formulae (2.17) have been used throughout. However, we may point out some problems 
where we are interested in knowing the field at the forerunner. This will be the case 
if we wish to estimate energy losses caused by the transient propagation of those 
modes which.are exponentially daniped in the steady state. In that case exact 
solutions discussed in Part I should be used. 

It should be added that all conclusions concerning the behaviour of the function 
u, remain valid for a function Р (t,z) describing any kind of one dimensional wave 


motion governed by the equation 


1 оф, 
cd ur Ба ЕКА. 
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with the boundary conditions 
Е (t, 0) = 1 (t) cos wot, 
F (0, z) = 0, 
F,(0,z) = 0. 


Examples of such motion are: propagation of cut-off trains of de Broglie waves 
(Rubinowicz, 1950) and propagation of cut-off pulses in the ionosphere (Ginzburg, 
cf. Alpert, Ginzburg, Fainberg, 1953). 


The author wishes to express his gratitude and indebtedness to Professor A. Ru- 
binowicz for encouragement and for his most helpful advice. 


APPENDIX 


We wish to find asymptotic expression for the integral 


TA 1 g(o)t do 
а= ] = E 


where 
g(0)=—i (o — 2 Мо? =a) 


valid for sufficiently large values of t. We shall use the method of steepest descent. 


Saddle points corresponding to g’(w) = 0 (prime denotes the derivative =) are 
о 


= Om Om 


w, = ———. ОЕ === | 
1 үл — 8° ; w = үт fi (А.2) 
so that on the path AB (Fig. 1) we have approximately 
~ =. 
g (w) = (93) + 5 8” ©) (o ~e (A.3) 
where 
КІСЕ т 
6 (@1) = — іо, М1 2; в Gham toa - (A4) 
Thus 
ix is сз ”ي‎ @) (o za 


I= о, (А.5) 


TEN 
AB 
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where 7 = w,, t V1 — В. Substitution of о = à, into the denominator of (A.5) leads 
to formula (2.9) which fails in the vicinity of t = ty. To avoid this difficulty, we con- 
sider c variable along AB. This assumption yields 


+00 
I E ere) | o (A.6) 
V2 лі VA Op V2 yer a, 


8" (61) 15? = =- 2°, (A.7) 


қ ig 
= : doe Pl ME "NM (A.8) 
m: z —V g” (61) (Фу — о) 


I I en = 42 
КОРУ |Р дез зар т UN (А.9) 
Vom B 
E í : 


Substituting 


— BP)’ t 1 7] Wo ery 
n= FE (y — og: a= ЕТРЫ ze а s 


(А.10) 
we get | 
: л л а dz 
ғыл кізе; м پچ‎ AM). 
= 16 2 fe 2и а | ( 
P : 


; у ‘ ( 
We evaluate this integral using a method analogous to that used by Rubinowicz Dur 
to find the diffraction field near the edge of the shadow. The idea of this Se = 
to multiply and divide (А.11) by such a factor as to obtain in the dpa же E. 
integrand a factor z? — iu?, equal to the denominator of the integrand. Mu tip ying 


л 


(А.11) by е 52 “ and езі” we get 


i | л л л а? $ .. еце | dz A.12) 
то [ейенен а 
ae І z2—iu l 
| p 


2 
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Our next step is to isolate the integral. Thus we have 


z d 


; zx ut Gr Eu dz 
pit es t) v ge ug 
п 


22 — iu? 


(A.13) 


In order to get rid of the singularity in the integrand we differentiate both sides with 


respect to |d|. Integration with respect to 2 then yields 


д та el ОЕ 
2414 4542 = odin ae | (А.14) 


This can be integrated with respect to |d|, the result being 


I; EY ME Jii а — * (А.15) 


With the help of (2.20) we get the final result 
E d (9 (Al FR i (n + х ар + =) А (А.16) 


d Eld) («++ е) 160 


= — а ya d үз 


КРАТКОЕ СОДЕРЖАНИЕ 


Р. Гаевский. Радиация диполя в волноводе в неустановившемся состоянии (11). 


Диполь, помещённый в волноводе любого сечения, начинает вибрировать, причём 
его момент зависит от времени через функцию М 1 (t) sin ct. Выведены асимптотические 
формулы для электрического и магнитного поля разных типов, важны, вдали от 
фронта первого возмущения. Анализ процесса распространения импульса приводит 
к выводу, что, ввиду размыва импульса во время прохождения через волновод, сущест- 
вует теоретический предел для сокращания длины импульса. 
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THE ROLE OF THE -EMITTER SHAPE IN COLD EMISSION 
RESEARCH : 


Ву Rvyszanp MECLEWSKI AND LucjAN WoJDA 
( Received June 29, 1956) 


Observations were made of the changes of electron patterns from iron and nickel 
points, these changes taking place under the influence of changes of temperature and 
voltage. It was established that the emission of electrons from weakly heated points takes 
place from the corners and edges of the monocrystalline point. 

Аз a result of surface migration or plastic deformation, there occurs on the corners 
a piling up, which causes a local increase in the electric field intensity. As a result, the 
distribution of the current in the electron pattern maps the shape of the surface of the 
monocrystal forming the end of the point. 

The observed effect permits one to assume that the different values for the work 
function obtained from the surfaces of the monocrystal by various authors result from 
various degrees of rounding of the point. 


Introduction 


The first quantitative attempt to investigate cold emission and also the first 
work whose aim was to check experimentally the Fowler-Nordheim theory already 
met with difficulties because of the impossibility of exactly determining the electric 
field intensity at the surface of the emitter. The existence on the cathode of an irregular 
submicroscopic unevenness did not allow the electric field intensity to be calculated 
from the macroscopic geometry of the system. After a number of attempts by different 
investigators, this difficulty was overcome by F. W. Müller. For the cathode he used 
a metallic point, but in contrast to others, he subjected it to strong heating in vacuo, 
as a result of which it became, in his opinion, atomically smooth (Müller 1953). 
However, some authors (Benjamin and Jankins 1940, Haefer 1940, Asadullin and 
Shuppe 1954, Spivak and Gelbert 1954) showed that even after being ireated in this 
manner, the point does not remain perfectly smooth. Thus, also in this case, although 
to a sinaller degree, there appear local increases in the electric field intensity. The 
lack of smoothness of the point is caused by the formation of corners and edges on 
the small monocrystal at the end of the point when there is a strong electric field at 
its surface. Under these conditions, the density of the emitted current represented 


(25) 
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by the Fowler-Nordheim formula 


3J 
i = АДР) Ete- PF 


(iis the current density in A/cm?, E is the electric field intensity in V/cm, W 
is the work function, А = 6.2 х 10-9A/V B = 6.85 х 10*V-^A-1, КИМ) = 
== DM NUN ү-і 
(п БИ) Wa 

u is the electron energy on the Fermi level) depends on three factors. Firstly, it 
depends on the work function of the electrons from the metal, this being a function of 
the crystallographic direction; secondly, it depends on the purity of the point surface 
(adsorption changes the work function); and thirdly, it depends on the existence of 
strong local electric fields’ оп the corners and edges of the monocrystal of the emitter. 
In the first stage of the investigations, when it was assumed that the emitter is perfectly 
round, the current distribution was explained primarily by the first two of the factors 
mentioned above. In later work, when it was seen that the point is not smooth, it was 
shown that the electron pattern (current density distribution) to a large degree (An- 
dreev 1952) and sometimes exclusively (Meclewski and Wojda 1955) depends on the 
third factor, that is, on the shape of the point. The hypothesis of Smirnov and Shuppe 
(1952) that different authors obtained different values of the work function from the 
same crystal planes as a result of different degrees of rounding of the emitter should 
be regarded as correct. In order to determine the work function without error, it is 
necessary to know the electric field distribution at the surface of the emitter, and for 
this, in turn, it is necessary to know the shape of the emitter. With this purpose in 
mind, investigations were begun to determine the shape of the point under different 
conditions by means of an electron-lens microscope (Haefer 1940, Benyaminovich, 
Smirnov, and Shuppe 1953). 

The influence of the geometry of the point on the electron pattern is difficult 
to observe for a metal with a high melting temperature. In order that this influence 
be the predominant one in this work (which is an extension and supplementation of 
the letter of the authors published in 1955) the points used were made from metals 
having a low melting temperature, namely, nickel and iron. 


Equipment 


The work was conducted with the Müller microscope which is commonly used 
for this purpose. The bulb was made from heat-resistant glass in order that it may 
be heated to a high temperature for degassing. A point of 1 to 2 mm. in length was 
spot-welded to a molibdenum loop, the top of which was thinned to about 0.15 mm. 
The use of a molibdenum loop permitted the point to be heated to a temperature 
close to the melting point. Both the nickel and iron points were electrolytically etched 
by alternating current in а weak alcohol solution of hydrochloric acid. By means of 


adia 
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an optical microscope, the radius of curvature of the etched points was found to be about 
1 micron. The work was conducted chiefly in a dynamic vacuum because of the ease with 
which the points would be changed. In addition, obserwations were made with sealed off 
bulbs because of the good vaccum which сап be obtained-in them with the help of getter. 
А dynamic vacuum of an order of 10-6 mm. Hg was obtained by means of a three- 
stage mercury diffusion pump with a freezer. А better vacuum was obtained in the 
sealed off bulb than can be measured with an ordinary ionization manometer. The 
entire bulb was heated several hours, before use, to a temperature of 500 °С. The 
loop, together with the point, was heated to a temperature close to the melting temper- 
ature of the point material. The voltage necessary to obtain electron patterns (Fig. 2a, 
3, 4) was, depending on the radius of curvature of the point, 5—6 kV. The current 
which then flowed was 3—5 wA. The last stage of the patterns (Fig. 2f) was obtained 
with voltages higher than the initial voltage by 300—500 V, the emission current 
increasing to 30—40 LA. 


Results of the Observations 
А. Nickel 


Nickel has a face-centred cubic network. The simple model of a monocrystal 
of such a network, with planes of the lowest indices, i. e. (001), (011) and (111) 
is shown in Fig. 1. Comparison of the electron photographs (Fig. 2a—f, 3, 4) with 
this model permits a one-to-one correspondence to be established between the areas 
of emission and the corners and edges. 

Changes of the electron pattern’ which arise as a result of changes of voltage and 
temperature of the point are shown in Figs. 2a—f. The point, heated for several 
hours at temperature of 1200? C, was not cooled to room temperature, but to about 
500 to 550? C and the voltage was applied very cautiously. The electron pattern 
arising under these conditions is shown in Fig. 2a. A very small, gradually increasing 
voltage causes changes in the patterns as shown in Figs. 2a—f. Similar pattern changes 
can be obtained with a fixed voltage and with a corresponding increase in the temper- 


_ ature of the point by increasing the current flowing through the loop. In Fig. 2a is 
outlined the edges surrounding the centrally situated plane —(111). The brighter 


spots correspond to pairs of corners of the monocrystal situated near each other (these 
are the pairs of corners A—B, C—D, and E—F shown in the model in Fig. 1). In 
Figure 2b the edges and corners are now entirely distinct. In Fig. 2c can Бе зееп 
similar pairs of further peripheral regions of the monocrystal (these аге the pairs of 
corners G—H, I—J, K--L, etc. in the model in Fig. 1). As seen from Figs. 2d—f 
the corners become more and more distinct. Fig. 2f represents the final stage of the 
pattern, and a further increase of voltage can cause the destruction of the point. In 


_ this photograph the corners which in the previous patterns formed pairs are now seen 


to be separate (e. р. corners А and B are now mapped separately). 
| D 
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During the observations, a ’drifting’ effect was frequently noted, that is, when at 
a certain moment the continual rise in voltage was stopped, the pattern 
continues to change, as if under the influence of a very slow increase in voltage. 

Numerous experiments (about 50 points were investigated) showed that the 
orientation of the monocrystal forming the end of the point is random, unlike the 
case for tungsten. Thus in Fig. 2a—f it can be seen that the central plane is of the 
type (111); in Fig. 4, it is a plane of the type (001); and in Fig. 3 none of the planes 
is centrally situated. A statistical study of the negatives showed that the probability 
that a given plane occupy the central position is the same for all planes. The manner 
in which the pattern forms is the same, independently of the orientation of the mono- 
crystal: first there arises the central part, and then, after an increase in voltage, the 
peripheral parts of the pattern appear. 

The large bright spot in Fig. 5 was sometimes obtained on the screen if the 
voltage applied to the point was somewhat greater than that used to obtain the fully 
built-up pattern (е. g. Fig. 2f). After the voltage is removed and the point reheated 
for about 15 minutes, the pattern obtained is the same as that obtained before the 
spot in Fig. 5. arose. Sometimes, however, the pattern is distorted to a greater or 
lesser degree. 6 

If the applied voltage is too large, then the pattern either vanishes immediately, 
or there appears, for a short time, the above-described spot, which subsequently 
also vanishes. In the first case, it appears that the point is destructed;in the second, 
its end is melted as a result of the Joule heat. Observations under the (optical) micro- 
scope show that in the latter case, at the end of the point, there is a sphere of a large 
radius of curvature, larger than the radius of the base of the point. Attempts to obtain 
a pattern from such a damaged point in the usual manner, with voltages up to 12 kV, 
did not give any positive results. Various attempts were made to regenerate the point. 
It was sometimes possible to do this. The voltage should be applied to a rather strongly 
heated point (about 1000? C, and when there appears on the screen a large bright 
spot (similar to that in Fig. 5) the point should be rapidly cooled by shutting off the 
heating current. Next, the point should be strongly heated without the field. After 
the temperature is lowered to about 500? C and the voltage applied, a pattern is 
obtained in the majority of cases. This pattern is, in general, distorted, which also 
sometimes indicates that the monocrystal is oriented in a different way, or that 
several patterns are superimposed on one another. 


B. Iron 


The points used for the investigations were made from steel, carbonyl iron 
(99.9%) and from carbonyliron with 1% addition of silicon. As in the case of nickel 
points, the iron points were heated to a temperature of 1200? C and then cooled to 
about 500? C before the voltage was applied. 

The electron patterns obtained from pure iron (Fig. 6) were not stable. The point 
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is damaged by a relatively small increase in voltage or an increase in temperature, and 
that is the reason why there could not be observed any vanishing of the emission from 
the edges as in the case of nickel. Schleicher (1952) drew attention to the difficulty in 
obtaining electron patterns from iron. He is the only author who systematically in- 
vestigated the emission from carbonyl iron. He did not, however, obtain good pictures, 
since, as he showed, the electron pattern does not map the entire monocrystal, but 
only a certain small part of it. 

Using a carbonyl iron point with an addition of silicon, we obtained more stable 
patterns. The electron pattern mapped the entire sphere of the monocrystal (Fig. 7). 
Despite the greater stability of the electron pattern, here, too, we were not able to 
make the emission from the edges vanish in order to obtain emission only from the 
corners. 

This effect was obtained only by using a steel point. The changes in the electron 
pattern obtained by gradually raising the voltage are shown in Fig. 8a—c. From Fig 8a 
it can be seen that the emission takes place from both the edges and the corners (the 
monocrystal model is shown in Fig. 9). In Fig. 8b can be seen the partial — and in 
Fig. 8c, the total — vanishing of emission from the edges. With the steel point, we 
did not always obtain the mapping of the entire monocrystal sphere; we succeeded 
sometimes with particular care in preparing the point (lig. 10). 

As in the case of nickel, the phenomenon of “drifting? was evident. Also, as in 
the case of nickel, it was possible to regenerate a damaged point. This, however, was 
now more difficult to do. 

А wider range of photographic material will be published in the Zeszyty Naukowe 
Uniwersytetu Wroclawskiego — Wroclaw University Scientific Reports — (yearbook 
1957—B). 


Discussion 


The above described changes of the electron pattern as well as the observations 
of other authors (Benjamin and Jenkins 1940, Becker 1951, Andreev 1952, Вепуа- 
minovich, Smirnov, and Shuppe 1953, Dyke, Barbour, Martin, and Trolan 1955) 
indicate that with a raised temperature, and in the presence of a strong electric field, 
a monocrystal previously rounded off by heating takes on a polyhedr shape. From the 
changes of the electron patterns, it can also be seen that this shape depends on the 
temperature of the monocrystal and the intensity of the applied electric field, that is, 
with an increase in temperature and electric field intensity, the corners undergo 
a ’sharpening’. As a result, the electric field intensity attains a greater value at the 
corners and edges (smaller radius of curvature) and smaller values at the planes. 
This, in turn, causes the emission to originate practically only from the corners and 
edges. For these very reasons, with a further increase in voltage and temperature, 
the emission from the edges disappears so that only the emission from the corners 


remains. ` í 
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On the basis of the above, it may be supposed that the monocrystal rounding 
by heating in the absence of an electric field is not in a state of stable equilibrium. 
If a strong electric field acts on such a rounded off monocrystal the latter 
takes on a polyhedr shape. This phenomenon was observed when the monocrystal 
temperature was increased. It may be supposed that this also takes place correspond- 
ingly more slowly at room temperature. 

The fact that the central part of pattern forms first, and that the peripheral part 
forms only later, is explained by the fact that the point possesses a non-spherical 
shape (e. g. hyperboloidal) as a result of which there is a strong electric field in the 
central part. 

The mechanism of the transition of the monocrystal from a rounded state into 
a polyhedr shape is not completely explained. The ato mswhich make up local points 
on the corners may arrive there in various ways. The most probable ways is surface 
migration or plastic deformation. : 

The migration of adsorbed atoms on the surface of metal was observed under 
similar conditions by Becker (1951), Sokolskaya (1952), Gorbaty and Shuppe (1955). 
They found that the atoms migrate on the surface of the monocrystal of the point in 
areas where the electric field is the strongest. The explanation of this phenomenon 
involves certains difficulties. Changes of the electron pattern in Figs. 2à—f and 
8a—c may be explained by this kind of surface migration of the nickel and iron atoms. 
After the point has been heated at a high temperature, examination of it under the 
electron-lens microscope does not disclose the existence of corners and edges. 
The observations of. Andreev (1952) indicate that they arise, however, even without 
the participation of external factors, after the previously rounded point remains for 
some time in the apparatus. On this basis, it may be supposed that there exist, on the 
point rounded off at a high temperature, corners and edges of submicroscopic size 
lying below the resolving power of the electron microscope, that is 'embryos', as 
it were, of corners and edges. Haefer (1940) observed such 'embryos' oma tungsten 
point previously heated to a temperature of 2400? K. He established that while the 
raising of the temperature increases the smoothness of the point, it is nearly impossible, 
however, to achieve complete smoothness. It should be supposed that perfect smooth- 
ness is achieved only near the temperature of melting. When the voltage is applied 
and the electric field intensity becomes greater оп the embryos, the surface migration 
of the metal atoms takes place in direction of these embryos. This results in the for- 
mation of corners and edges and the subsequent sharpening of them. The sharpening 
takes place more so on the corners than on the edges, since the electric field is stronger 
on the corners. К 

The above-described mechanism explains also the phenomenon of drifting’. 
Despite the fact that the voltage is raised no further, the patters changes, since the 
migrating atoms sharpen the corners, and by this fact the electric field is further 
increased. ” | | hi 

After the formation of local points on the corners, the emission takes place praéti- 


^ 


# 


The Role of the Emitter Shape 31 


Дм в 


cally only from these points. The current density can achieve very high values, and, 
as a result, produce a local increase of temperature and cause the evaporation of 
the atoms of the metal. The rate of evaporation depends on thé temperature of the 
heating loop, on the Joule heat produced by the emissiori current, and on the heat 
lost as a result of radiation and conduction. Because of surface migration and evapo- 
ration, a certain state of dynamic equilibrium is established, which governs the 
dimensions and shape of these local points. When the speed of the migrating atoms 
is low (small electric field intensity, low temperature) the cerners are not 
yet sharpened. Under these conditions, the radius of curvature of the corners is not 
much smaller than the radius of curvature of the edges; the emission arises both from 
the corners and the edges. Corners lying near one another are mapped in the shape 
of a single spot. When we increase the speed of the migrating atoms by increasing 
the voltage or the temperature, the sharpening action is greater on the corners than 
on the edges. The emission from the corners increases very rapidly while the emission 
from the edges is very small. The pattern now presents sharp contrasts emission now 
being observed only from the corners. As a result of the large degree of sharpening 
each of two corners situated even very close to one another is now mapped separately. 

The formation of the spot in Fig. 5 is explained by the mechanism of surface 
migration, i. e. one of the corners is sharpened to a greater extent than the other. 
Hence practically the entire emission arises from this particular point. Strong heating 
of the entire monocrystal without an electric field will cause it to become rounded 
again. 

The regeneration of a point whose end has melted, or of a point whose end has 
been pulled off, is explained, although with some reservation, by the mechanism of 
surface migration. With a large increase in the speed of the metallic atoms migrating 
on the surface, on of the corners at the end of the point with a large radius of curvature 
begins to build up in the above-described manner. Since this phenomenon takes 
place very rapidly (strong electric field, temperature about 1000* C) the atoms can 
pile up chaotically. The strong heating of the point without the field causes a mono- 
crystal with a small radius of curvature to form from these piled-up atoms. 
Sometimes two or more such monocrystals may arise, and we then observe several 
patterns superimposed on each other. The above-described mechanism of the regene- 
ration of the point gives rise, however, to certain doubts. The pattern from the 
regenerated point is obtained with voltages nearly as high as in the case of undamaged 
point. Since the hypothetical, local piling up of atoms arising through the regeneration 
of the point, should have a smaller radius of curvature, then the pattern should arise 
at lower voltages. 

А second possible mechanizm for reducing the monocrystal to polyhedr shape is 
plastic deformation by a strong electric field. Since practically the entire emission is 
localized in the corners, the current density there is very great. The generation of Joule 
heat causes local temperature to rise to the melting point. The mechanical action of the 


electrical field on such a local liquid surface causes liquid points to form on it by 
£ 
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means of plastic deformation. A similar formation of points on the surface of mercury, 
under the influence of a strong electric field, was observed by Beams (1933), Quarles 
(1935), and Moore (1936). The intensity of the electric field, the surface tension, and 
the rate of evaporation determined the shape and size of the point. The regeneration 
of a melted or pulled-off point may also be explained, using the mechanism of plastic 
deformation, by the creation of one or more local points. During heating in the absence 
of an electric field such local points become monocrystals with small radii of curvature. 

It is difficult to establish which of the above-described mechanisms play the 
more important role, although both doubtlessly exist. 

It should be supposed that at smaller currents, the surface migration mechanism 
predominates, while at larger currents, especially during the regeneration of a damaged 
point, the plastic deformation mechanism predominates. 

The points made from materials that melt at low temperatures are doubtlessly 
covered by a layer of adsorbing gas residue, since they cannot be heated strongly 
enough. However, since the differences in the density of the emitted current caused 
by changes in the point shape are many-times greater then the differences caused 
by changes in the work function as a result of adsorption, this latter factor does not 
play the more important role. 

We wish to thank Professor Jan Nikliborc for his interest in and concern over 
this work. 


KPATKOE СОДЕРЖАНИЕ 


P. Менцлевский и Л. Войда, Роль форты эмитёра при исследовании холодной 
эмиссии. 


Произведены наблюдения над изменениями электронных картин с никелевых 
и железных острий, возникающих под влиянием изменения температуры и на- 
пряжения. Установлено, что эмиссия электронов со слабо накаляемых острий 
исходит из наугольников и краёв монокристаллического острия. 

Вследствие поверхностной миграции, либо пластической деформации, на на- 
угольниках образуются накопления, которые вызывают местные повышения на- 
пряжения электрического поля. В результате разложение густоты тока на элек- 
тронных картинах воспроизводит форму поверхности монокристалла, составляю- 
щего конец острия. 

Замеченное явление дает возможность предполагать, что получение отдель- 
ными авторами разных величин работ выхода из тех же самых стенок моно- 
кристалла было вызвано разною степенью округления острия. 
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POINT SOURCES FOR X-RAYS 


Ву JULIAN AULEYTNER 
Institute of Physics, Polish Academy of Sciences 
(Received July 3, 1956) 


An X-ray tube with adjustable focus of the order of 20 to 70 u was built. The electrod 
beam is collimated by means of an electrostatic lens. The different focus sizes are obtained 
by changing the position of the cathode filament with respect to the grid diaphragms ann 
by changing the grid potential with respect to that of the cathode. The tube wave has an 
original design. 


The purpose of this work was to buil an experimental X-ray tube with a focus 
size adjustable to the order of several tens of microns. 

The fine focus X-ray tube is finding wider application in investigations of structure 
in spectroscopic work and in defectoscopy. 

The first publications on the subject of the design of these tubes and-their use to 
investigate structure appeared in the years 1943 and 1947 (Gunier, Devaux 1943; 
Bezverkhii 1947). To obtain good resolving power, which is especially important for 
the precise determination of the position of interference lines, the shape and structure 
of both the interference lines and the Laue spots, a diaphragm of very small diameter 
or a very narrow slit is usually used. In this way, a considerable amout of the initial 
beam intensity is lost, and, as a result, the exposure time is increased by several times. 
To increase the intensity the effective power of the tube per unit area of the focus 
should be increased. One of the paths leading to this is based on decreasing the dimen- 
sions of the focus by a suitable collimation of the electrons on the anti-cathode. The 
use of point or linear foci permits the first diaphragm or slit to be eliminated from the 
collimation arrangement. 

Point or linear sources of X-rays of high power per unit area of focus are particul- 
arly important for investigations of small-angle scattering, where the exposure, using 
ordinary X-ray tubes and very narrow slits, are frequently of the order of hundreds 
of hours. | 

Small-angle scattering of the Bragg type is characterized by the appearance of 
numerous interference lines in the region of small angles of the order of less than 
20 minutes. In this case it is highly advantageous to use point sources of X-rays in 
order to obtain rapidly interference patterns of high resolution. 
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Fine focus tubes play an important role in investigations of the size of crystallites 
whose dimensions are of the order of from one to several tens of microns (Keller, 
Hirsch, Thorp 1950; Keller, Hirsch 1951, Pines 1955). Samples of such a crystalline 
structure do not yet show a widening of the reflexes caused by the small dimensions 
of the crystallites (as the number of diffracting planes decrease, the angular width 
of the deflected beam increases), however, they already give continuous rings. 
Rings obtained in this case with the help of beams of diameter of several tens of 
microns (35и) show a non-continuous structure, 1. е. they consist of small interferences 
spots. If we know the number of spots on the ring, the main dimensions of the crystal- 
lites may be determined. By investigating the shape of the spots, we can obtain some 
information as regards the stresses in the crystals. The widening of the interference 
rings may be caused by the small dimensions of the crystals as well as by the presence 
of stresses in it (Pines 1955). The use of microbeams makes possible analysis of the 
intensity profile of the lines in order to distinguish these effects. The investigation 
of microstresses in polycrystalline samples requires the irradiation of individual parts 
of the sample, the surfaces being of the order of 0.01 mm?, and the determination 
of the widening and shape of the reflexes. For these purposes as well, fine focus X-ray 
tubes are being used ever more frequently. The use of the tubes permits the investig- 
ation of certain types of microdefects іп monocrystals (Gunier, Temerin 1958, Lamber 


et al. 1953). 


In spectroscopy, the technique of fine beams of electron rays is used to excite 
the emission of X-rays from very small grains of various minerals on the surface of 
suitably prepared polished specimens. With the use of this method the electron beams 


collimated by an electromagnetic lens attain a diameter of 1 micron (Borovski and 
Ша 1955). 


The use of tubes with point foci has a big future both in defectoscopy and in 
medicine. The smaller the dimensions of the focus, the sharper the picture of the 
details of the structure of the irradiated object, this being due to the decrease of half- 
shadows. In order to obtains uniform sharpness in all directions there should be used 
such a focus as would result in a circular effective cross-section, i.e., the projection 
of focus on the plane perpendicular to the direction of observation should be circular. 
The power of the X-ray tube, and therefore, the intensity of the emitted radiation are 
limited by the melting temperature of the anode material. The tube power per unit 
focus area can be increased by decreasing the linear dimensios nof the focus and 
increasing the density of the electrons responsible for the production of X-rays. 
With a decrease in the ratio of the focus area to the cooled area of the anode, the 
conditions for conducting away the heat improve. In a state of balance between the 
amount of heat produced on the flat anode surface and the amount of heat lost as 
a result of radiation and cooling, the permissible tube-load is given by the relation 
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where 1 is the linear dimensions of the focus 
Q is the quantity of heat lost 
О, is the quantity of heat lost by cooling 
О» is the quantity of heat lost by radiation 
Т, is the permissible anode temperature (at the focus) 
Q = Q +0, 
О, = лх(Т,— То) (it is assumed that the focus has a circular shape) (1) 


Q, = ай 


mx(T,—TQli-ral  mx(T, — T 
= ( a EN n 


where @ is a constant, 
ж —is the heat conductivity. 

The smaller linear dimensions of the focus, the greater the permissible load of 
the tube per unit surface area of focus. In ordinary X-ray tubes used for analysis of 
structure, the permissible load is 50 to 200 watts/mm?. 

In fine focus tubes intended for continuous work orders of from a few to several 
hundred kilowatts per mm? are used (Taylor 1949; Thorp 1949, 1950, 1951; Schall 
1950; Drenck and Pepinsky 1951; Bublik and Pines 1952; Cosslet 1953). 

The tube built in our institute fulfilled the role of an experimental station for 
testing various types of electron collimating arrangements. It worked also a source 
of X-rays for investigations of structure. In working out the design, the tube built 
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Fig. 2 


by Ehrenberg and Spear (Ehrenberg and Spear 1951) was taken as the starting point. 
The electron-collimating arrangements used by these authors ron of an electro- 
static lens C, cylindrica) shield E and an arrangement for changing ihe le nal 
distance (Fig. 1). 
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The electron beam is affected by: 

a) the potential of shield E with respect to lens C 

b) the potential of shield E with respect to filament K 

c) the potential of grid C with respect to filament K 

The authors report that the most favourable conditions of operation for the tube 
are obtained when the potential Й = Vp = — 50 volts. A negative grid potential 
has little influence on the focusing of the electrons but primarily limits the emission 
from the side walls of the filament. By a change of distance d of the lens C from the 
anode, different sizes of focus are obtained. The minimum effective linear dimensions 
obtained for the focus by the authors was 40u. This size was determined by photo- 
metry of the line photographs, the lines being the image of the focus. The image was 
obtained by using a narrow slit (Ehrenberg and Spear 1951) between two smooth 
metal cylinders of 4 mm. diameter. The witdth of the lines in this case was measured 
at half the height of the maximum shading. 

In Ehrenberg and Spear's arrangement, the voltage is smoothed by using a con- 
denser and a kenotron connected in parallel with the tube. The tube is stable at a voltage 
of 50 kV and anode current of from 0.3 to 0.5 mA. Allowable maximum load is 


b) 


Fig. 3a. Change of shape and size of focus caused by application of negative grid bias: left to right —OV, 
38V, —6 V, —9 V, —12 V, —15 V, —20 V К 
Fig. 3b. Shadow of braid of thin copper wire obtained by means of X-rays: a) —3 V, b) —6 V, с) —12 d 
Fig. 3c. Photograph of finger; detailed structure of bone can be seen. 
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11 kW/mm?. The anode is cooled with carbon tetrachloride by means of a circulating 
pump. 

| The constructing of the tube at the Institute of Physics of the Polish Academy of 
Sciences can be divided into two stages. In the first stage an electron collimating 
arrangement, the same type used by Ehrenberg and Spear, was built. In the second 
a new collimating system was constructed. The minimum effective linear M ERE 
of the focus obtained in the first case was 45ш, and the shape of the projection of the 
focus on a plane perpendicular to the direction of observation was elliptical. When 
photographs of the focus were made, a strong background originating in the entire 
surface of the anode was obtained, this would indicate strong electron aberration. 
Nevertheless, the tube gave rather good radiograms. Fig. 3a, b, c. Fig. 4, 5, 6. 


Fig. 6 
Fig. 4. Debey-Scherrer rings obtained from very thin layer of PbS dusted on mica. Laue spots from 
monocrystalline mica. Exposure 45 min. Chamber radius 1,8 cm. 


Fig. 5. Small-angle scattering of X-rays by activated carbon 
Fig. 6. Interference patterns with copper wire of 0.1 mm. diam. Chamber diam. 6 mm. Exposure 12 min. 


Тһе electron collimator built in the second stage of the work differs in the following 
details from the arrangement used by Ehrenberg and Spear (Fig. Т}; 

1) Diaphragm D has a different shape; its opening is two and half times smaller. 

2) During operation the distance from the lens to the anode is constant. The 
distance between the cathode filament K and diaphragm D is adjusted under vacuum. 

The external surface of diaphragm has a radius of curvature equal to the distance 
of the diaphragm to the centre of the anode target. There is also considerable difference 
in the constructional details. The tube is completely demountable. The grid is insulated 
from the shield by means of a double ground joint of plexiglass F. Terminals Z facilitate 
the rapid replacement of the filament. The lamp has two windows one of which is 
of mica. This allows direct observation of the focus spot on the anode. The anode, 
which is under high voltage, is cooled with tap water through a 15 m. glass pipe wound 
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Fig. 7. А anode, W anode cooling, О tube window, E tube shield, G rubber seal, F ground plexiglass 
joint, J glass insulator, 52 vacuum bellows for changing cathode distance, M knob with micrometer for 
adjusting cathode to anode distance, K cathode filament, D lens diaphragm, S lens, Z terminals for 


attaching filament. 


in a spiral. Fig. 3 shows the schematic diagram for the electrical equipment energizing 
the tube. 


Operating Characteristic of the Tube 


The lamp works at 50 kV. The anode current depends on: 

1) the potential of the filament with respect to the shield E 

2) the potential of the grid S with respect to the filament 

3) the filament current and anode voltage 

4) the distance between the cathode filament K and the grid diaphragm.D. Fig. 8 
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shows the variation of the anode current as а function of the positive potential of 
filament K with respect to shield F. 

The maximum anode current is obtained with the potential of the shield at about 
— 110 V. After this potential is exceeded, a sudden fall of anode current occurs. If a de- 
finite potential difference between the grid and cathode filament is not maintained, 
the stability of the anode current is very poor. The stability improves considerably 
if a small negative potential (with respect to the cathode) is applied. Fig. 9 shows 
the variation of anode current as a function of the filament current for two different 
grid potentials. At a potential of —9V the stability is considerably better than at a 
potential of —6V. Fig. 10 shows the difference of the curves as a function of the grid 
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potential with different filament to diaphragm distances. The anode current decreases 
with an increase in the negative grid potential and increases as the cathode filament 
is moved towards the diaphragm. 

The application of a negative grid potential of between —1 to —3 volts affects 
the increase in anode current intensity with respect to the measured intensity with 
zero grid potential (in Fig. 11 this increase is denoted by a dotted line). 

Fig. 11 shows the variation of anode current as a function of grid potential for 
one position of the filrment with respect to the diaphragm. In this case, this dependence 
is expressed by the equation 


i= = V? 4- 35V-+ 520 in the range — 17 |V| > — 15 volts 


As the filament to diaphragm distance is decreased, the anode current increases 
(Fig. 12). Beginning from some particular position, this increase becomes linear. 
The curve representing the function i = i(4 1) depends on the grid potential. 

The quality of the lamp is characterized by 
the effective linear dimensions of the focus and 
ЖА its efficiency. It should be mentioned here that 
by effective linear dimensions of the focus’ is und- 
500 erstood the projection of the real linear dimen- 
460 sions of the focus on the plane perpendicular 
420 to the direction of observation. 
The efficiency of the tube is characterized 
50 by the ratio of the power flowing through the 
340 tube to the real area of the focus. 

Figs. 13, 14, and 15 are drawn on the basis 
300 of a series of measurements made for a constant 
260 value of potential of the grid shield апа anode 

voltage and filament current. A change of the 

Í8Ü ^ cathode filament to diaphragm distance within 

14) Фе limits 0.3 mm. causes а visible change of the 

effective linear dimensions of the focus (Fig. 13). 

100 Тһе smallest linear dimensions of the focus is 

- obtained at some definite position of the filament 

Кр Оа БЕП (40.45 scale divisions). As the filament approaches 

` Fig. 11 the diaphragm, the anode current increases (Fig. 

14) and the X-ray radiation attains a maximum 

(Fig. 15). The efficiency of the tube is a function of the anode power and focus area. 

The real shape of the focus is usually elliptical (among other factors, this is due to 

the shape of the filament). Knowing the angle of observation, the ratio of the axes 

of the elliptical-cross section of the beam in the plane perpendicular to the direction of 
observation, we may determine the efficiency of the tube. 
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The efficiency changes witch the distance of the filament from the diaphragm. 
In the case already discussed, where the size of the focus is 24и, the efficienco 
is of the order 2.3 kW/mm?. 

The maximum efficiency achieved with this lamp is 8.30 kW mm? (at a filament 
current of 3.6 A). By changing the filament to diaphragm distance, we can change 
the shape of the ellipse in order that the focus have a circular shape in the plane of 
observation. 

А very convenient method for checking the shape of the focus is to photograph 
a fine mesh. 

In the case of a beam with a circular cross-section, the horizontal and vertical 
wires are uniformly sharp (Fig. 16a, b). The sharpest picture is obtained with a focus 


b) 


Fig. 18 


Fig. 17 


Fig. 16. Picture of mesh obtained with various distances between ud K and diaphragm D. The | 
HE sharpest prone is obtained when the shape of the focus is curcular. Go 
p Fig. 17. Picture of bone З m 
Fig. 18. Radiogram of aluminium wire 
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size of 25u (V, = —6V). А change of V, requires an adjustment of the filament 
position. Fig. 17 shows the structure of bone (У, = —3V). The distance from the 
preparation to the film is 15 cm. 

If we want to obtain Laue spots uniformly sharp in all directions, the position 
of the cathode filament relative to the diaphragm should be adjusted to give the beam 
a circular cross-section. 

In order to obtain sharp Debey-Scherrer rings it is not necessary to use a very 
small focus or a small collimator diaphragm. А good photograph may be obtained 
with a focus size of the order of 60u. 


Fig. 20. Interference pattern obtain ed with aluminium wire. Chamber diameter 36 mm. 


А photograph of a tungsten wire 0.15 mm. thick, obtained by using an ordinary 
collimator of 0.8 mm. diam. in a 57.6 mm. chamber with a focus size of 60u, shows 
` very good line resolution (Fig. 19a, b). Fig. 20 shows the interference pattern obtained 
with an aluminium wire; it is characterized by very sharp lines. The time of exposure 

was 20 minutes. The chamber radius was 1.8 cm. (20' radius 0.9 cm). 

In discussing the results obtained, we should note that by changing the filament- 
to-diaphragm distance and by decreasing the dimensions and shape of the diaphragm 
opening, we can obtain a focus smaller than that obtained by Ehrenberg and Spear. 
The shape of the focus may also be varied in this way. The efficiency of the source, 

‚however, is smaller. The cause of this may be the insufficient heating of the cathode 
` filament and a slight decentring of the cathode tip in relation to the centre of the 
diaphragm (a deformation of the filament may occur during heating). 

Application of a small negative grid bias with respect to the cathode enables us 
to obtain a rather stable source. 

During our measurements the anode current Ec within the limits of 0 to 
5%; with a poor vacuum this fluctuation was from 5 to 10%. E 
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The radiation from the tube is relatively rich as regards a continuous spectrum. 
This allows it to be used for investigation of monocrystals by Laue's method. 

In operating the tube, filament should be centred as accurately as possible and 
the earthed part of the tube should be adequately air cooled. 


KPATKOE СОДЕРЖАНИЕ 
IO. Аулейтнер. Точечный источник рентееновских лучей. 


Сооружена рентгеновская лампа с регулируемой величиной фокуса разряда 
20—70 ц. Пучок электронов подвергунт коллимации при помощи электростатите- 
ской линзы. Разные величины фокуса получаются, как посредством перемещения 
нити катода по отношению к диафрагме сетки, так и посредством изменения 
потенциала сетки по отношению к катоду. Лампа имеет своеобразное конструк- 
ционное решение; 
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ON THE INDUCED REVERSE RESPONSE OF THE OXIDIZED 

ALUMINIUM G—M COUNTER TO VISIBLE LIGHT, THE JOSHI 

EFFECT, AND THE GEIGER COUNTER SENSITIVITY TO THE 
AIR-FLUX 


Ву B. SUJAK 


Institute of Experimental Physics, Bolestaw Bierut University, Wroclaw 
( Received July 7, 1956)! 


А study was made of the induced emission of photoelectrons from aluminium oxide. 
It was found that the aluminium counters with an oxide layer on the cathode, heated to 
500°С for few hours in atmospheric air, showed a counting rate run up effect under vi- 
sible light illumination after a preceding 'excitation'. This effect is ascribed to the generation 
of colour centres in the aluminium oxide during the continuous discharge ‘excitation‘ 
period and to the positive charge (ions) on the surface of the oxide layer. 

After the counting rate run up, the counters show a very high background in darkness 
and a reverse response to light. This reverse response to light is discussed on the basis 
of the Joshi effect mechanism. It is held that the Joshi effect is also (at least parily) due 
to the generation of colour centres in the surface layer of thc dielectric (glass) cathode. 

The spectral sensitivity of the induced reverse response to visible light is given for 
the aluminium counter with the momentary photocathode. This light sensitivity does not 
depend linearly on the intensity of the illumination. à 

Ап explanation is also given of the mechanism of the reverse response of the counters 
to light on the basis of the counter's recovery time. On hand of this mechanism is shown 
that the atmospheric Geiger counter with a negative point may be operated so as to display 


sensitivity to an air-flux. 


Introduction 


It is well known to the physicists working with the not light tight counters, that 
after a period of rapid counting or continuous discharge, a temporary increase takes 
place in the response to light?. This is especially noted when working with the 
aluminium-walled counters. Such an 'excited' counter responds temporarily even 
to red light when operated a short time after a fast counting rate or continuous discharge: 


1 This paper is a conglomeration of four notes sent successively by the author in the 
june 2, 1955; June 18, 1955; September 29, 1955; December, 30, 1955. 
2 The phenomenon was often utilized by many authors as a sensitization process of the photon 
counters in an analogy to the historical sensitization process applied by Elster and Geitel (1911); compare 
also other works dealing the problem of counter sensitization (Kreuchen 1935; Scherb 1948; Grotowski 
1952; Grotowski, Hrynkiewicz, Niewodniczanski 1953; Sujak 1952, 1956a). 

a (49) 
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(Veksler 1949 рр. 293—294). Another phenomenon which till now did not seem to 
be fully explained was the reverse response to light shown by some photon counters, 
as noticed by Duffenbach and Morris (1942). These authors noticed that some of 
the photon counters constructed by them showed in darkness a very high counting 
rate, which decreased when light was admitted and ceased altogether when fully 
illuminated by a light bulb or by daylight. They described this ав а peculiar effect 

of the alkohol-hydrogen tubes' and made 


S no suggestions concerning the mecha- 
5 5 & nism involved the spectral sensitivity ог 
N oec Cid the sensitivity to the change of light 
AO rs atuminium d 4 

< S 3^ — 500°C intensity. 
£ © 2 air Analysis of the two phenomena 
300 | together with the problem of the ’indu- 
| | ced emission of electrons' suggested that 
250 | | the origin Нез in the production of some 
200 | о-1-%%--9-----о-- 9,1240, 5 inthe surface layer ofthe coun- 
| | ter cathode. These centres seemed 
150 ° 1 to have a rather low effective work-func- 
x 1 1 = 
100 j? tion. To begin the study on the cen 
xx tres, the author covered the cathode 
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Fig. l. Dependence of the counting rate run-up 

efect on the time and light wave length for the 

excited counter with an oxide layer on the alumi- 
nium cylinder. 


with alkaline halide layers, in which the 
well-known colour centres may be ge- 
nerated. The first note published on 
this problem (Sujak 1953) suggested 
that the induced emission of photoelec- 
trons as noticed by Kramer (1951) is 
connected with the existence of F-cen- 


tres in the layer on the counter cathode. 

Bohun (1954) also showed in a series of experiments that when colour centres are 
` present in the alkaline halides, an emission of thermo-and photoelectrons into the 
air atmosphere of an opened Geiger counter is observed.? 

In a further study on this problem, in connexion with Kramer's emission of 
'exoelectrons' of the alkaline halides, the author found an unusually high counting 
rate run up effect in some aluminium-walled experimental counters possessing an 
oxide, layer on the inner side of the counter cylinder. The cylinders were specially 
heated in an atmosphere of air at 500°C for 2—3 hours. When, after excitement’, 
the counters were illuminated by visible light, they could run up into a continuous 
discharge-like counting rate at normal operating voltage. This high "counting rate 
‘continued even after the light was removed. Only a decrease of the operating voltage 
to a value near the Geiger threshold could stop this phenomenon. 


E 


3 I am very thankful to Dr. A. Bohun for his kindness in forwarding reprints of his papers. 
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When the counters were triggered by y — rays in full darkness, the run up 
effect was not observed for intensities giving up to 200 counts per second. 

Fig. 1 contains several graphs illustrating the time dependence of the counting 
rate run up effect of the 'excited' counters, the effect being due to the visible light 
of various wave lengths. That the temporary sensitivity of aluminium to the visible 
light is due to the colour centres generated in the oxide layer seems to be accepted 
by the other workers too (Bohun 1955, Grunberg and Wright 1954, 1955, Nassenstein 
1955, Seeger 1955). It therefore seems that the run up effect under illumination is 
due to the colour centres generated in the oxide layer by the gaseous processes in the 
counters (compare — Sujak 1953, 1955b, 1956b) and to a strong electric field generated 
by the positive ions covering the layer of oxide (Granovskii 1952 pp. 196—208). If the 
concentration of the generated colour centres were constant in time, it would be 
expected that the light sensitivity of the ’excited’ counter would also remain constant 
in time (Apker and Taft 1950, 1951, 1952 pp. 246—260; Sujak 1957, in printing). 

In this paper I shall show that there exists a simple method for continuously 
generating the centres and, at the same time, measuring the intensity of the photo- 
electrons released from them by external light. This may lead to the construction of 
a phototube with a "dynamic photocathode’ showing a reverse response to the light 
as noted by Duffenbach and Morris. This method is to utilize the mechanism of the 
Joshi effect (Joshi and Deshmukh 1941). 


The Joshi Effect and the Induced Emission of Photoelectrons 


By the Joshi effect one understands the instantaneous and reversible photo- 
variation (enhancement Ai and diminution — Ai) of the current i of the gas discharge. 
This was mainly observed for low frequency electrode-less discharges. By electrode- 
less discharge it is meant that the electrodes are covered by a dielectric layer. In most 
cases this was a layer of glass in form of a glass Siemens' ozonizer. The effect was 
observed for a number of gases and vapours both of pure elements and compounds 
over a wide range of operating conditions as regards e. g. the applied potential V, 
the intensity J of the irradiating light, the continued exposure to discharge (aging), 
the surface dependence etc. (cf. for example, Mohanty's paper in Z. Phys. 140, 370 
(1955)). ; 

The current intensity і is measured with an indicator in darkness (ip) and in 
light (i,). The net effect is given by 

+ Ai = і — ір 
Joshi's theory of the effect treats the emission of photoelectrons from activated 
adsorption-like ionic and molecular layers on the electrodes (mostly glass) as primary. 
Capture of slow photoelectrons by the electronegative atoms or molecules leads to a 


negative space charge and the negative effect —A i. When the electrons escape capture 
and cause an increase in the current i, the positive effect +4 i is observed. 


59 В. Sujak 


The surface dependence of the Joshi effect was shown by many experiments on 
this problem made by Mohanty and co-workers (Mohanty 1953, 1954, c, d, e; Mohanty, 
Jajaraman and Krishna Rao 1954, Mohanty, Nagarjunan, Srinivasan and Srinivasan: 
1955). Furthermore, the work of these authors demonstrated that an inversion of the 
Joshi effect (from +A i to — A i) can also be observed by increasing only one of two 
parameters: 

1) the illuminating light intensity, 

2) the frequency of the illuminating light, 
when a suitable voltage is applied across the electrodes covered by a dielectric layer — 
soda glass (Mohanty 1954a, b; Mohanty and Srinivasa Rao 1954-55a, 1955c; Mohanty 
and Kamath 1949). 

Mohanty also observed that the positive Joshi effect --Лі changes to the 
negative —Ai under continuous exposure to discharge (aging). This phenomenon 
of inversion is described by Mohanty as due to an increase of the cathode emissivity 
with the discharge (Mohanty 1955—56c, 1955a, b). Recently, I suggested that one 
is obliged to take into account not only the electrode surface-gas interface as the seat 
of the residual electrons deposited by the preceding half-cycle of the voltage, but also 
some of the generated imperfections in the insulating layer (Sujak 1956с). Recent 
work on the induced emissivity of electrons from metal oxides (Bohun 1955с; Bohun, 
Karpiskova and Duskova 1955; Grunberg and Wright 1954, 1955; Kusterer and 
Bruna 1955; Schmid and Lintner 1955; Seeger 1955, 1956; Sujak 1955b), and 
dielectrics (Bohun 1954a, b, c, 1955a, b; Lintner and Schmid 1956; Nassenstein 
1954, 1955; Sujak 1953, 1955a, 1956b; Sujak and Lewowski 1957) confirm the 
suggestion that the increasing emissivity can be partly regarded as due to a progressive 
increase in the concentration of imperfections in the excitet cathode. Mohanty's 
'residual electrons’ seem to form on these imperfections some colour centres like 
Е,Е? — ог M centres, at least in the cases in which the alkaline halides were inside 
the Siemens' ozonizer (Mohanty 1953, 1954d, 1954—55e; Mohanty, Nagarjunan and 
Srinivasans 1955). The colour centres generated during the gas discharge, in the first 
stage, are a source of additional free electrons (the primary of the Joshi effect) when 
photoionized. In this way the colour centres take part in the observed Joshi effects. 
The quantum yield of the emission of photoelectrons from the colour centres generated 
in the surface of the alkaline halide layer depends on their concentration (Apker and 
Taft 1950, 1951, 1952; Sujak in printing). Therefore, when a Siemens' glass tube 
containing finely powdered potassium chloride is used, colour centres are generated 
in the KCl layer spread over the glass walls, and the Joshi effect, along with its inversion 
from positive to negative, is more pronounced. To confirm this point of view, the 
progressive increase of the photoelectric emission from a KCl layer due to the colour 
centres generated by a continuous discharge in a photon counter (compare Sujak 
1953, 1956b) is shown in Fig. 2. - 

The point not clear as yet is why a linear dependence of the Joshi effect on the | 
light frequency is observed. On the basis of my observations of the emission of photo- 


| 
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electrons from the excited KCI layers (compare Fig. 2a) one should expect a selective 
dependence of the Joshi effect on the light fraquency, at least in the cases in which 
KCl powders were used. The same contradiction occurs in the observation of the 


emission of photoelectrons from an abraded aluminium surface (Grunberg and Wright 
1955) and in my observations when the 


Joshi mechanism is evaluated. I shall re- T! 
turn to this point later. У ur J КО 
When clean glass ’electrodes’ are used, [ КОО 
the imperfections also seem to be generated VO (from water solution) 
in the surface layer of or on the glass, as 200 
shown by observations on colour centres іп 280 
glass (Przibram 1953 рр. 73—76; Yokota 260 
1952). 240 
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Fig. 2. Progressive change in the induced emission of photoelectrons from КС] layer, the change being 


due to the progressive excitation by a continuous discharge, in the counter. The KCL layer is on a brass 


cylinder. The exciting discharge was interrupted when photoelectrons b реа under the influence 
of light with a wave length of about 5700 А. 


Fig. 2a. Spectral sensitivity of the excited KCl layer when the emission of photoelectrons к reduced 
to an equal light intensity J. The selectivity decreases with the time of irradiation (compare Sujak 1953, 
_ 1956b). 


In searching for the induced emission of electrons from the glass surface, Sees 
together with Lewowski that such an emission. may be observed, but that the БЕ 
electric field of the Geiger counter is essential. Therefore, the glass probe had to 2 
put inside the counter in order to lay near the sharp point of the counter anode Сал 
and Lewowski 1957). „То show the emission of electrons or photoelectrons, the 
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glass surface had to be abraded or bombarded by a high frequency discharge. This 
emission is much smaller than in the case of alkali halides. 

All this seems to confirm the point of view that in the Joshi effect the colour 
centres on the electrodes must be taken in account. 


Reverse Response to Visible Light of the Aluminium Counter with an Oxide Layer 
on the Cathode 


My experiments with a gridded counter showed that the run up effect of the 
counters with aluminium oxide on the cathode as mentioned in the "Introduction' 
may be strongly influenced by the external electric field. A positive potential of a few 
volts applied to the aluminium cylinder (the grid acts as an earthed counter cathode) 
may stop the run up effect at once, without any decrease in the counter’s operating 
voltage. Furthermore, I observed that a high negative potential applied to the alumi- 
nium cylinder (of the order 10? V) acts to stabilize the large background attained by 
the excited counter when run up with visible light. These facts confirm that the effects 
I observed have their origin in the cathode (oxide layer). 

The gaseous processes accompanying a high pulse sequence in the counter, 
generate colour centres in the oxide layer. A high pulse sequence can be obtained, 
for example, by using a radioactive source or by utilizing the phenomenon of excitation 
of counters. Dynamical stability of the concentration of the colour centres in, and the 
positive charge on, the oxide layer may be reached when a suitable (accelerating 
potential) negative potential is applied to the aluminium cylinder. In my experiments 
I utilized a process which may be described in short as follows. To begin with, the 
counter was excited’ by a continuous discharge of short duration. When the voltage 
was decreased to the normal value on the working plateau, the background was larger, 
but decreased with time. When a high intensity of illumination from a light bulb 
was admitted, the run up effect began. After one or two minutes of such irradiation 
the counter was run up and showed an enormous background in darkness, of the 
order of 10* counts per second. 

Such a counter tube then showed a reverse response to the entire region of visible 
light. The very high counting rate in darkness decreased when light was admitted 
(observed with a rate meter) and ceased altogether when fully illuminated by a light 


bulb. If the counting rate measured with a counting rate meter is denoted by N 5 
2 ый t 


the background of the run up counter is denoted by (7) and the counting rate 
ar 
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registered by the rate meter when light is admitted is denoted by (7) Tid the dif- 
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may be regarded as the measure of the response to light by-way of analogy to the 
measure of the net value of the Joshi effect. This is evident from Fig. 3. 

In Fig. 4. a spectral characteristic of the response to the source of light used is 
shown (curve A). The spectral response to light, when reduced to an equal value of 
light intensity, is given by curve C. 

It is not clear to what extent the real induced spectral sensitivity to light of the 
oxidized aluminium is influenced by the gaseous processes in the counter. In the 
region of the light spectrum examined no selectivity was observed, contrary to the 
recent observations of Grunberg and Wright (1954, 1955). This seems analogous to 
the fact that no selective dependence 


on the light frequency is observed in 2 N 474702 
the experiments on the Joshi effect, #3 И Ра а? a 
: J 
аз already mentioned. Furthermore, i ров : Ait 
it is interesting to note that Schmid Ж ` aluminium / РЕ 
о é К 
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Fig. 3. Reverse response to light. When light is admitted, the counting rate registered by the rate meter 
decreases. The picture of the oscilloscope screen as it may be observed is given schematieally. 
Fig. 4. Curve A (—x—) depicts the response to light of the oxidized aluminium counter, when working 

in the reverse manner: 


A = SRON 


Curve B (—.—) depicts the light intensity with respect to the wave lenth: J = F(A). The intensity was 
measured in arbitrary units with a thermopile and a galvanometer. The monochromatisation of the light 
is given at the top. | Е 
Curve С (—o—) shows the response to light when normalized to equal light intensity of the illuminating 
light, regardless of the fact that the dependence is not linear: 


N 


Е Ф (1). The response to light 4 х is given in the scale of the counting rate meter without any 
J t 


correction for dead time 
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aluminium, but in the near UV region. The results of Schmid and Lintner are 
reproduced for comparison in Fig. 5. 

One must take into account the rather strong electric field due to the counter's 
anode and the adsorbed gaseous ions, when the oxidized aluminium counter's reverse 
response to light is observed. The 
influence of the electric field on the 
induced photoemission from abraded 


aluminium is already evident from 


(erue the experiments of Grunberg and 
%0 Wright.5 The electric field influences 
by E Schmid not only the intensity of emission, but 

К Linlner 


also the spectral characteristic. This 
agrees with the suggestion that the 
counter’s field strongly influences 
the observed, induced emission of 
3400 3800 330 Ф) А electrons (Grunberg and Wright 1955, 


Sujak and Lewowski 1957). This 
Fig. 5. Emission of photoelectrons of extended alu- "ee 
minium as a function of the wave length of the light 


(reproduced, from E. Schmid and К. Lintner — com- 


agrees with Bohun's sug- 
gestion that the mechanism of the 


pare Anz d. Österr. Akad. d. Wiss. math-nat. Kl., 164, Malter effect should be revised | ) 


158 (1955). (Bohun 1954с). 


Mechanism of the Counter’s Reverse Response to Light 


In an analogy to the change in the amplitude of the high frequenty part of the net 
values of the negative Joshi effect (Mohanty 1950), my observations on the oscilloscope 
screen showed that when light is admitted, the amplitude of the counter pulses decrea- 
ses. This is registered by an electronic device (the pulse rate meter) as a reverse response 
to light by the examined counter. The dependence of the counter's response to the 
light intensity, when operated in the described manner was not linear (Fig. 6.) This 
is fully analoguos to the Joshi effect (Mohanty and Srinivasa Rao 1954— 55a). 

The diminution of the amplitudes of the counter pulses is ascribed to the "recovery 
time' of the counter and the RC constant of the electronic device. Trost (1941) already 
showed that the electronic device registering the pulses of a counter shows a dimi- 
nution of the registered counting rate when one increases the already large intensity 
of the particles then triggering the counter. When one reproduces the curves given 


* Tam very thankful to Professor E. Schmid and Professor, K. Lintner for their kindness in 
forwarding reprints of their papers. 


5 T also thank Dr. К. Grunberg and Dr. K.H.R. an for kindness in sending their papers on the 
'exoelectrons' problem. 
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by Trost (Fig. 7.) the reverse response 
to light is immediately clear. When 
the counter's background reaches the 


х 
value of M ‚ фе admitted light 
crit 


will be registered as a decrease in the 
counting rate. 

In the cases observed by Duf- 
fenbach and Morris (1942), as well 
as in those described in this paper, 
the additional triggering particles 
arethe photoelectrons released by the 
light from the continuously ’excited’ 
walls of the counters. The reverse 
response of the counters with alumi- 
nium oxide on the cathode, even to 
the red light, must be considered in 
connexion with the entire electronic 
equipment, in analogy to the Joshi 
effect (Mohanty 1950), and especially 
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Fig. 6. Response to the light intensity of the gridded 
counter used in the experiments. The light intensity was 


changed with the aid of metallic grids 


the amplifier and pulse discriminating unit. This may lead to a simple changing 
of the sensitivity to the light intensity. This manner of operating a light counter as 
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Fig. 7. Registered counting rate as a function of the counter discharge frequency (reproduced from 


4 А. Trost — compare 2. Phys., 117, 257 (1941)). 
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а phototube with a simple 'momentary cathode’ may be especially useful when the 
counter is filled with a non decomposing mixture of a low operating voltage. 
The used experimental counters were of the same dimensions as those already 


described — Sujak (1955b). The grid had 400 meshes/cm? and was about 1,5 mm. 
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Fig. 8. а) Тһе Siemens' ozonizer, as used by Mohanty and others in the study of the Joshi effect (repro- 
duced from S. В. Mohanty — compare J. d. Chim. Physique, 52, 816 (1955)). 

b) Diagram of the gridded aluminium counter with an oxide on the inner surface of the counter cylinder, 
as used by the author (compare Sujak, Acta phys. Polon., 14, 489 (1955)). 


from the aluminium cylinder. The counters were filled with the standard mixture 
of alcohol and argon or air and alcohol. 

The schematic view of the Siemens' ozonizer used by Mohanty and other workers 
in the study of the Joshi effect and of the gridded counter used by myself are shown 
in Fig. 8. for purpose of comparison. One may see that the difference lies only in the 
principium of the current measurement and in the materials of the electrodes. 


The Air-flux Sensitivity of the Atmospheric Geiger. Counter 


Mohanty and co-workers (Mohanty, Nanda and S. C. Mohanty 1954) also studied 
the Joshi effect in stationary and flowing oxygen and observed that the negative effect 
—Ai is smaller for flowing than for stationary oxygen. They also showed that it decreases 
progressively with the drift velocity. These workers described this as: 'due to the viscous 
drag of the flowing gas, wall-adsorption is less in flowing than in stationary oxygen. 
The emisssion of photoeléctrons is consequently less. Further, the negative space 
charge, which is immediately responsible for — Ai is attenuated. These effects are 
more marked, the greater the drift velocity’. The oscillograms obtained also show that | 


the drift velocity influences the high frequency рад of the discharge in oxygen at 
a pressure of 75 mm. Hg. 


i 


On the Response of the G-M Countesr to Light 59 


Another worker (Werner 1950) showed earlier that the glow current from а point 
of fine Pt wire in air is slightly sensitive to the velocity of the air. He measured a small 
decrease in the glow current when air flowed toward the Pt wire and suggested that 
this is due to the cooling effect of the discharge space near the Pt wire. 

In an analogy to the light sensitivity of the oxidized aluminium counter, the idea 
arose to investigate this problem with an open atmospheric Geiger counter operated 
in the continuous discharge region (Sujak 1956 4). 

The negative point-to-plane corona, glowing in air at atmospheric pressure, may 
be looked upon as the result of a large number of selfquenched electron avalanches 
per unit time (English 1948, Kapcow 1953). The selfquenching effect occurs only 
when an electronegative gas component is present. It is ascribed to a cloud of negative 
ions strongly influencing the electric field in the vicinity of the sharp end of the cathode. 
When a sufficiently high potential difference exists betweer the cathode and a flat 
anode, the current pulses follow each other at a high frequency and the negative ions 
greatly decrease the amplitude of the voltage pulses, occurring across the resistance R 
placed in series with the point cathode (Kapcow 1953 p. 381; Loeb 1950 pp. 645— 
653). This also describes the mechanism of an Geiger counter with a negative point, 
working in atmospheric air in the continuous discharge region. When the applied 
voltage is chosen so, that the glow current in the 
Geiger counter is I, (compare Trost’s curves іп 7 


~8/msec 174 тес 
Fig. 7.) one should expect a decrease in the current 10 
under the air-flux, as shown by Werner. But this 77 
means that the pulse amplitudes must increase in 9 
order to be registered as can be seen from the 4) 
curves given by Trost (1941). 7) 
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Fig. 9. Diagram of the open Geiger counter working in the atmospheric air. The discriminating unit 
includes an pre-amplifier. The counter shows an airflux sensitivity, when operated in the continuous 
discharge region. 

Fig. 10. Graphs showing in arbitrary units the counting rate of the registered pulses as a function of the 


Я 0 $ 
air velocity. The air velocity in terms of pressure differences: Др = E: v? as measured with a Prandlt 


tube. Each curve was measured with the amplitude discriminator set to a different discriminating factor. 
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sisted of a sharp cathode and a metal face-ring. № was used together with a high 
voltage stabilizer, a cathode follower, a discriminating unit and a counting rate 
meter (Fig. 9.). The influence of the air-flux on the pulse amplitudes could also be 
observed on the oscilloscope screen. When an appropriate voltage V (V, < V < V3) 
is applied to the counter, the light phenomenon of the negative point corona can be 
observed at the sharp end of the cathode. The amplitude of the observed pulses is 
very small. When these pulses pass through a discriminator the counting rate meter 
shows no counting rate at all. But if an air stream is directed into the discharge space, 
the amplitude of some of the pulses increase. The counting rate meter now registers 
the increased pulses. Fig. 10. shows a graph of the registered counting rate due to 
the air velocity. Each of the curves is measured with the pulse discriminator set for 
a different discriminating factor, but for a constant voltage applied to the counter. 
The oscilloscope observations show, furthermore, that there exist critical voltages 
similar to those described by Mohanty for studying the Joshi effect. 

The following approximation describes the change in the amplitude of the voltage 
pulses as a result of the negative ion concentration or pulse frequency: 


I Т 
ЛЕТТІ Т езу 


A X: AA 


for М, < N « №, and V, < Vy. 


Here A is the amplitude of the voltage pulses, N their frequency, V, and V, are the 
values of the voltage, N,, N, are the frequencies of the pulses at the critical voltages V}, 
V,. К denoted the negative ion concentration in the vicinity of the cathode when the 
frequency of the pulses is N. 

The values of N or K may be changed by: 

a) blowing away a number of negative ions per unit time; this results ina ‚ decrease 
(—AK) of K, 

b) cooling the vicinity of the sharp end of the cathode; this results in a decrease 
(-АУ) of N the number of primary electrons initiating the avalanches. (The same 
may be due to a local and momentary increasing of the air реши at the sharp end 
of the cathode). 

This suggests in general, that a change of N or K may be measured by means of 
the change in pulse amplitude. In our case this was done by discrimination of those 
pulses whose amplitudes were smaller than a certain fixed value. 

As may be seen, the G—M or the Geiger counter may be used to measure non- 
radioactive events (compare also Anderson and Hertz 1955) when one utilizes the 
changes of the gaseous discharge parameters only.$ 


6 Added in proof: Some results of my experiments dealing this problem, and showing, 
that a Geiger counter may be used for studies of surface charges generated on insulators will 


be published in the Zeitschrift fur angewandte Physik. 
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KPATKOE СОДЕРЖАНИЕ 


Б. Суяк. О стимулированной обратной чуветвительности счётчика G— М на видимый 
свет, о явлении Joshi и чувствительности счётчика Гейгера на струю воздуха. 


Исследована стимулированная фотоэлектронная эмиссия с окиси алюминия. 
Установлено, что алюминиевые счётчики, нагреваемые в атмосферном воздухе 
при температуре 500° С в течение нескольких часов, проявляют состояние разгона 
под влиянием видимого света, если предварительно подвергнуто их процессу „сти- 
мулирования”. Выступление этого явления приписывается образуемым в окиси 
алюминия цветным центрам, а также положительному заряду (ионы) устанавли- 
вающемуся на поверхности слоя окиси, что происходит во время ,,стимулирова- 
ния” счётчика непрерывным разряжением. 

После разгона, счётчики обнаруживают в темноте очень высокий фон и обрат- 
ную чувствительность на свет, состоящую в том, что под влиянием света, система 
регистрирует меныпую частоту импульсов, нежели в темноте. Эта обратная чув- 
ствительность на свет разбирается на основе механизма Joshi. Считается, что вы- 
ступление являения Joshi следует тоже (по меньшей мере частично) приписывать 
образованию в поверхностной прослойке диэлектрического катода (стекле) цвет- 
ных центров. | 

Измерено спектральное разложение обратной чувствительности алюминевого 
счётчика с динамическим фотокатодом. Чувствительность эта не зависить линейно 
от силы падающего света. 7 

Механизм описанной обратной чувствительности счётчика на свет объясняет- 
ся тоже при помощи его времени восстановления (мертвого времени). Исходя из 
этого механизма показано, что атмосферический счётчик Гейгера с отрицательным 
острием можно употребить так, чтобы он проявил чувствительность на струю 
воздуха. 
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DECAY OF PHOSPHORESCENCE OF RIGID SOLUTIONS 


Ву Mieczystaw FRĄCKOWIAK 
Photoluminescence Laboratory, Institute of Physics, Polish Academy of Sciences, Тогай 
(Received July 23. 1956) 


Ап investigation was made of the decay curves of the intensity of phosphorescence 
components parallel and perpendicular to the direction of vibration of the electric vector 
of the linearly polarized light exciting gelatin phosphors dyed either with acridine 
orange or with acridine yellow. The decay curves appeared to be, in agreement with 
Jablofiski's theory, a superposition of exponential functions with different decay 
constants. Different lifetimes correspond to different groups of luminescence centres 
with different limiting degrees of polarization. Changes of the degree of polarization of 
the total phosphorescence during the further stages of the decay (longer than 10-8 
sec. after the interruption of the excitation) are conditioned by the emission of light 
from successive groups of centres playing ап essential part in the phosphorescence. 


Introduction 


The investigation of the changes of the intensity of luminescence during 
the emission of light from a luminophor after the interruption of the excitation 
permits the most direct determination of the probability. of the process leading to 
luminescence. If the emitted light is polarized, then it is useful to investigate se- 
parately the components of the luminescence intensity parallel and perpendicular 
(denoted by I(t) and I, (t), respectively) to the direction of vibration of the electric 
vector of the linearly polarized excitating light. This type of investigation also permits 
the determination of those processes which-lead to changes of the degree of polari- 
zation of the emitted light. Such a process may be the migration of the energy of exci- 
tation from the initially excited molecule to the initially non-excited molecules, 
the orientations of which may be assumed to be d'stributed isotropically. 

Although any fluorescent process is simpler than a phosphorescent one in accor- 
dance with the level diagram of Jablonski (1935), the investigation, however, of the 
fluorescence decay curves presents much greater difficulty. This is connected, on the one 
hand, with the necessity of investigation a process of such short duration (~ 1058 to 
10-9 sec), and, on the other, with the production of very short (atleast 10-19 sec), uniform, 
and reproducible impulses of excitation radiation. In the already mentioned case of 
phosphorescence, the requirements are less acute and it was possible to use for 


_ (63) 
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this type of investigation a mechanical phosphoroscope of special design. In general 
any experimental arrangement used for this investigation should have the following 
properties: 

1) It should have a large time base in relation to the duration of the luminescence 
in order that the phenomenon be spread out in time аз much as possible and that 
large resolving power be obtained. | 

2) It should maintain a constant intensity of excitation light (or uniformly shaped 
pulses of uniform intensity). 

3) It should not have any depolarizing element between the luminophor and 
the light receiver. i 

4) The geometry of the system (source of excitation light, polarizer, luminophor, 
analyser, receiver) should be invariant in time in order to-have the certainty that the 
observed changes of light intensity are connected only with processes arising in the 
luminophor. 

5) The time of irradiation by the excitation should be small enough in relation 
to the duration of the luminescence. 


Experimental 


The work described was limited to the investigation of the phosphorescence 
decay curves of organic dyes with a gelatin base. The diagram of the phosphoroscope . 
designed for these investigations is shown in Fig. i. 

The light source is a 1000 W/110 V bulb operating on direct current. The light from 
the bulb, after passing through water filter W and polarizer Р, (transmitting vibrations 
parallel to the length of slit S,), is concentrated on slit S,. The time of excitation depends 
on the width of the slit S, and the linear velocity of the luminophor behind the slit. 
The phosphor, in the shape of an 'driving belt’, moves between slits S, and Sg, which 
shield it from outside light for a space of 12 cm. The luminescence decay is investigated 
in this space. Slit S, has a special slot which separates the space in which the phospho- 
rescence is observed from the nearest edge of slit S, by a lining 5 mm. wide. The light 
passing through slit S, and the phosphor falls into a special trap K which protects 
the phosphor from scattered light. The phosphor is set in motion by a motor 5 whose 
speed of rotation is controlled by voltmeter V, and tachometer T. The phosphor is 
supported on rollers R on edges of its surface only. The phosphorescent light is photo- 
graphed by means of а ’Contax’ camera on ’Ultrarapid’ film. The objective of the 
camera is equipped with a special attachment which permits the diaphragm of the 
objective D to be set accurately. The diaphragm is calibrated so that a definite position 
corresponds to a definite and known reduction of light intensity falling on the 
photographic film. The calibration is made by means of a rotating sector disk with 
openings cut out in the ratios 1:2:4:8:16:32:64 and with the light colour (from a double 
Hilger-Müller monochromator) as close as possible to the investigated phosphorescent 
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light. In this manner Jablonski’s method (Jablonski 1953) сап be used to determine the 
sensitivity curve of the film without resorting to any external attenuators, which in these 
investigations are harmful because of the necessity of making the investigations under 
uniform conditions of excitation and without introducing additional depolarizing 
elements. The above-mentioned method is based on the photographing of the same 
phosphorescence band twice, but with different settings of diaphragm D and (giving 
a known attenuation of the light beam falling on the photographic film). In this manner 
we obtain two mikrophotometer curves, one of which is attenuated. By finding the points 
corresponding to equal density on both curves, we obtain the sensitivity curve. 

The camera is connected with the phosphoroscope in such a manner as to permit 
the changing of the diaphragm setting and the film without changing the setting of 
the camera. To determine the magnification of the image, a transparent centimetre 
scale is mounted at the phosphor between slits S, and Sg. It is illuminated through 
a special window at the front of the plate of the phosphoroscope and photographed. 
The distance between the lines of the image of the scale is determined very accurately 
by a comparator. At the same time the image of the scale serves to control the correct 
setting of the camera and the sharpness. Ап analyzer P,, which can be set in two po- 
sitions, is placed in front of the objective. The appropriate setting of P, with 
respect to P, permits 11(4) ог Г, (#) to be observed. 

’Ultrarapid’ films were used to make the phosphors. Unexposed rolls of film, 
about 2 m long, were fixed and dried. They were then placed in a bath of a solution 
of the given dye in water for an hour and then rinsed for 10 min. in running water. 
The phosphors were dried and kept in total darkness. The measurements for each 
phosphor were made when the same amount of time had passed as counted from the 
time of removal from the bath. This time was 4 days. Phosphors of different concen- 
trations were obtained by using baths of different dye concentrations. ’Ultrarapid’ films 
showed birefringence. Although birefringence of the base of the.film does not play any 
role, since the gelatin phosphor lies on the side of the camera, the induced birefringence 
of the gelatin may however be harmful. To investigate this, the gelatin was washed 
away over an area of an unexposed and fixed film. The border between this are: 
and the gelatin was marked by a notch. After inserting at the place of the phosphor a roll 
of film so prepared and introducing the notch which separated both parts of the roll 
between two crossed Nicol prisms, not the least difference in the brightness of both 
parts was visible (the slit was clearly visible). The proper operation of slits 5, and Sg was 
controlled in the following manner. An unexposed and-fixed film roll was placed at 
the phosphor. The roll was put in motion and after it had attained the speed used n 
the measurement and with the same illumination of slit S}, it was photographed, the 
time of exposure being twice as long. If no blackening was noted on the photographic 
film, it was assumed that the slit was operating properly. Such a check was made before 
and after each series of measurements. i i Ў 
The photometric measurements were made with a Zeiss microphotometer. The 
concentrations of the dyes were determined less accurately by means of a visual Zeiss 
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colourimeter with the assumption that the absorption power of the given dyes in water 
and gelatin solutions is the same. The light was made monochromatic by means 


of a Hilger-Müller monochromator. 


Results of Measuerements 


The measurements were made on gelatin phosphors dyed either with acridine 
yellow or with acridine orange. They were made in the following manner. After the 
scale was photographed and the proper functioning of slits 5; and 5. checked, the 


Fig. 1. Diagram of Phosphoroscope. 


V4, Г, — Voltmeters 
Ly, L4, L4 — Lenses 
P,, Р, — Polaroids 
51, Se, Sg — Slits 

K — Trap 

Я — Water filter 


R — Rollers 

F — Phosphor 

D — Diaphragm 

B — Photographic film 
T — Tachometer 

Z — Lamp 


analyzer P, and diaphragm D were 
set in the required positions. The 
motor was turned on and after 
it attained the proper speed, the 
phosphorescent light was photo- 
graphed. The time of exposure was 
measured with stop-watch to an 
accuracy of 0.1 sec. During the 
measurements the voltmeter and 
tachometer readings were checked. 
Next, a measurement was made 
with a changed setting of P, this 
procedure being repeated. In this 
manner it was attempted to elimi- 
nate accidental errors in setting the 
position on analyzer P,. Each of 
the components /j(t) and I (t) | 
were photographed three times. 
To determine more accurately the 
sensitivity of the photographic 


film, ae I, (t) ит was photographed three additional times with а changed 

setting of diaphragm D. To obtain the same conditions during the developing and 

fixing of the film, all photographs for the same dye were made on a single film. 
The conditions during the measurements were the following: 


а) the measurements were carried out in room temperature; 


b) the speed of the phosphor was 1476 + 60 cm/sec; 


c) the time of excitation was 1.35 - 10-4 sec. 


. d) the decay curve was observed from the instant 3.38 - 10-4 sec after interruption 
of the excitation to 5: 10-3 sec; 
e) time of exposure was 20 sec; 
f) the excitation light was white and linearly polarized. 
. As already mentioned, the photometry of the photographs was performed by 
means of a Zeiss microphotometer with a slit width of 0.2 mm.; the image of the film 
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on the slit was magnified 24 times. The galvanometer deflection was noted every 0.2 mm. 
of the image of the phosphorescence band; this corresponded to a decay time interval 
of 7.08 10-5 sec. The "beginnings! of the blackening were determined with particular 
accuracy, since according to the predictions of Jablonski’s theory (1935, 1936, 1955), 
a mutual shift of the maxima of both components of the luminescent light intensity 
should be expected. However, in these experiments, such shifts were not noted. The 
. mean value of the photographic density for three measurements was taken. Then the 
decay curve was decomposed into exponential curves. The curves were drawn on sheets 
that were large enough to permit the points to be located with an accuracy of 0. 005 on 
natural logarithmic scale. In this manner it was possible to represent the decay curves 
аз а superposition of the exponential functions J, (t) = C,e-"* . The error of measure- 
ment of 1 (1) and Г, (t) did not exceed 9%. The increase of the temperature of the 
phosphors at the places of excitation was estimated as not exceeding 5 · 10-5°С. 

Fig. 2 shows In/, (2), In! | (t) and the corresponding exponential functions, obtained 
from the analysis. The degree of polarization of the phosphorescent light P(t) as а 
function of time for acridine yellow in gelatin is shown for three different dye intensi- 
ties: C, = 51.7 g/l, C, = 67.8 g/l, and C, = 99.4 g/l. The degree of polarization P(t) was 
determined directly from the value of the intensity curves I(t) and Г, (t) asa function 
of t. As may be seen from the curves, in all cases, the decay curves I(t) and Г, (t) 
may be represented in the form of a sum of a few exponential functions. We may 
therefore write the following formulae: i 


In (t) = > a, e * = » Lie (0 (1) 
k=1 k=1 

0 = Ў bpe = n 0) 
k=ł k=1 , 


It is worth-while emphasizing that in equations (1) and (2) the same decay constant 
а, appears. An exception is curve Сз (Fig. 2) in which the decay constants of the first 
functions differ considerably. However, this exceptional departure may be explained 
by the. fact that they represent some function which is the resultant of two or more 
exponential functions whose powers differ so little that with the phosphoroscope used 
for measurements it has not been possible to separate them. We may therefore assume 
that in this case as well; the decay curves may be described by equations (1) and (2). 
We assume (see below) that the phosphorescent light described by the functions 
а, е7% and b, е-% is emitted by a group of centres (k group) with similar properties 
(like molecules which neighbour the excited molecule have the same confi- 
guration). 

‘The total degre of polarization of the luminescence emitted by all groups of centres 


is occording to the definition: 


Ji — 1.0 E К 


ЕТО 
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Fig. 2. Acridine yellow. Natural logarithm of the intensities versus the total degree of polarization 
--ө-- Parallel components and their exponential functions. 
--О-- Perpendicular components and their exponential functions. 
—x— Total degree of polarization. 
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We define the limiting degree of polarization for the K group of centres as: 


wy () — 1,600 


Auer Isa (t) + 14, (t) кое e 
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Table I 
Acridine Yellow 
Concentration 
g/l Component T sec P, (оо) % 
2 : 
Ji 220 LOSS 
| 1 
Zt 2.29 x 1073 
= |9244 
Jj 3.20 x 1074 | 
| P, = 382 
Zu 3.20 x 10-4 
I 2.27 x 1073 
| 210 
Tif 2.27 x 107? 
2 m 
29 102% 
С, — 678 ^i EM i P, = 29.0 
Ji БОХ 10-- 
Л М АИ s 
1 P, = 432 
8.72 x 10-5 


1.585 10-2 


Р, 152 
ilt se MO 
531 1025 
Р ж 197 
5.37 x 1074 
С; == 99.4 
2.96 xX 10-4 
a P, = 34.9 


` 2.36 x 1074 


3.54 x 105 
4.00 x 1075 


The decay curve of the total luminescent light intesity, can be described, because 


of the relation 


I, () + 21,6) = 10. : © 
by the equation 
Fit) = >, асы 75 (0; (6) 
іті k=1 
where 
ere dı чн д 


Particular values of the decay constants, as well as the limiting degree of polari- 
zation, determined directly from the difference of the logarithms of the intensities 


Ije (t) and J, , (0) ps listed in Table I. 
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Fig. 3. Acridine orange. Natural logarithm of the intensities and total degree of polarization versus time. 
—e@— Parallel components and their exponential functions 
—O— Perpendicular components and their exponential functions. 


---x--- Total degree of polarization. 
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Table II 


Acridine Orange 


71 


С ЗЕ е 
а-я ion Сорос ee P,( 00) % 
2 -4 
ЛІ 8.85 х 10 P, = 28.0 
n 8.85 x 107% 
Су - 50 priv 1 
-4 
Jii 1.10 x 10 Pie 120 
"um 1.01 x 1074 
3 -3 
Ji 1.54 x 10 Po= 210 
N 1.54 x 10-3 
2 2.36 x 107% 
NS Jii IDE РЁ 
2 74 2.36 x 107* 
1 7: 1075 
Jj ays P, = 28.0 
ТАТ 7.09 1055 
3 -3 
Jj Tent ЩЙ PP ToT 
A ПОШ < 1O 
2 —4 
Jj 2.93 x 10 Р, = 21,6 
Co ==12.6 2 2.80 10-* 
218 
1 m 
Jj 6.45 x 10 P, = 353 
T 6.25 x 1075 
3 
Jii ps = 15.6 
Ji 
- 3 
Ji Ба = 21.8 
C, = 16.3 iG 
1 -5 
Jj 6.76 x 10 320 
b» 6.76 x 10-5 


To check whether equations (1), (2), and (6) have a more general character, a fully 


analogous experiment was made in which acridine orange was used. The investigation 


was made for four different dye concentrations C, = 5.0 g/l, C; = 7.2 gll, Сұ = 12.6 g/l, 
6.16.3 g/l. | 
І The results are presented іп Fig. З and in Table П. As may be seen from Fig. З . 


and Table II the decay curve, in this case too, can be represented by equations having 
F 
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the same form аз (1), (2), and (3). It thus follows that the equations have a general 
character and are not closely connected with the type of dye, but rather with the type 
of processes arising in the luminophor during the emission of light. 


Interpretation of the Results 


From among the various theories of the decay photoluminscence (Feofilov and 
Vavilov 1942; Frank and Vavilov 1931; Jablonski 1935, 1936, 1954, 1955; Shishlovskii 
and Vavilov 1934; Vavilov 1942, 1943, 1944, 1950), Jablonski’s theory best agrees 
with the results of the present investigation. 

We take into account the fact that the intensity of the light emitted by the k-th 
group of centres is proportional to the probability of a transition with phosphorescent 


emission y, and to the probability of finding a given (k-th) centre at the instant the 
5-24 


y 
excitation is interrupted (initial moment) (1 — ei) e? N . Then, ma- 


king the assumption that y, == a, (see below), the theory for the simplified model 
of the luminescence centre gives the following expressions for the decay curves: 


һ () Я [Aue 2, еее = DM: G 
до = > [Dj +: — B, e (Fee) | = xu 11.0, (8) 
k=1 БЕЗІ 


ғу 


M сз Жаға рес. 453 xc E M 
% д von 2410 DL SER) we, = 9) 
where BE А 2D 
апа 
on АГА 2(1 — Qo) zs 

i (ее Е а EE 

er жа (1—0) (E — 1 
ар 


er суйем ln 
Dee A E лырга от 
JT NE) ELE ( er aane), 


хь (1 — ета“), 


where v = n'v, 

v — volume of the 'active sphere', 

n — number of molecules causing luminescence per cm? of solution, 

0, — degree of polarization, 

&,— sum of the probability of internal quenching, probability of quenching by foreign 
molecules and by like molecules, and probability of a transition y, with 
photoluminescent emission for the centre belonging to the k group, 

и — probability of migration of the excitation energy, 


— қаларлық чуч» 
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t — time, 
“ару Raber" 
a, (1 — е7) apr — probability of findling a centre belonging to the k group 
at the instant 2 = 0 (Smoluchowski — Poisson distribution !)- 


The limiting degree of polarization for the k group of centres may be expressed 
by the formula: 


Lec ызга с 
k(3 — р) р. 


where p, denotes the so-called ground polarization. 

The similarity between equations (6) and (9) is striking. From a comparison of 
formulae (4) and (6) with theoretical formulae (7) and (8) we see that the first pair. 
of formulae lack the expressions corresponding to the second term in formulae (7) 
and (8). It should be noted that these terms contain the probability of migration of 
the excitation energy, which is responsible for the depolarization of the luminescent 
light. However, equations (4) and (5) represent the same shape as (7) and (8) if we 
assume either that ku + x, 2» а}, which means that the probability of migration of 
the excitation energy is very much larger than the total probability of a transition with 


(10) 


emission of phosphorescence, or that u = 0, which is equivalent to the assumption 
that the migration of energy does not occur. The second alternative would, however, 
be contrary to the present experiment as well as contrary to work by other authors 
(e.g. Bowen and Mikiewicz 1947). As may be seen from Figs. 4 and 5 and Tables I and 
П the partial degrees of polarization p, are,smaller than the ground polarization p, 
(we may assume the latter to be of the order 45%). 

There thus existsa factor causing depolarization of luminescent light, the probability 
of which is different from zero. Since in the phosphors used in the investigations, 
the migration of energy from the primary excited molecules to the primary unexcited 
molecules may be the chief cause of the depolarization of the luminescent light (see 
below), we must assume the first alternative. In the conditions under which the 
experiments were conducted it is sufficient to assume that ku + a, 108 secl. 
However, in the experiments described here, it was impossible to determine the pro- 
bability of migration of energy because of the rather large amount Of time elapsed 
from the beginning of the excitation to the moment of observation. 

The agreement of the theoretical expressions with the experimental results allows 
us to accept the model of luminercence centre for which these expressions were 
derived. By (simplified) model of a luminescence centre we shall understand 
the excited luminescent molecule together with the layer of solvent surrounding it, 
this centre having a definite radius for all centres. In the layer, besides the molecules 
of the solvent, there can also be found molecules of unexcited luminescence centres. 
Depending on the number of unexcited luminescence molecules populating the 
layer (active sphere’), all centres are divided into groups. Thus, the centres belong to 
group k = 1 if, in addition to excited molecules, there are no luminescent molecules 
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in the "active shpere'; the centres belong to group Ё = 2 if there is one luminescent 
molecule, etc. The centres belonging to different: groups have different probabilities 
of migration: of energy and decay. The total intensity of light emitted by a definite 
‘group of centres, for example, the k-th, is described by I,(t) from formula (6) or (9). 
If the light emitted by this group is polarized, then the intensity components of the 
emitted light are described (as has been found to be the case) by the functions а, е % 
and bye” *#', or, if we take into account the previous disscussion concerning the value 
of the probability of migration of energy, by Г! ,(t) and 7,,(4) as given by eqs. (7) and (8). 
Thus, every pair of associated exponential functions (i.e. having equal decay constants) 
of the parallel and perpendicular components of the phosphorescent light, represented 
in Figs. 2 and 3 by straight lines, describes the components of the light emitted by 
a definite group of centres. The entire process of light emission by the luminophor 
is treated as a process consisting of the emission of light by the centres belonging to 
different groups. It may be noted that such an interpretation of the phenomenon of 
-phosphorescence is not in accordance with Bredel's interpretation (1955). The pre- 
sent experiments disclosed only those groups whose centres have a sufficiently large 
probability of occurring in solution at a given concentration of luminescence molecules. 
Later on we shall indicate some characteristic properties which allow us to assume that 
the above-mentioned probabilities are described according to Jablonski’s theoretical 
predictions (Jablonski 1954, 1955) (see formula 9). 

Fig. 6 represents the experimental curves I(t) for acridine yellow at a concentra- 
tion of C, = 67.8 g/l and acridine orange at a concentration of C, = 12.6 g/l along 
with the theoretical curves described by formula (9), where only the first three terms 

- of the sum are taken into account. This is equivalent to assuming that, in this case, 
only three groups of centres, k = 1, k = 2, k = 3 play an essential role in the phospho- 
rescent process. The exponents are put equal to those found experimentally 
(Table I — С, and Table II — C,), while the constant у was so chosen as to fit the 
theoretical curve as close as possible to the experimental points. 

As may be seen the results obtained are in full agreement with experiment. 
Finding the values by the method indicated above, and knowing the molecular weight т 
of the dye under investigation, and its concentration C, we may calculate the radius r of 
the ’active sphere’. 

For the curves shown in Fig. 6 the following values were obtained: 

Acridine yellow: 


v=8.0 т = 223 тої С=6.15- 10-2 ст? г=222А 
Acridine orange: | 
»—85 ` m= 265 вши —C—126-10-3g/com r=313A 


The molecular weight m was taken from Beilsteins Handbuch der Chemie (1935). 
Figs. 4 and 5 show the limiting degree of polarization p, versus the mean lifetime 
of a given group of centres. It should be stressed that according to Jabtonski’s theory 
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(Jabłoński 1955) the degree of polarization of different groups of centres tends in 
the limit ¢ — со to various limiting values different from zero. The limiting degrees 
of polarization for successive groups of centres are indicated by continuous lines 
parallel to the axis of abcissae, have been calculated from the theoretical formula (10) 
with the assumption that the ground polarization p, = 45%. With the help of these 
curves we can determine approximately the kind of groups for a given phosphor. 
The lines parallel to the axis of ordinates give those mean values of the lifetimes about 
which the experimental points fall. Figs. 4 and 5 as wellas Table I and II indicate the 
existence of a certain correlation between the ground polarization of a given 
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Fig. 4. Acridine yellow. Limiting degree of polarization versus mean lifetime. 
===®== (б ol SINE 
mA Ory — 61.951. 
—vy— Сз = 99:4 g/L 
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Fig. 5. Acridine orange. Limiting degree of polarization versus mean lifetime. 
== Сі == 5.0 g/l. 
ue 304 9-172 g/L 
—s¢— Сҙ-- 12.6 g/l. 
—y— Сз = 16.3 g/l. 
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group of centres and Из mean lifetime. This means that there should exist some con- 
nection between the number of luminescent molecules contained in the centre and 
the mean time. The greater the probability of migration of excitation energy in the 
centre belonging to a given group, the longer its mean lifetime. These relations should 
not depend directly on the concentration of the luminescent molecules, but only 
on the type of centres. 

As seen in Figs. 4 and 5 the points are distributed considerably far apart. This 
might be partially due to the fact that besides the migration of energy there exists 
an additional depolarizing factor connected with the changes of the orientation in space 
of the axis of the luminescence molecules. This is probable because the investigations 
were made with gelatin phosphors. 

As already mentioned, an analysis of Figs. 2 and 3 discloses some features characte- 
ristic of the Poisson distribution function. With an increase in the concentration of the 
luminescent molecules the number of groups of centres noticeably increases. If we 
assume that there exists a connection between the number of molecules contained 
in a centre and the mean lifetime of the group to which the centre belongs, then the 
existence of an increase in the probability of finding a centre of greater population with 
the increase in concentration can be confirmed. However, the shifting of the 'centre 
of gravity’ (max. of the Poisson function) in the direction of increasing k with the 
increase of concentration cannot be confirmed. This is probably connected with the 
fact that the difference of concentration used to investigate the phosphors was too 
small. If we assume a Poisson distribution, then, because of the above mentioned corre- 
lation, we face a new problem, namely: in what manner do the luminescent molecules, 
contained in the ’active sphere’ influence, according to their number, the probability 
of a transition S from the level M to the level F. The levels M and F denote, following 
Jablonski’s (1935b) diagram of levels of the phosphorescent molecules, the metastable 
level and the initial level of the transition connected with the emission of fluo- 


rescence and phosphorescence (Fig. 7). We can assume, following Jablonski (1935a) 
and Randall (1945), that the probability S is ee 


S=S,e” 7 (11) 
where S is a constant, E the energy difference between the M and F levels, Е Boltz- 
mann’s constant, T the temperature in °K. 

Тһе luminescent molecules contained in the 'active sphere’ may influence either 
So or E, or both values simultaneously. The experiments of Levshin (1927) show that 
in general the concentration of luminescent molecules has very little influence on 
the position of the fluorescence band. It therefore seems that it is S, that undergoes 
change rather than the position of the levels. Further investigations will be able to 
provide a definite answer. 

Figs. 2 and 3 show the total degree of polarization as a function of time. As may 
be seen, the changes of P(t) are conditioned by the emission of light from successive 
groups of centres. Changes in the degree of polarization of light emitted by the indivi- 
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Fig. 6. Total intensity of phosphorescent light versus time. 


Acridine yellow ® experimental curve 


С, = 67.8 g/l. —x— theoretical curve 
Acridine orange 0 experimental curve 
Сз = 12.6 g/l. —х— theoretical curve 


dual groups of centres were not perceived; only the limiting degrees of polarization 
were observed. From the conditions of our experiment we may conclude that the process 
of depolarization terminates in less than 10° sec. 


Summary 


: 

By means of a phosphoroscope of special construction the decay curves of the 
parallel and perpendicular components of polarized phosphorescent light were 
investigated. 2 ў 

The investigations were made for gelatin phosphors at room temperature for two 
dyes: acridine yellow and acridine orange in two concentrations. We were able to 
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Fig. 7. a. Fluorescence, b. phosphorescence c. persistent fluorescence. 
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represent the decay curves as a superposition of exponential functions in accordance 
with Jablonski’s theory. The entire process of the emission light was intepreted as 
a process consisting of the emission of light by the centre belonging то a difinite 
group. The results obtained indicate that there exists a certain correlation between 
the limiting degree of polarization for a given group and its mean lifetime. It was 
shown that the variation of the degree of polarization of the total light of phos- 
phorescence is conditioned by the emission of light from the successive groups of 
centres. 


I express my thanks to Professor Jablonski for the theme of this work and for 
his valuable advice. 


KPATKOE СОДЕРЖАНИЕ 
Фронцковяк М. Кривые затухания фосфоресценций растворов. 


Исследовано отдельно затухание компонент фосфоресценции параллельной 
и перпендикулярной элекртическому вектору линейно споляризованного возбуж- 
дающего света. Употреблено желатиновые фосфоры окрашенные желтю и акре- 
диновым оранжем. Кривые затухания, согласно теории Яблонского, являются 
суперпозициями экспонентных функций с разными постоянными затухания. 
Группам люминесцирующих центров, с разными предельными степенями поля- 
ризации, соответствуют разные средние продолжительности жизни. Изменения 
степени поляризации полной фосфоресценции во времени, обусловлены эмиссией 


очередных групп центров, принимающих существенноое участие в фосфоресцен- 
ции. 
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THEORETICAL CONSIDERATIONS ON THE POSSIBILITY OF 
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In this paper a new way for separating mixtures and isotopes is proposed. This 
method is based оп the superposition and mutual intensification of the thermal diffusion 
(the boundary conditions being different from the usual ones) and distillation processes. 
The phenomenological theory of proposed method is set up. A numerical example based 
on the theory presented shows that there exists a possibility of increasing by several 
times (mean — 15 times) the degree of separation in an isotope mixture at equilibrium 
as compared with the separation obtained under similar conditions when using thermo- 
gravitational or distillation methods alone. 


1. Apparatus and Principle of Operation 


The apparatus proposed in this paper appears to be suitable as the basis of a new 
method for separating mixtures and isotopes. This apparatus will be composed of two 
concentric vertical tubes. А reservoir is provided on top to collect the gas mixture. 
Another reservoir at the bottom enables liquid mixture to be boiled. One tube will 
be cooled to the temperature at which the vapour of the liquid mixture condenses, 
the secorid will be heated to a higher temperature. 

Аз a result of the temperature gradient a thermodiffusion mass flux in the gas 
phase will be present. This flux and the convection current in this phase (caused by 
the dependance of the density on temperature) will produce changes of concentration 
along the tube. Simultaneously, the presence of a thin liquid film on the “cool“ wall 
and movement of the entire of the gas phase relative to liquid film will produce a distil- 
lation process. Thése two separation processes, thermogravitational and distillation, 
if taking place in the same direction, ‘will be superimposed and will TE су one another 
on every elementary area of the tube. | 

It can also easily be seen that even if in the proposed apparatus no distillation 
were to take place (relative volatility — 1) thermal diffusion alone would proceed at 
different, significantly more favourable conditions than usually encountered thermo- 
gravitational method. Indeed, owing to presence of the liquid, the film on the *cool* 
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wall, the enriched constituent can be transported inside the gas phase by the thermo- 
diffusion mass flux in distinction to usual boundary conditions in thermal diffusion, 
where there is no flux through the walls. 

In order to investigate and predict theoretically the possibilities of the proposed 
method, the remainder of this paper will be devoted to setting up a phenomenological 
theory expressing the dependence of the degree of separation onallthe variable para- 


meters. 


2. Basic Physical Assumptions 


(a) We assume that thermal diffusion and distillation take place ın the same di- 
rection. Thus if a component, because of thermal diffusion, flows in the gas phase 
to the warmer wall, we assume that the same constituent flows to the gas phase in the 
distillation process. This assumption is usually confirmed for isotop% mixtures, where, 
in the main, the lighter isotope gathers in the gas phase at the warmer wall and flows 
in the distillation process to the gas phase. There also exists a rather broad class of 
nonisotopic mixtures for which this assumption is true. 

(b) We assume that on the phase boundary there exists an equilibrium between 
vapour and liquid. This assumption seems reasonable owing to the fact that at atmosphe- 
ric pressure and above the rate of interphase mass transfer is considerably greater 
than rate of the diffusion processes transporting the mass inside the gas phase (Rose 
1951, Westhaver 1942). Thus for example Westhaver's theory of the distillation column 
based on this assumption is well confirmed by experiment (Rose 1936). In our case 
this assumption would not be limited — at least not above atmospheric pressure — by the 
greater flow of heat through the phase interface than is usually encountered in rectifying 
columns. Indeed, investigations on the rotary thermal rectifying column where there 
exists an additional heat flow through phase boundary, show that at atmospheric pressure 
and above these columns perform like conventional contactors (Byron 1951). 

(c) We assume that in the thin liquid film flowing along the cooler wall there exists 
a uniform concentration in the direction perpendicular to liquid's surface. The theory 
of the rectifying column in the form of an empty tube based on this assumption has 
been well confirmed by experiment (Rose 1936). Also, the direct calculations of 
W. Kuhn (1954) show that in columns without packing one can neglect the diffusion 
“resistance“ of the liquid as compared with that of the gas phase. | 

(d) We neglect the vertical diffusion in the liquid phase as compared with that 
in the gas phase. The reason for this is quite evident when we take into account the 
order of the diffusional coefficients in both phases. 

(e) We assume that the mean velocity of liquid phase is small as compared with 
that of the gas phase. This assumption seems quite reasonable because of the difference 
of densities in both phases. 

(f) The geometrical construction of proposed apparatus will be such that the 
height between the two reservoirs will be large as compared with the distance between 
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the two plates. The mass of a given constituent contained in the working space inside 
two plates (or concentric tubes) will be small as compared with the mass of the same 
constituent contained in the reservoirs. Owing to this construction and to the fact 
that in the gas phase the horizontal concentration gradient is caused by the thermodif- 
fusion mass flux, (horizontal changes of concentration will be small) we can use two 
assumptions made in Furry, Jones and Onsager's (1939, 1946), well-known theory 
of the thermogravitational method. We assume, namely, that the horizontal concen- 
tration gradient is independent of the height and that in the working space a quasi- 
stationary process takes place. 

(в) We assume further that (а) we have the same thickness of liquid film at every 
point on the “cool“ wall (В) in the gas phase a laminar movement occurs (y) we are 
dealing with a mixture of two constituents (д) and the process takes place periodically, 
so that there is no total flow through the horizontal cross-section of the apparatus. 

(h) We assume finally that the temperature gradient employed and the difference 
in concentration are fairly small, so that changes in the diffusion coefficient, thermal 
diffusion constant, viscosity and density are insignificant. We can easily free ourselves 
from this assumption in this theory, but this would only lead to some complications 
in the final equations. However it appears (Furry 1946) that in the theory of the thermo- 
gravitational method in gases this assumption is usually a good approximation, and 
use of the mean values of these quantities gives good results even when the temperature 
difference is not small. 


3. The Basic Equations and Boundary Conditions 


We take our system in Cartesian coordinates and direct the x-axis perpendicularly 
to the parallel plates, the z-axis vertically upwards and place the origin on the warmer 
plate. We denote the breadth of the plates by b and by ¢ the distance from the surface 


.of liquid to the warmer plate. For the “lighter“ constituent (i.e. flowing to the warmer 


wall) the x-component of the mass flux as caused by thermal diffusion will be (Furry 
1939, 1946) given by 


дс 1 OT 
5,0 + Dac (1 ©) т эх ® (1) 


Jr =—D 
and its z-component by 


VE ا‎ EFT Do 3P vco, (2) 


where c denotes the concentration, D the diffusion coefficient, œ the thermal diffusion 


constant, T the temperature, ọ the density, v the convection velocity. 
In the gas phase there are no sources for the components, and the concentration 


will satisfy the continuity equation 


9 (oc) ‚зу 3 
— + div J 0. ( А 
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As the boundary conditions for (3) we have 


lim J, = 0, (4) 
x0 
and owing to assumptions (b), (c) 
ыны сен ы даа d (5) 


ERE 
xt === = с 


where c; denotes the concentration іп the liquid phase and where the relative volatility € 
is regarded similarly to assumption (h), as constant. Since motion takes place only 
in the direction of the z-axis, the equations of the viscous, stationary and incompressi- 
ble! motion of the gas phase have, in our case the following form 


en ep. | 
d? v, 9 
n 0E (7) 


where p denotes the pressure, n the viscosity, v, the z-component of velocity. 
The boundary conditions for equation (7) are (for simplicity we denote v, by v) 


(0), о = 0; (8) 
by assumption (e) we have 

lim v = 0 (9) 

х= 


Since there is по total transport of the mixture through the cross-section of the 
apparatus, the flow through the gas phase is 


b 


PRES == [ой zb, 
0 


00 
from which we obtain the third boundary condition for (7): 


è 
j vods = ô. (10) 
0 
The convection velocity in the gas phase is caused by the change in temperature 
along the x-axis. Therefore, at this point we must take into account the dependence 
of the density on temperature. For this purpose we differentiate equation (7): 


аз - dT 
Кожа gp di (11) 


1 Vertical pressure changes in the column, of order of several mm Hg per 1 т, (Rose 1951), 
are insignificant. 
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where 
до 
8— 5% 
since 
dT 
since only the heat conductance is taken into account in the gas phase, we have 
AE AE 
c Sta EY: 
finally, we obtain 
ае | — ВАТ 
dx? nd з 


In other parts of the theory o will act as a factor “calibrating“ the concentration 
in appropriate units, hence as a result of assumptions (h) we may treat о as a constant, 
mean quantity. 


4. Introduction of а С Function 


From (3), owing to our assumption that the process is quasi-stationary we have 


divj — o ке 
Or 
9 Jx 9?c дс 
SER et = а 14 
0x De 222 ence dz Y = 


д°с А 
Let us neglect, for the moment, the vertical diffusion part Dos nen we obtain 


from (14) 


Qui ads oos gu 
QR ELA (19) 
We hypothesize that there exists a function G(x, z) such that 
9c №. n (16) 
s G(x, 2) = 7, 
From equations (15), (16) it follows, since 9%c/dxdz = 0, that 
= ER (17) 


therefore 


gue t ү ovdx + fiz), | ш 
0 
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where f(z) is any function of the z-coordinate. 
From boundary condition (4) and equation (16) it follows that 


G(0, =) = f(z) = 0 (19) 


Hence G depends on the x-coordinate only. | 
From equations (12), (17), (19) it follows that the С function satisfies the equation 


аза ағас 
ee — 20 
dx? a dx ) ш. (20) 


As the boundary conditions for equation (20) we have from (4), (8), (9), (10) 


(G),-0 = 9, (21) 
OT 0; (22) 
ас 
ES — 0, (23) 
dG 
(E) «a 2 


Integrating equation (20) four times and taking into account equations (21) to (24), 
we obtain 


G(x) = — E: ох? (x — C)? — 5 x? (36 — 2x) (25) 


where, in accordance with assumption (h) and the last sentence of the section 3, we 
regard о as a constant, mean value. 


5. Final Form of the Equations for the Degree of Separation 


The mass of the constituent with a concentration c transferred in a unit of time 
through the horizontal cross-section of the apparatus is given by the equation 


a 


Ё а t 
m = b | J;dx + b Јах = — Db 292 dx — 0 MU dx + 
| Oz Oz 
0 t 0 


с 
4 a 
+6 f ocvdx + b | отсдх, (26) 
0 
where the index / denotes the liquid phase. Taking into account assumption (d), we 
have 
i t t a 
д 
т = — Db у. б dx 4- b | ocvdx + b | оплсах (27) 
0 б | 


с 
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Let us investigate more closely the integral 
с с 
ore 
"ne Ах = — [| 
fi остах | 2x cdx 
б б 
Integrating by parts, we obtain, using boundary conditions (21), (22) 


4 
P= (6+ с2а 
0 


Let us calculate дс/дх from equations (1), (16) 


9c 9c 1 1 ӘТ 
3,0 = — Des, + Dac(1—9 т, 0, 
therefore 
Ic 1 дс 1 oT 
айна no Gann hen Оро, 
so that 


с 
I WM 1 [G? дс 
0 0 


Using assumptions (В) and (c), we obtain 
t 


8c f “тх бе TEL 
т = —Dbo 52 ах + cb | omıdx — b(Gc),-; + «b Ge (1 — o) т ae dx — 
9 с 6 
t 
ЖЕ » 
ary; G 22 4 (28) 
5 


As there is no total flow of the mixture through the horizontal cross-section of 


the apparatus, we have | 
; : ы 
jl ovdx = — f op,dx = д у; (29) 
0 t 


We note that the concentration in the gas phase changes slowly with the x-coordi- 
nate. By putting c outside the integral and writing сС(2) instead of (сб), We obtain 
an equation expressing the dependence of the mean concentration, relative to the 
x-coordinate, on the height. 


[4 
de. «AT PANDI dc 2dx 30 
= —Dbot =e + cbó — cjbó — ЕТ bc(1 — E Gd. Do di. G | (30) 


where Т' denotes the mean temperature. 
4 : 
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Boundary condition (5) may be written in a more convenient form. It can be shown 
(Westhaver 1942) that for 10 < e < 1.1 (which holds for all isotopic mixtures and 
a broad class of other mixtures) instead of condition (5), one can write with an error 
smaller than 5% 


(с=с). = dne (cl ©). (31) 


Using (31), we obtain from (30) 


t 
n=— [ee + р, | |a Шеш (32) 
0 


where, since С < 0, we use |G|. 


Let us introduce the notation 


Ку = Пьос, (33) 
(5 
ке icis | (34 
оза Do. X, ) 
Н, = böln e, (35) 
(Б 
«ATb 
Hr E (36) 


0 
Calculation of the respective integrals gives 


a(AT)*Bgobl? , «ATdb 


ze Tn IT SD, 
_ 6°ВҰАТ)2оь7 — gBATbO¢* 135662 (38 
: 9! Dr? 6! 7D 35рр ` 
For the gas one can put 
——— Cem 
B 25% T , 
we then obtain, finally 
є a(AT)%go «АТ | 
2552 goAT 
Же ( (АТ) во 7 80 4 Й 
Dbio © OT DT 9! Б?Т? Suc DT ® BED 520; (0) 
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| Denoting by Г, the capacity of the top reservoir and by Г, the volume of the 
liquid in the bottom one, we obtain the following equations for the rate of change of 
concentrations: 


de, m à; 
асаба OTA 
den = 01V; : ao) 
The initial and final conditions for equations (39), (40), are 
hm с = бы 
1-0 
lim т = 0, ш 
>00 


where c, denotes the initial concentration. 
The dependence of the rate of change of concentration in both reservoirs on all 
the variable parameters is given by equations (39) to (41) together with the equation 


d 
m——K^ + Hc — 9); (32^) 


we note here that c(0, t) = cj, c(h, t) = c, where h denotes the height of the column. 

The individual parts of equations (32^), (37°), (38) can be given direct physical 
significance. The term Hc (1—с) is positive; it favours the separation. The first term 
of H gives the separation by distillation; the second by thermal diffusion; the third 
term represents the separation caused by thermal diffusion arising from the new boun- 
dary condition, viz. the existence of a liquid film on the *cool* wall. The negative 
term — Kdc/dz, counteracting the separation, denotes the remixing caused by longi- 
tudinal diffusion, convection currents and motion of the liquid film. 


6. Discussion 


In order to obtain an equilibrium separation we put m = 0 in equation (32^). 
We then have 
de H 
— = (1—0). (42) 
dic KW ( 9 
om reservoirs, respectively, we 


Denoting c,, с, the concentration in top and bott 
(42), a formula for the equi- 


obtain, after carrying out the integration in equation 
librium separation factor q 


DIE: 


сал бул» ACH 
в 1-9 
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If c, < 1 the equilibrium separation factor reduces to the ratio of concentrations 
H 
= (44) 
Cb 
Suppose that in the proposed apparatus we distil in a unit of time smaller and 
smaller amounts of liquid mixture. In the limiting case д — 0. Thus if there is no 
liquid on the “соо]” wall, the apparatus will work as an ordinary thermogravitational 


column. 


Indeed, if 6 — 0 then by (377) and (38’) we obtain 


Hi => 0, 
a (AT? go? 5b 
Haee a 
(ДТ)? g? 03 C" b 
OLDg* T3. 5 
Now we see, with the help of our theory, that in the limiting case д — 0 one obtains 
the fundamental formulae of the wellknown theory of the Furry, Jones and Onsager 
(1939, 1946) thermogravitational method. 
On the other hand, if the difference in temperatures between the walls were 
diminished, in the limiting case AT — 0 the proposed apparatus would work as 


К, > 


a rectifying column. 
Indeed, let us assume that we have an equilibrium and that the difference in 
temperature between the walls is zero: 
m = 0 
АТ-0. 
Thus 
Н —bólns, 
К > P 520% + Обо. 
35 Do 
The height equivalent to the theoretical plate, i. e. the height J of the column, for 
which 4 = &, is given by the formula 
Н 


к [= ше; 
thus 
12 ¿ô? 
35 Do + 220 
[peu mnm 
д > 
but 
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therefore 
13 902 
De 

Thus in the limiting case AT +> О, we obtain for equilibrium without thermal 
diffusion the fundamental formula deduced by Westhaver in his theory of the 
rectifying column (Westhaver 1942, Rose 1951). 

One can deduce from (15), (16) and (10) that the G function defined above 
indicates the existence of a mass flux through the surface x = £. This may be explained 
by assumption (b) which says that the velocity of the establishement of equilibrium 
on the boundary of the phases is much greater than the rate of transport to the inner 
part of the gas phase. The flux J, is able to transfer the component of mixture from 
the boundary to the inner part of gas phase, condition (5) being satisfied quite well. 

In all the formulae for enrichment there appears the distance between the wall 
and the surface of the liquid. Measurements usually give us only the distance between 
the walls. But from the Navier-Stokes equations the thickness of the free liquid film 
on the “cool” wall may be determinated (we assume that this thickness is independent 
of the height) if д is known. Indeed 

18] = fee dx = v, 0, (a — 6). 
[4 
For laminar motion of the liquid film we have (Lewicz 1952) 


l= 


E в(а—)*%о 
4-0 UNE е 


37, 6\” 
-а-( ne)" Я 
801 


In the theory given above we сап free ourselves from assumption (e), but the 
final formulae are then complicated. Suppose, we assume as an exact boundary condition 
that the tangential component of velocity hangs in a continuous manner when we 
pass from one phase to a second. 

We then obtain instead of (24), the following condition: 


AES 
ах ENT = 


where x is the velocity of the surface of the liquid and is equal to v,(£). Hence 


adi 36 (ee ji 
2 01 376 


The G function then takes the following form: 


3 х? 
6) = e аа 01—020 Fa — (3 8 — xe) а. 
4 | 
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To the quantities H, К one must then add the supplementary terms 


‚ — &АТхобь 
Ны а 
‚206%. аВАТохбь ПІ xoc? b 
= 71051) HO 105 D 


It can easily be seen, however, that usually H' << Н and K' << К so that the 
supplementary terms may be neglected. It is also readily seen that when 6 — 0 then 
x — 0 and we obtain, as above, the Furry, Jones and Onsager theory as the limiting 
case of our formulae. In order to investigate the degree of separation that one can 
expect for an isotopic mixture with the proposed method as compared with that 
obtained under similar conditions when the thermogravitational or distillation method 
alone is used, we shall consider a numerical example. We shall employ the following 
values ordinary met in thermogravitational columns 


Тізе 073: Kod ME VOOR аа Т- 5239, б=5х 10-lem. 
озо LOFT poise;! 04 O 108* оош О оО CME gen an EDT. 


As the rate of distillation, we take bô = 3 x 10-3 б зес-1, i.e., а mean of 10 g per hour; 
for an isotopic mixture we take the values of а, e to be: 


@=О IO ves 102 


From equation (34) we then obtain, for thermal diffusion only (д = 0): 


H 
9 5 -2. 
K 2.2 X 1015; 


and for distillation only (AT — 0): 


H 
25 -2. 
K 2.6 x 1072; 


while in the method proposed here: 


= 4.2 х 10°. 

Hence, it may be expected that, by using the method proposed in this paper 
the height equivalent to a thereoretical plate will be half the height used in a thermo- 
gravitational method for an isotopic mixture under similar conditions. 

Let us consider the equilibrium in a 1.5 m. column when c, < 1. In this case 
the numerical example for isotopes which was considered above shows that by using 
the proposed method one can expect to increase the ratio с/с, by about 20 times as 
compared with that obtained under similar conditions by the thermogravitational 
method and by 11 times as compared with distallation method alone. 
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It is quite evident that for non isotopic mixtures the separation would be much 
greater owing to greater values of the quantities œ, e. 


In a future paper we shall consider the possibility of separating azeotropic 
mixtures and of continuous operation by the proposed method. 

The autor wishes to thank Professor К. Gumiriski for his interest in this work, 
Professor B. Sredniawa and Docent J. Janik for their critical remarks and to Mr 
А. Golebiewski for his valuable comments. 


KPATKOE СОДЕРЖАНИЕ 


A. Витковский. Теоретические рассуждения над возможноствю нового метода 
разделения смесей и изотопов. 


Был предложен новый способ разделения изотопов и смесей. Этот метод 
основан на наложении и взаимной интенсификации термодиффузии пробегающей 
при отличных от обычных граничных условиях с процессом перегонки. Выдвинуто 
феноменологическую теорию предлагаемого метода. Численный пример рассчитан 
на основании предлагаемой теории свидетельствует о возможности увеличения 
в несколько раз (в среднем 15) степени разделения изотопной смеси — в сравнении 
с разделением в смежных условиях посредством перегонки или термогравита- 
ционного метода. 
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NEUTRON POLARIZATION AND ANGULAR DISTRIBUTIONS 
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The neutron polarization and angular distribution of the neutrons produced in 
а (р, п) reaction are calculated with the assumption of the direct mechanism of the 
reaction. The cloudy crystal ball potential modified by a spin-orbit coupling is 
assumed as the neutron-nucleus interaction. The Coulomb effect of the incident proton 
wave results in a flattening of the angular distributions and gives smaller absolute values 
for the polarization. The model is applied to the reaction ?9Si (p, n)®®P for E, = 6 and 
6.5 MeV. The maximum negative polarization ~~ — 65% appear for the c. m. scattering 
angle 0 ~100°. Further extensions of the model are discussed. 


$ 1. Introduction 


А) General Situation in the Theory of Two-Nucleon Reactions 
of the Surface Type 


Some of the features of the process of inelastic scattering of nucleons by medium 
and heavy nuclei have been explained on the basis of Bohr's (1936) compound nucleus 
model. In this model it is assumed that when a nucleon of moderate energy strikes 
a target nucleus it may be absorbed to form a compound system of some stability. 
Subsequently this compound system can decompose emitting a nucleon and leaving 
the target nucleus in an excited state. The theory supposes that the formation and 
the decay of the compound nucleus may be treated as independent processes, so that 
the mode of decay of the compound system depends only on its energy, total angular 
momentum and parity, but not specifically on the way in which it. was formed. Hauser 
and Feshbach (1952) have considered this theory in some detail! and the results 
of their calculations account for most of the observed features of inelastic scattering 
for low neutron energies. i 

However, at medium energies (~ 10 MeV) angular distributions for the two- 


nucleon reactions of the types (п, п’), (p, р’), (n, p), (p, п) have been found to be 


strongly peaked. Further, it was stated that the nucleon groups, which loose little 
energy in the reaction, are of much greater intensity than it would correspond to the 


1 Some of the assumptions of this theory have been recently criticized (Eisberg and Hintz (1956). 


; (95) 
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evaporation from the compound nucleus. This led to the statement that in many cases 
the dominating mechanism of the reaction is not the compound nucleus formation, 
but rather one of several types of direct interaction. This is especially marked for 
the cross-sections measured far away from the pronounced resonances of the excitation 
function. Which type of direct mechanism dominates in a given reaction depends 
very sensitively on the particular configuration of the given nucleus. 

For the inelastic scattering of neutrons or protons by some nuclei the most 
important part of the excitation results from the direct nucleon-nucleus interaction, 
both the initial and the final nucleons being identical. This means that for certain 
configurations such as е. g. double magic nuclei the incident nucleon preserves its 
individuality in the process so that the exclusion principle may be neglected in the 
description of the incident and the outgoing nucleon waves. The nucleon-nucleus 
interaction is then represented well by the optical model potential. For some even- 
even nuclei such as 24Mg (where rotational excitation plays an important role) and 
for heavy nuclei the additional direct interaction is the Bohr-Mottelson coupling 
between the incident nucleon and the surface deformations?. This was shown in the 
theory of Hayakawa and Yoshida (1955) and Yoshida (1956). Recently Yoshida (1956) 
was successful in explaining the experimental results of Fisher (1954) and Greenlees 
and Souch (1956) for the ?*Mg (p, p’) ?*Mg reaction. Brink (1955) was successful 
in his theory of the inelastic scattering of neutrons from heavy nuclei. 

For many nuclei at energies — 10 MeV the most important mechanism for the 
reaction seems to be the knocking out of one of the surface nucleons by the incident 
‘nucleon. This is especially important for (p, n) or (n, p) reactions with nuclei whose 
configurations include one or two loosely bound outer nucleons. The idea of the 
direct mechanism of these reactions as suggested by Mc Manus and Sharp (1952) 
formed the basis of the fundamental theory of Butler, Austern and Mc Manus (1953) 
hereafter referred to as B.A.M. This theory developed in impuls approximation 
describes both the incident and the outgoing nucleons in terms of free plane waves, 
the interior of the nucleus being excluded from the integration as in the case of 
Butler's theory of stripping reactions. This means that the nucleus is black for the 
incident nucleon (any deeper penetration of the target at the energies in question 
would result in absorption and compound nucleus formation). The crudeness of this 
approximation becomes apparent especially at higher energies (— 15—30 MeV), 
where the nucleus becomes fairly transparent. Therefore the angular distributions 


of the type Уа. Jall k, Е |Ro) as given by the B.A.M. theory give only a crude 


agreement did experimental data. Neverthless Benveniste et. al. (1956) obtained 
good agreement with the B.A.M. theory for their results on the inelastic scattering 
of 31 MeV protons by ?Be. Satisfactory agreement with the B.A.M. theory was achieved 


2 The same mechanism of reaction was discussed for inelastic scattering of &-particles from 


24 Mg by Gugelot and Rickey (1956) 
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by Schrank et. al. (1954) for the 59Fe(p,p)9Fe reaction with 17 MeV protons and 
not bad agreement was achieved by Cohen (1955) for (р, п) reactions on various 
targets for 23 MeV protons. The results of Eisberg and Igo (1954) indicate that for 
the inelastic scattering of 31 MeV protons nearly the entire process takes place 
by direct interaction and the compound nucleus formation is unimportant. 
Thomson’s results (1956) on the 83Cu(p, n) 6371 reaction at 18 MeV also indicate the 
importance of the direct mechanism of the process. 

In spite of the crudeness of the B.A.M. theory its comparison with experiment 
nearly always allows the derivation of spectroscopic data concerning the orbital angular 
momenta of the final bound nucleons. The parameter L characteristic of the peaks 
of the B.A.M. angular distributions, satisfies the relation I; + 1, > L> Il; — 14. 
There is no interference between different L; thus the position of the first maximum 
determines the lowest vàlue of L, i. e. |; — 1. Hence l; = l; + Lmin Therefore 
information concerning. spins and parities of the final states may be obtained in 
a manner similar to the case of deuteron stripping. 

The Born approximation (with the integration over the entire volume) corresponds 
to the transparent nucleus idea. An example of the Born approximation is the paper 
of Selove on the ?H(n, p) reaction. Demeur (1955) presented a Born type theory of 
(n, p) reactions with Bhatia's approximation. 

Satchler (1955) developed the theory of the (n, p) angular correlations following 
the B.A.M. black nucleus idea. 

The more eorrect theory requires the introduction of the proper boundary 
condition on the surface of the nucleus (instead of the black nucleus.concept) as 
proposed by Tobocman (1954) for (d, p) reactions. This is equivalent to replacing 
the outgoing plane waves by waves distorted by a certain average nucleon-nucleus 
potential. Such a theory is already able to give not only qualitative crude results, but 
good agreement with the data in many cases. Angular distributions for such reactions 
as e. g. 9Li(p, p’)®Li, 12C(p, p’)!2C, ??Ne(p, р)%Ме are reported to have been 
explained with the help of such a type of theory (see e. g. Levinson (1956)). Particu- 
larly good agreement with the data of Schrank et al. (1954) was obtained by Kajikawa 
et al. (1956) for the %Ее (p, p^) Fe reaction at 17 MeV. The distorted waves method 
was the first to give a good explanation of the data of Schrank et al. (1954). 

Very frequently the experimental data were compared with the B.A.M. theory 
and from the comparison conclusions were derived concerning the contribution of 
the direct mechanism in the reaction. Since the distorted waves method makes serious 
modifications of the angular distributions, such a procedure involving the B.A.M. 
theory is not always justified. As it is shown in the present paper the Coulomb effect 
often results in the considerable flattening of the angular distributions. 

Another type of theory for the reactions of this type at higher energies 

| (--20--30 MeV) was presented by Hayakawa et al. (1955). Their semi-classical calcu- 
lation is based on the Fermi gas model. The mechanism of the reaction is the knock 
out of one nucleon of the target by the incident nucleon in direct collision. 
# 
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The same mechanism is also important for some reactions with complex particles. 
Huby and Newns (1951) suggested it for the (d, d") reactions. It is also suggested in 
several papers for (р, «) and (n, œ) reactions. Demeur (1954) gave the theory of the 
(3He, 2p) reactions and of general (x, y) reactions (1955). 


B) Polarization in the Reactions of the Sripping Type 


In the last few years a great deal of experimental and theoretical papers were 
published on polarization of nucleons elastically scattered from nuclei. Polari- 
zation of nucleons from nuclear reactions of the exchange type, on the other 
hand, has been very little investigated thus far. The main reason are very serious 
experimental difficulties in double scattering experiments for these reactions, which 
have rather small yields. 

Theoretical investigation of polarization in nuclear reactions affords various 
advantages: 

1) it provides additional information for nuclear spectroscopy, e.g. in some cases 
the sign of polarization indicates which value of j corresponds to a given / of a state 
involved in a reaction; 

2) it predicts particularly efficient sources of strongly polarized nucleons; 

3) it provides additional information concerning scattering potentials for the 
final state nucleons, and in particular concerning spin-orbit couplings and other 
noncentral couplings. Polarization is sensitive to such couplings; 

4) it may provide additional information concerning the mechanism of a given 
reaction. 

The first theoretical papers were published on stripping. On using various 
assumptions quite contradictory results were obtained. Experiments were necessary 
to choose between them. 

The first of these papers (Newns 1952) is a semi-classical picture of the polarization 
in stripping reactions. Horowitz and Messiah (1953) carried out the distorted waves 
calculation with the assumption of the hard. sphere scattering of final state nucleons. 
Yoccoz (1954) and Grant (1955) performed calculations without the distorted waves 
but taking into account the Coulomb effects. Cheston (1954) used the nucleon-final 
nucleus potential in the form of the optical model potential of Weisskopf et al. (1954) 


Table І 


Neutron Polarization P (in %) in **Si (p, n)*P for E, = 6.5 MeV for: b) no proton scattering; с) ab- 
sorption of proton partial wave: Lp = 0.1 


6 (degrees) 


КЕГЛЕЛЕЛЕЛЕЛЕЛЕЛЕЛЕН 
ШЕТЕНЕНЕЛЕЛЕЛЕПЕЛЕЛЕЛЕДІІ 


с) Р(%) 0 -7.3 | -15.2 | -26.7 | 41.5 | -57.9 | -76.8 | -39.7 | 65.2 


55.6 
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modified by an additional spin-orbit coupling. This type of potential was successful 
in describing elastic scattering of polarized neutrons (see Adair et al (1954), Okazaki 
(1955), Thomas (1955), Culler et al. (1956)). Cheston assumed stripping mechanism 
for the !?C (d, p) 13C reaction at deuteron energies of ~ З MeV. 

Recent experiments of Jentschke (1956) and Hillman (1956) backed up Cheston's 
results?. This was a new confirmation of the optical model potential with spin-orbit 
coupling. Both the sign and the order of magnitude of the measured polarization 
were in agreement with Cheston's predictions and in contradiction to all other papers. 

The only two-nucleon exchange reaction for which the polarization was investi- 
gated was the "Li(p, n)'Be reaction. The polarization was measured and calculated 
(Adair et al. (1954), Adair (1954), Okazaki (1955)) for a rather low energy on and 
near a pronounced resonance. The assumption of the compound nucleus formation 
was used. 

On the other hand it seems interesting to calculate the polarization for the direct 
mechanism of a (p, n) reaction off the pronounced resonances of the excitation function 
in ше medium energy range. 


C) General Remarks on the Example Used in the Calculations 


In the present paper the ?9Si(p, n)?9P reaction will be investigated on the basis 
of the distorted waves theory. The assumed knock out mechanism seems to be particu- 
larly reasonable in this case. The outer odd neutron in the 2951 nucleus is quite 
loosely bound to the double half-magic core 851 (the binding energy of the odd neutron 
is ~ 8 MeV while the average binding energy of a neutron in the #Si core is ~ 16. 
MeV). Thus the contributions from the tightly bound internal nucleons are much 
smaller at energies < 10 MeV. The incident protons “зее” mainly the neutron tail 
of the greatest extent. 

The neutron-final nucleus (distorting) potential is assumed to be in the form of 
the optical model potential with the spin-orbit term. Its imaginary part corresponds 
to a partial reabsorption of the beam of the outgoing nucleons by the nucleus. 

Since the direct mechanism is assumed for the (p,n) reaction the extent of 
disagreement between the predictions and the eventually forthcoming experiment 
might indicate the extent to which the direct mechanism is unsuitable at the proton 
energies considered below, i.e. — 6 Mev. | 

‘Least suspect are the results for small scattering angles, since the direct interaction 
neutrons are strongly forward peaked and the compound nucleus neutrons play 
a minor role at small angles. 

Since the odd proton in the ??P final nucleus (В+, 4, 6 sec) is in a S i5 state, as will 
be shown below, the neutron polarization in this reaction will be entirely due to the 


феде рғ тан Сі A EN a re IT T CUm 
3 however, Cheston's paper contains certain errors (see Sawicki (1957), which might affect 
numerical results considerably. . 
$ 
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spin-orbit coupling in the neutron-final nucleus potential. If the final state were 


1, > 0, polarization different from zero would result even if the LS term were absent 
an polarization is zero for free plane waves only, i.e. no final state interaction of the 
outgoing particle with the residual system provided J, > 0). Therefore the present 
29Si(p, n)?9P reaction seems to be of particular importance, since it might give the 


most direct informations concerning the spin-orbit coupling. 


§ 2. Notation 


J, М, — the total angular momentum of the free particle and its projection 

L, M, — the orbital angular momentum of the free particle and its projection 

Ј Mp Jp и; — the total angular momenta of the initial and final bound states and 
their projections 

Hy lp — the projections of the neutron and the proton spins 

(x), пу (x), h(x) — the usual spherical Bessel functions 

Е, (n, x), Gp (n, х), Н, (п, x) — the respective usual Coulomb wave functions of the 
Gamow parameter n (j,(x) > е1 Е, (n, x)/x etc) 

Е, (x) = М? (ix) — the spherical Hankel function with imaginary argument 

W,g(y)—the Whittaker function 

F 1, д; 2) — the confluent hypergeometric function 


Е, k, — the wave vectors of the neutron and the proton 
V,, — the neutron-proton interaction 


У (1, m) — the single particle orbital 
i 
(2 > Мұ, D m) — the usual Clebsh-Gordan coefficient 


TACAT Ku, (о) — the neutron and the proton spin functions 
Ф, (&) — the core internal wave function 
а, (Г, M,) = Ani Yrm, (О) — the neutron wave expansion coefficient _ 
b* (L,, M,) = Алй» 47 м, (Qz p) — the proton wave expansion coefficient 
Yim, (©) — the E ете as defined e.g. by Blatt апа Weisskopf (1952) 
1 
у.м) = Y (2 5; Мь M, 
Mr, M; 
gular momentun eigen function of the neutron | 
V, (J, Г, М,) = const. Гу, (Е, г) — B (Е, Л АО (Е, 7)] ЕТЕ (0), (LM) 


= const. [jr (kpr) — В (Lp) А (E, r)] Yr, m (25) 
В (Г, J), B (Г) — the о parameters 


5 м, ИИ Хм, (0) — the total an- 


M; M; 

i €i & 

i = ыз = = 1 __ ће wave numbers of the bound states 
M;, М; — the reduced masses of the initial and final fet nucleons 


а 8/-- T respective binding energies 


E 
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ж ж ” 
M,» М, — the reduced masses of the free neutron and the free proton 


* 
25 Ze M; h GGG И) Я 
nee c г — the “internal” Gamow parameter of the bound proton 
Ze M; с „ ^ 
пр = ES — the “external” Gamow parameter of the free proton 


Г(=) — the usual -function 
o,(n,) = arg L (L + 1 + in) = Im In I(L + 1+ inj) 


$ 3. Calculation 


А) Incident Proton Wave without the Coulomb Effect 

The nuclear shell model is assumed for the initial and final bound states in the 
reaction. Following Austern et al. (1953) we assume that the neutron, initially the most 
loosely bound with the spin zero core, in a state of total angular momentum ]; orbital 
angular momentum [; and its projection т; is knocked out from the surface 
of the nucleus by the incoming proton which is captured into a state of total 
angular momentum jp, orbital momentum l; and its projection ту. In the free wave 
(B.A.M.) approximation the matrix element for the (p, n) collision of particles with 
spin in the “outside” region ((г> rj) becomes: 
Ma = f do, da, | Е dr, dr, [ш] e^ д (о) V, Pi [us] едә д, (с) (1) 


(rn >") 


1 ` 
Y, (1 = » ( Б „Mi, Mn 


тр Un 


where: 


% и) У (Ll, т) x,: (с,) Po (8) (2а) 


1 ; 
P; [ug] == Уу (» ә т, Ир Ју» и) V (ly, my) Хи, (0,) Ф, (5) (2Ъ) 
тр ГА 

A simple derivation of Мр for spinless particles is given in Appendix I. The 
generalization for particles with spin is straight-forward. 

As mentioned in § 1 the improved theory consists in replacing the plane waves 
іп Mg (with the “Butler cut-off" at rọ) — by the distorted waves corresponding to 
the proper boundary conditions at rọ appropriate to the scattering problem. The 
derivation of the scattering amplitude corresponding to the distorted waves method 
is similar to that of Tobocman's (1954) for (d, p) reactions*. The distortion parameters В 
related simply to the average reflection factors 2] are determined by the scattering 
potential of the given particle in terms of the logarithmic derivatives of the wave 
functions of the particle at the boundary гу (see Appendix IT). It may appear strange, 
at first sight, to use for the knocked out particle the boundary condition of the scatte- 
ring of free incident nucleon. However, the reciprocity theorem for nuclear reactions | 


justifies this procedure. 


4 The derivation for J — dependent forces will be presented in Appendix IV. 
ПШ " 
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We assume that the outgoing neutron is scattered in a spin-orbit potential. The 


expression en’ y, (o,) occurring in (1) must then be replaced by: 
n 


lul = >, У, 0.M) e ‚Мр, 


Г, Мі, J, My 


5 м) PLM) (9 


To simplify the calculation the incident proton wave is assumed to be distorted 
by a nonspin dependent potential. This simplification seems not to be very significant. 
In general, the influence of the incident proton nuclear distortion on the neutron 
angular distribution and the neutron polarization is shown below to be small (this 
was first stated for the angular distribution only by Austern et al. (1953)). The case 
of an eventual LS coupling for the proton is considered in Appendix Ш. Accordingly 
ep p Xu, (o,) in (1) is replaced by: 

p (шы = 21 b* (L,, M,) Ф (Ly, M,) Xup (с,) (4) 
Lp Мр 

The Coulomb effect of the incident proton wave (4) will be neglected for the 
present. The rather small polarizations and the corrections to the angular distributions 
for (d, n) reactions due to Coulomb interactions as fount by Yoccoz (1954) and Grant 


(1955) ericourage the use of the latter assumption for the proton energies considered 
below. In the next section we shall see how far these assumptions hold. 


If the transition operator Z,, does not contain апу spin dependent terms such 


as 6," 6, we obtain for the partial matrix element?: 


== H (J, L, М; ф (1, my) Kup (6,) | p | Ф Cr М,) Хи, (,) ф (l, т) Xu, (с,) = 3 


1 Я 
== (2 9 5 Мі, и, |J, м ди, < (J, L, Му) v (бту) Ше Ф(1,, Mp) № (l; т) > 
> (5) 


The reaction amplitude corrensponding to the u, component of the neutron spin 
is then: 


M (и; Ир > Ир My) = 22 {J, L, M,, Г» Мы Һа а (1, Mj) b* (Lp M,) х 


il 1 T 
X (1,3 $ Mis m |J, Mi + a) (мень e 


J, M, + 


1 1 Tee 
+) (45 $ Mg — Up Mp |р н) UI En is =) Е МЕ 


= Еа” My) ф (lr, my = Ш; — Ly) Ше Ф (Lp M,) v (1,, тиза се ш) > (6) 


5) This assumption is not essential for /; = [+ = 0. In that case the polarization and the angular 
distribution remain unchanged if V,,, is modified to be Г, (г) + Г, (r) On * бу. Similar result for 
(d,p) reactions and for Born approximation only (P=0) was recently reported by Cheston and 
Werntz (1956). ` | 


uw 
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If the complete spin dependent reaction amplitude is: 


1 
МХ = м( —À Ил = + i qu, (on) + М | m Ln = — the- Y, (os) (7) 


then the polarization vector P is given, as usual, by: 


ic {< My |с, | My >} av 
= {< My | My >} av (8) 


where { } denotes the averaging over the the initial state quantum numbers Ji, 1, 
and summation over jj. The denominator of this expression determines the angular 
distribution. 


Now let us choose a system of reference in which only one component of P, say P, , 
is different from zero. The symmetry properties of the polarization problem indicate 
that the quantization axis z must be then directed along the vector k, X k, i.e. perpendi- 


cular to the reaction plane. On directing the positive x-axis along the vector k, we 
then obtain; 


B * л . л 
5* (Lp Mj) = Anite Үг, м, (2 ; o) : a (L, Mj) = Алі Yim, Е D (9) 


First we shall show that in this system of reference the components P, and P, 
vanish. b*(L,, M,) and a(L, Му) given by (9) have the property that they are zero 
provided that L, + M, # even and L+ М, = even. Let us assume that 


M ( eam + 1) is different from zero. From Eq. (6) її is seen that the selection rules 
1 1 
фот p, = 1 are І+1І,+1, = еуеп, M, = M, + т — ш + 5. Then 


1 К 
both L+ М; апа L, + My + т — т. 2 must be simultaneously even at least 
for one component of the sum occuring in (6). Thus we get : L+ М + lg +l; + 


1 ў 
+ m, — + = = even, i.e, lz т — p; + СЕР even. Correspondingly, for 


the М bs „= 25 amplitude we have (the summed in (6) quantum numbers 
, ) , 

being now denoted by primes): 1/ -- L-+ L, + l= even, M, = M; + ту — i — 

— 5 . Then if M | Un = — 2) is to be different from zero one gets: L' + MAE 


= even and L, -+ M,+ My — и;— чу = even, Le. L' + мМ, + I, + lp + mp— ш 


1 
Ja dana: even. Thus we would obtain J; + J; + m; — uj 5 = even. However, we 


2 
previously had J; + и m; — Ш; — 52 — even — 1 = odd. Hence either L' + М ENS 


£ 
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Р 2 1 1 4 ° 
= odd or L, + М; ty KM o — odd. In any case M mL zero if 


M E + i — ( and vice versa. 


Thus the polarization in the reaction plane is always zero: 


1 
Pip DD 


v2 


RE 


| 2 
bre) 


The degree of polarization is: 


It is easy to show that always P (0 = 0°) = P(0 = 1807) = 0 

Let us now pass on to the integrations in Eq. (6). The single particle orbi- 
tals may be written аз: y(l, т) = const k (t; ^ fre s (P2 and w(l, m) = 
= const f. у (trp) Y, e (0). ). The radial function fy (tr) =W n, 727: 1 (2t, r,)|r, takes 


account of the Coulomb effect of the final state (bound) proton. Лу (terp) was fairly 
well approximated by some equivalent Ak, 7 (t; r,) in the important region of ГА 
А and t, are defined in such а way that Л, (t; rp) and АЁ, y (ty r,) are equal for two values 
ofr, іп the neighbourhood of ro (т, = 70 and тр = 2ге, for example). When /; increases, 
у 09 rapidly to the value t;, the Coulomb Е being less important if there 
is a strong centrifugal potential. The only parameter is t and UA, 27) [1 , (tj ro) = 
= ky, (t; 2ro) Fy (t; ro) may be calculated with the help of the well known series develop- 
ment or integral representation of the Whittaker function (see Whittaker and Watson 
(1940)). t; is always larger than ty. 

| In the energy range considered below it is reasonable to assume И,, in the form 
of the zero range n-p interaction Vip = Год (г, N) similarly as it was done by 
Demeur (1955), Satchler (1955) and Kajikawa et al. (1956). This approximation is 
the better the less the energy. In the last paper it was shown that V,,, assumed in the 
form of a Yukawa well introduces very small changes to the zero range results (the 
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calculations were made for the ®Fe(p, р”)9Ғе reaction at 17 MeV). In this approxi- 
mation the partial matrix element occurring in (6) may be written as: 


< v (J, Г, M, + и, — ha) v (lp, т) |У, np) Ф (Ly, Mp) v (Ll t, — Hn) > = const х 


X û (M, + и, M, + Hn т) (Ге dr ki, (t; г) fig (ter) [Л (5,7) — BL, J) hy? (ka г)] х 


То 


X Uro (hp 7) — В (Ep) Ру» (Е, DIS (12) 


In the next stage we shall confine ourselves to the l; = /, = 0 case. In this case 
the reaction amplitudes may be written in the form: 


M (m >. =2) = у L umm (1.3 E (M "WE 2) 
1, М; 
(9.2 la) eo Mot Qs Mie Dg Mis MEX) 
Ge Meat = UE H (oj ا‎ (13a) 
м(ш->ж--4) = 3 КАЛУ (L, Mj; — 1) (z 5 My, 


1 
+ a (L, МІ) b* (L, Mz) (2 2 ; Mr, 


: оа В |в 13b 
ps. хм, —4] (o oer әз I ( ) 


where: 


Ry -= fr dr ko (fo (tr Li gt) —B (La J) AD (ka Li 06,0 —8 (2 MPG]. (14) 


To 
чы xac Pi mw Жул шш = 


Lh ET RER TE m c MM ee en 
* In taking into account the effect of the nucleus recoil one should replace both Ё, and kp by к= 


A—1 owing to the Coulomb 


in the zero range approximation. However, 


2 7. / . li it 
effect and the scattering terms ~ В one should have kp, k, rather than kp, ky For the sake of simplicity 


ESS i ў 
modification in all expressions. 


A 
we neglect the 
£ 
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In our system of reference (Oz || k, x k) both P, and P, vanish. The polari- 


zation 13: 
ive IE Түр yu 
пм 
(15 

where: 


L,Mr, 
І — Мі 8 
16 
N (16a) 
1 ж л L+ 1— ML + 
L,Mr, 
L+ ML x 
+ (GEM) re} "ic Маш 
where Куу m = Кг. fbe 15 coupling equals zero, then ЁЎ = В; = В, and 


Р = 0. 
Only with L< kr, В gives significant contribution to the М” s. Thus only the 


first several B's must be retained. Therefore it is worth while representing М in the 
form: 


Moo (и) = const | d? ko (1,7) fo (tyr) [e + fit (ke + f* (Е, 0D 


(r>ro) 


where 
Fi (быт) = 
2 

У: а (Г, Mr) (z = Мі, Ал | J: МІ, + 2 (-В (Г, J) hy? (б, r)) Үімі ( 2’) (18a) 
L,ML;Jj 
and г 

f* (hyp t) = У) 6" (1, Mj) — 8 (2 AP Yam (0) (18b) | 

L,Mr, 


The distortion coefficients В (Г, J) will be calculated for the optical model potential : 
in the form of a complex square well: 


=, РОА) 18 -г<в 0: (9 
=) 1 D о 
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The "Butler radius" го was arbitrarily chosen to be the nuclear radius 
Ко = 1.45 413 10-13 ст. В (Г, J) are given by гу as В (Г, J) = = ‚ where 
2 


nış =l i. e В (L, Ј) =0 if no scattering takes place. We may write: 


ги, n-—V04 i9 — 5 Ls (20) 
l. е. | 
1 A = per Pe ДА 
(=1) =i = — VU tiL (21а) 
1 u = A [pu 
“(ь2-:-4)-ш = ҮП) +5 E+ (21b) 


Тһе logarithmic derivatives of the wave function of the scattered neutron inside the 
nucleus is 


fir = Еи; (игу 0), 8, = 14 XGLNAXGIM DX ЛІ 62) 


where 
2M; {ў 
Х* (Г, J) = (kin Ro)? + Xû + it) + 73 | al (23) 
а х2 LM, Y gà Then: 
an 0 — Fa 0: еп: 
m -i UR AOD IERE 
gie d) Ke 


where 6, = arc tg (— Jr (Е, Ro) | т. (kn Ко), Ap +35, 1+ А, МУ (E, Ro) | М? 
(k, Ry; M(L,J)= S, — т Лиз, N (b J) = — A, + Re fi. 


Weisskopf et al. (1954) presented recurrence formulas for the calculation of 7] L in the 
case of no LS coupling, To use these formulas for the calculation in the LS case, 
“опе has to go back to the calculation of fy [X(L, J)] for every new L, J. However, because 
of the complexity of the formulas for higher L, it is worth while computing fry [X(L, J)] 
in this way, at least for the first few L. 


To take into account the incident proton scattering, the first few L, waves are | 
| У i 
assumed to be absorbed by the nucleus | В (L) = EU ) all other partial 


waves being uneffected (similarly as was done by Austern et al. (1953). This should 
overestimate the effects of the proton nuclear scattering. 
The integrals А; у are performed by a numerical method with the use of the Tables 
of Spherical Bessel Functions (1947). 
£ 
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B) Coulomb Effect of the Incident Proton Wave 


To take into account the Coulomb effect of the incident proton wave one has 
to replace the incident proton function occurring in A) by their Coulomb analogues i.e. : 
e»? by F (ky, 1) = e- 2? T (1 + in) e” F, (— inp 1, i (kp r Ё, 1) 
ЛС») by етт») F L(%p> kpr) [kyr 
iz (y) — ВЧ) AY (y) һу ец") [F; (np, ky) — 8° (Ly) Н, (nj, Ат) ky 
Because of the complexity of the numerical computations we set 8° (L,) = 0 (no proton 
nuclear scattering). The computations were carried out with the help of the Tables 


of Coulomb Wave Functions (1948) and Tables of Ln T(z) (1953). The expressions 
analogous to Eq. (17) — (18) were used and the integra: occuring in them, namely: 


Ie = | ат) foltyr) e Е (hy. 9) (25) 
(r2r9) 
was calculated by the method of interpolation polynomials given by Yoccoz (1954). 


If one uses for F (k,, r) the well known integral representation опе obtaines: 


Г sem rn (1 + їл) (2л)! ф (—u) ="?! (1 — uy"ef(u,0)du (26) 
where 


Stu, 0) = | dr keg (tr) fo (tj) ебе n gr (27) 


("> то) 


Ҳи, 0) is an analytical function of u, without singularities. The integration path in 
Eq. (26) us a closed loop starting from и = 1 and encircling и = 0 in the positive 
sense. The path can be chosen as close as we wish to the segment (0,1) on the real 


axis. Therefore we approximate f (и, 0) by a polynomial P,(u, 0) = У) 8,(0)u such 
5-0 
that P,(u;, 0) = f(u;, 0) for a given set of real и; between О and 1. 
If we define the quantities 4,; in such a manner that В,(0) = 3. A,;f (uj, 0) 


then we get: 


ED ( As K) f (u; 0) (28) 


i=0 s=0 
wW hes 


= (= Г + "Гб + 1) Га Ей in, — s] * = (—* іп, ux —]).. 
.. (ing — 5 + Ms). 


For the purposes of computation we use u; — i/n with n — 4. The method of Yoccoz 
was found to be fairly good. 
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84. Numerical Results and. their Discussion 


The ?9Si(p,n)??P reaction (ground states, l; = lp = 0, О = — 5.503 MeV) 
was employed; the incident proton energies were 6 and 6.5 MeV. The binding energy 
of the neutron in 2951 is e; = 8.194 MeV (see e. р. Holt and Marsham (1953)). The 


final state neutron scattering potential was taken to be: 
V. (т MeV) = —40 (12-1005) —2Z-S ~r <R, (29) 
= 0 —т> Ro 


and Ry = 1.45 4"3 10-13 cm. The value of V = 19 MeV (see Adair et.al. (1954) 
and Cheston (1954) is less realistic. 


А) Incident Proton Wave without Coulomb Effect 


Three cases are discussed: а) В(Г, J) = 0, B(L) = 0 (no neutron and no proton 
scattering i.e. the B.A.M. case), b) В(Г, J) corresponding to Eq. (29), В(Г) = (no 
proton scattering), c) B(L, J) corresponding to Eq. (29), В(0) = 8(1) = 12 all other 
B(L) = 0. 


А 3 
For E, = 6 MeV weget B |1, | = В (1 —)=0.0245 — i 0.0097 and В (1, A 


= В(1 +) = 0.0874 — i 0.0124. АП В(Г, J) for L > 1 can be safely neglected. In 
` this case the polarization P is: 

sind (16 E A 

` 000) — cos 0 (167 Re (A*J,(28 (1 +) + 8 (1 —)1 + 


— sin 20 К Лт іш +) (5% ао сеа zl 
аа PEREAT ИЩИ 


| (30 
+ 00820 (609 L7 Pe |8 @ +) (ran +a | 


where 


oo (0) = 224* + 2 (62)* | Л 


1 | 2 : 2 2 1 2 
1ва++@в@—/ АО 


ji 
A = A (0) = I, (0) — 4n (в (o а Io + f (0) Io (i) — В (б. а В (0) ы) - 
— 12 л В (l)Iı cos б, (32) 
where J, = ВО & and 1, (0) = f ko (j) f е" 


("2то) 
со 


I, = f 7 dr ko (tj) fo (y7) М? (hn) jz (Fp 7 


То 
# 
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со 


I, (i) = f rh dr ko (t7) fo (ут) М ЭЛ. (kun): 


E, = f ® dr ko (фт) fo (tr) МЭ (т) AD (kr) 
The angular distributions for E, — 6 MeV are presented in Fig. 1. The only important 
effects of the scattering terms are: the zeros of the angular distributions disappear, 
the first minimum for b) and c) is displaced with respect to a) towards smaller angles, 
the first minimum for c) is very slightly shifted with respect to b) towards smaller 
angles. The curves for E, — 6,5 MeV have been found very similar in shape but steeper. 


d6 


( Relative) 


ОРЕК Nom GOS 2, 
Fig. 1: Angular distribution of neutrons from 29$} (p, n)??P for E, = 6 Mev a) no proton and no neutron 


scattering, b) no proton scattering, c) scattering of neutrons and absorption of proton partial waves 
| Lp = 01 


The neutron polarizations for E, — 6 MeV are shown in Fig 2. The main feature 
of the curves are: 1) P is negative up to 0 — 90? 2) the first minimum of P is 
rather large (— 7 5%) and occurs near 0 — 75°, 3) for 0 > 80? P decreases very 
rapidly and changes sign, then increases in the positive direction until the second 
maximum of |P| is reached. The polarizations for E, = 6.5 МеУ are shown in Table I. 
All the main features of P for E, — 6 MeV remain preserved save that the rather 
high positive values of P appear at large angles. This is mainly due to the fact that 
the differential cross-sections decrease much more at large angles than at small angles | 
with an increase of E. 
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Cheston’s (1954) Р for the 12С(4, p)18C reaction at E, = 3.29 MeV is also negative 
for small angles (contrary to the P obtained by Horowitz and Messach (1953) and 


Newns (1953). Тыз feature of P seems to be characteristic of the LS coupling used. 


р 
05 


03 


o 
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Fig. 2: Polarization P of neutrons from 2951 (р, n)??P for E, = 6 MeV; b) no proton scattering c) ab- 
sorption of proton dartial waves L, == 0.1 


В) Coulomb Effects of the Incident Proton Wave 


The Coulomb barrier of the 2951 nucleus is about 3.3 MeV. For E, = 6 MeV 


Ze?M, 
the Gamow parameter п, = 2 
Ё 


_ 0.904. The Coulomb effect may be important. 


On setting f° (L,) == 0 we get the case analogous to the case b) of A). The respective 
angular distributions are compared in Fig. 3. The Coulomb effect results in a consider- 
able flattening of the curve and shifts the deep minimum towards the larger angles. 
The respective polarizations are compared in Fig. 4. The Coulomb effect results in 
smaller absolute values of P and a shift of the maximum of |P| towards larger angles. 
This shift is rather connected with the corresponding shift of the minimum of the 
angular distribution. The general character of both curves remains the same. 


85. Conclusions 


As regards the angular distributions it is to be noted that it seems rather un- 
reasonable to compare the B.A.M. curves for E, < 10 Me V (case a) of 4A)) with 
experiment and conclude from that about the contribution of the direct mechanism 
in a reaction. The scattering (distortion) effects and, above all the Coulomb effect 
can considerably deform the curves: e.g. the Coulomb effect flattens them. 
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Fig. 3: Angular distribution of neutrons from Si (p, n)??P for E, — 6 MeV, and no proton nuclear 
scattering: b) no Coulomb effect, b,) Coulomb effect included 
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Fig. 4: Polarization P of neutrons from Si (р, n)??P for Ер = 6 MeV and no proton nuclear scattering: 
b) no Coulomb effect. b,) Coulomb effect included. 


The LS coupling in the neutron-final nucleus potential results in polarization 
sensitive to the scattering angle 9 and which is negative in the range of 0:0 < 6 <90— 
—110°. The sign of Р determined by the LS term is not changed by the Coulomb 
effect and the general character of both the curves remains preserved. The maximum 
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absolute values of P are high and the reaction might serve аз a source of strongly 
polarized neutrons. 

Culler et al. (1956) calculated the polarization of 14 MeV neutrons elastically 
scattered from different. nuclei (the potential was used of the type given by Eq. (19). 
V'-the strength of the LS coupling was varied (0.5—2 MeV). Both the angular distri- 
butions and the polarizations were found not to be very sensitive to this square well 
parameter V’. Recently Thomas (1955) calculated the polarization produced in elastic 
scattering for a complex square well and a spin-orbit interaction which is effective 
only at the nuclear surface where it is of the form of a 6-function. The results for this 
surface spin-orbit interaction are qualitatively the same as those obtained for the 


uniform LS term. 


No other types of radial dependence of the LS coefficient have been investigated 
thus far for low and medium energies, even for elastic scattering. In the high energy 
region, on the other hand, the polarization of nucleons elastically scattered from nuclei 
was found by many authors (see e.g. Fernbach et al. (1955), Guliev (1956)) ratcher 


sensitive to the functional shape of the LS term. For a potential corresponding 
to a smooth density distribution Gulev (1956) found that sharp minimum resulting 
for the square well potential for 290 MeV neutrons elastically scattered from 12С 
disappeared in agreement with experiment. It is generally well known that the square 
well potential of the optical model gives too sharp angular distribution curves. This 
is due to the sharp edge reflection effects. 

It would therefore be desirable to perform calculations with a more realistic 
tapered-edge LS term as e.g. the LS coupling as introduced by Blin-Stoyle (1955) 
which is effective near the nuclear surface’. Such calculations, however, could be 
carried out rather only approximately (e.g. by the J. W.K.B. method) and would be 
very cumbersome. 

Extension of the present calculations to the additional coupling with surface 
deformations of the Bohr & Mottelson type seems also to be desirable. 

It seems that the experiment with the ?9Si(p, п) Р reaction would be the most 
desirable. Above all it might provide valuable information concerning the type of 


the LS coupling. 


Some of the results of the present paper were published previously (Sawicki 
(1955), (1956)). 


6 Bjorklund et al. (1956) obtained a'very good fit to the 14 Mev neutron scattering data by making 
use of an absorptie surfacve. 


7 Also different realistic types of tapered edge complex wells ГА (1 + i 6) i.e. [27 (г), С = C(r) 
r—R 


should be investigated with respect to P. The best tailed V = Pg e 4-6 а =) Қа) 193 10:22 
A'^ cm, а = 0.5 x 1071? cm was obtained by Saxon et. al. 111956) 


£ 
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Appendix I 


Simple Derivation of the В. A. M. Reaction Amplitude 
The Hamiltonian of the system is: 


HIT ITE p ар е” (1.1) 


where Т represents kinetic energy, V, and V, are the neutron and proton interactions, 
respectively, with a fixed centre of force. Now let us define yg and y4 by: 


(T, + T, + V, — E) рал —0 (1.2) 
(T, + T, + V, — E wy = 0 (1.3) 
We may now write for W — the solution of the equation (H — Е) ¥ = 0. 
| V = yo — GV, НУ, V (1.4) 
where С is the Green function 
(ТИ — E) 6 = = 8 (5 —7) 9 6, — 10) (1.5) 
i TR - ee! (Tp, A) Ф%% A) (1.6) 


Thus Eq. (1.4) may be rewritten: 
V — ра — G(V, — Ил — G(, + Vag)? (1.7) 


The last two terms give the reaction amplitude 4(6) = Ау(0) --.4,(0). The first one is: 
A) = fe» yh (ry, A) (V, — Vp) е” v, (Fas ғ) dr, dr, (1.8) 


and may readily be shown to be zero (similarly аз was done by Gerjuoy (1953) for 
(d,p) reactions). On replacing У іп A,(0) by y according to the Born approximation 
one obtains: 


Ago (0) = | dr, dr, ат еу (ry, 2) Vy + Veg) er? qi (rss г) (L9) 


The B.A.M. amplitude results with the additional assumptions: 1) V, is neglected, 


2) the nucleus is black i.e. V”? is neglected inside the nucleus: V, 


пр == О for r, < ro: 


————— 
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Appendix II 
Distortion Coefficients В (Г) 


If the wave function of a neutron with angular momentum L appropriate to the 
scattering problem is y, (г) then и, (г) = ғр, (r) may be written in the form: 


i, @ = au) — n, 490) | шл) 

where a = # +1 V2L +1 Т апд | 
49 (уйе = MP (kr) =}, (kr) + in, (®) (1.2) 
ЩЕ) = MP (kr) = jr (kr) — in, (kr) (11.3) 


Then the logarithmic derivative 4. Sle (йї@є (ОБ becomes 


Z MP (b) + m. 2 MP (hr) 

fp =1 +R — gg (11.4) 

TTT PE | 
д | —1 
| New, since ms YL == Me get: 

E hy (kr) — h т In yr. | 

doon, ——|.. (11.5) 

n M dms i „In yr|r | 


Tobocman’ s (954) distortion Socii nts BL) are: 
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Appendix Ш 


LS Nuclear Distortion of the Incident Proton Wave 
On assuming а LS term in the incident proton-initial nucleus distorting potential the 
incident proton wave (Eq (4)) should be replaced by: 


V^, [up] = 3 У b* (Lp, Mp) (>з 95 + Mp» Me | Jo» м) Ф (Jp, Lp» My) (HI.1) 


Jp Му, Lp» Mp 


As in the case of the nonspin dependent V, potential, it can be shown that for 
the 1; = 1, = 0 case in the same system of reference (Oz || k, xk, ,) the reaction am- 


Оден reduce to: 


Moy tn u = У Y, abe MV (L, Mi) (Lg 5 Meran Jo Mi + 2 
L,My J» Jp 
(2 E Мі, Up | Jp; Mr. + 2i R (L, J, Jp) (111.2) 
where: 
R(LIJ)= ТЕ dr ko (t; r) fo (tj 7) Uz (6,7) — 8 (L, J) AL? (km 7)] Ur (kp 7) — 
ў — В, (L, Jp) М? (kp 7)] (111,3) 


The polarization becomes: 
1 \ |2 
У (з) ме) 
Ёр 
ТШ? 
У | Moo (+ 3) aF | Moo (»- 5 


Hp 


| 
| 


On introducing the notation АЁ = К (z Jeep } az =L + z) one obtains: 


P= (Ш.4) 


1 ji | 
Moo Эр » моно ню RI” + 


L,ML 


LETE МЕИЕИ. m E D Mc 
ME [MER pec (o Ar ا‎ RL | (III.5a) 


TEC 1+ 1 
Moo (+ y EE == Т а (Г, МІ) b* (L, М) TE 
242 2 


L+1—Mz\ p++ М КЕСЕ МО * 
Gres | ч! ns A ОЛОТ eet 


«(жем (кезем) ac + (ЫМ) (Me) ас) anm 


2-1 2-1 2-1 2L--1 
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Moo ( m. Ц- ра a (E, мәз» мо (СЕМ) 
DA 


LP Y, 1, 
о = Me) оем. 


L MaN Ete Me My voor MEN DM ine 
== 2 mm | L AE ecu | (Ш.ос) 


po Е І--1- м.\? 
My == - У a (L, Mj) b* (L, М1) ке f INS 
„ML 


т Lo ый мне TEST 
+ ( ЕТЕП еш з иң 


If the f distortion coefficients are small enough so that the products A(L, J )8,(1, Jp) 
may be neglected then Муо(и,, /4,) may be approximately split into the form: 


Moll Hn) = M(u,) + A(uj), (111.6) 
where М is independent of u, and A of u,, respectively. Then P becomes: 


edente npe 0) 
ИЕА 
ele) een 


(ы) PR С 


+ 
In one numerically oversimplified case indication was obtained with the help of 


2 
| R; (115d) 


2 2 


JA Ta) 


+ 


Eq. (III.7) that the introduction of the LS term in И „ results in rather small deviations 
of P from the V, — V, (r) case in smaller absolute values. 


КРАТКОЕ СОДЕРЖАНИЕ 


Савицкий Г. Поляризация и угловые распределения нейтронов с реакции (р, п) 
а оптическая модель ядра. 


Поляризация и угловые распределения нейтронов в реакции (рп) вычислены 
на основании механизма „воздействия прямо” в реакции. Для описания взаимо- 
действия нейтрон-ядро основан потенциал оптической модели, видоизменённый 
путём сопряжения спин-орбита. Куломбовский эффект падающей волны протона 
даёт сплющение угловых распределений и уменьшение абсолютных величин по- 
ляризации. Модель применена к реакции ?9?Si(p,n)*P для Ер =6 и 6,5 MeV. Макси- 
мальные отрицательные поляризации со — 65% выступают для угла диссипации 
Өс100°. Продискутировано дальнейшее решение предъявленной теории. 

£ 
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THE TRANSITION CURVE OF THE ELECTRON-PHOTON COM- 
PONENT OF EXTENSIVE AIR SHOWERS IN LEAD ABSORBER 
OF THICKNESSES BETWEEN 0 AND 25 CM. 
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Institute of Nuclear Research, Department of Cosmic Rays, Polish Academy of Sciences, Cracow 
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Measurements of the absorption curves of extensive shower particles have been 
performed with an apparatus containing a normal extensive shower detector composed 
of three counter trays and a telescope consisting of three large area counter trays. There 
have been registered simultaneously 8 kinds of coincidences with good statistics for two 
positions of the lead absorber relative to the trays of the telescope. Тһе thickness of the 
absorber was changed from 0 to 25 cm. Analysis of the results has made possible the 
evaluation of the ratio p/e, i.e. the ration of the number of photons to the number of 
electrons in extensive air showers. The value found for this ratio is less than 
l and is equal to 0.70 + 0.15. This value has been determined with greater accuracy 
than in measurements of other authors. All results of measurements obtained by different 
authors in connection with the determination of the ratio p/e are analysed. It is shown 
that large numbers of penetrating low energy photons are amongst the particles leaving 
the lead absorber several centimeters thick, and that they may be registered by single 
trays of the telescope, but not by a telescope composed of two or three trays. The regis- 
tration takes place according to the mechanism suggested by Greisen. - 


Introduction 


+ 


The problem of the transition curve in lead of the extensive air shower particles 
of cosmic radiation has been dealt with by many authors. These earlier investigations 
were mainly concerned with the shape of the curve for large thicknesses of absorber, 
namely its slope beyond 10 cm. of lead. This slope has been explained independently 
by Zatsepin (1948) and by Greisen (1949) by taking into consideration the influence 
of the energy dependence of the absorption coefficient of gamma radiation on 
the development of electron-photon cascades in materials of high atomic numbers 
(e.g. for lead this dependence has a distinct. minimum in the energy interval 
2—7 MeV). $ 

, (119) 
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Lately, some authors have investigated the beginning of this curve for lead absor- 
bers, in connexion with the problem of the ratio of the number of photons to the num- 
ber of electrons in extensive air showers. Some of these authors, (Millar 1951, Bassi 
et al. 1952) obtained for that ratio a rather low value р/е < 0.30. Because of the difficulty 
in explaining the low value of this ratio, various authors have carried out further mea- 
surements concerning this problem (Milone 1952, 1954, Bruin 19521, Massalski 1954). 
They have shown that the value of this ratio is higher than that-given by Millar and 
by Bassi et al. This problem has been thoroughly discussed by Milone (1954). 

Some of the methods for measuring p/e entail the covering of one or more counter 
trays with the absorber. In these cases the quantity: 


R(t) = P, (t) + ple: P, (t), (1) 


may be evaluated from experimental data. P,(t) and P,(t) in formula (1) are the pro- 
babilities that an electron or photon réspectively falling on the absorber of a given 
thickness ż will produce beneath it at least one secondary electron which operates the 
counter tray covered with that absorber. The functions P,(t) and P, (t) have been 
evaluated theoretically by Arley (1943), and they are conventionally adopted by most 
of the authors (except Millar) in their work on the determining the ratio p/e despite 
the objections against Arley's theory (the assumed energy spectra of electrons and 
photons are probably not correct for electrons and photons of extensive air showers). 
Knowing the values of P,(t), P,(t) and R(t) опе can evaluate the ratio p/e from for- 
mula (1). 

In the anticoincidence method we make use of the counter telescope, in which 
an absorber is placed between its two trays, and we determine the rate of events when 
the lower tray of the telescope is operated while the upper one is not. The rate of these 
anticoincidences depends only on the expression {p/e * P,}, but does not depend on P,, 
and this decides the superiority of the anticoincidence method. 

The aim of this paper is to present the results of measurements of the transition 
curves in lead of extensive air shower particles. These measurements were made by 
the authors for absorber thicknesses between 0 and 25 cm. We have made some modi- 
fication of the anticoincidence method. Because of high rates of coincidences due 
to large areas of the telescope trays used in this method, the transition effect was 
measured with low statistical errors. Owing to large dimensions of the telescope trays, 
the side effects are negligible, and we could measure in this way the ratio p[e for 
a rather large range of absorber thickness: 0 to 2.5 cm of lead. — 

Since we have established the presence of the above-mentioned penetrating low- 
-energy photons beneath a lead absorber of such a large thickness, we investigated 
these photons for the entire range of thickness over which measurements of the 
absorption of the extensive air shower particles may be influenced. 


1 We know this paper only from references to it by Milone. 
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Apparatus 


We used in our experiment an apparatus which was a modification of arrangements 
used up to now. The detector of extensive showers was composed of three counter 
trays A, B, C, each of an area 5 = 0.468 m?, placed at the vertices of an equilateral 
triangle with sides of 3.4 m. In the centre of the triangle was placed a telescope composed 
of three counter trays D, E, F, each of an area s = г. S = 0.290 m?. Tray E was placed 


GZ ЕСЕСРЕРЗ ; 


(а) 


Fig. 1. Position of lead absorbers in two series of measurements 


right on tray F. A free space was left between trays D and E to allow room for an absor- 
ber 15 ст. thick (Fig. 1). The telescope was enclosed by lead walls 10 cm. thick and 
a lead base 5 cm. deep. All counters were made of brass tubes with 1 mm. walls. The 
measurements were made in Cracow at an altitude of 229 m. 


Principles of Measurement 


Two series of measurements were made. In the first series lead absorbers of thick- 
nesses 0 to 25 cm. were placed over the telescope DEF (Fig. 1а). In the second series 
lead absorbers of thicknesses 0 to 15 cm. were placed between the trays D апа Е 
(Fig. 1b). In each of this series there were registered simultaneously: threefold coinci- 
dences (4BC) — T, fourfold coincidences (TD), (TE), (TF), fivefold coincidences 
(TDE), (TDF), (TEF) and sixfold coincidences (TDEF). The block diagram of the 
apparatus is shown in Fig. 2. In the anticoincidence method we have used an arrange- 
ment in which one tray (D) has been placed over the absorber and a second (E) beneath 
it. The third tray (F) of the telescope was added to eliminate the influence of the low- 
energy penetrating photons with'an energy of a few MeV, which are present in relatively 
great numbers beneath thicker lead absorbers. These photons distinctly contribute 
to registered rates of coincidences when only one counter tray is to be operated under 
the absorber, but have no influence when two trays are operated in coincidence. These 
penetrating photons will be discussed separately. 

All the coincidences may be divided into two groups А and B. Group А comprises 
measurements of transition curves in which the absorber is placed over the telescope, 
containing one, two or three trays of counters. The rates of coincidences of group А may 


be expressed by the formula: 
| o6 
А= CREE mar: аы aa с oci ENDE I р E 
0 > : : 
= C: St {Wa + A, + ple P, + rr] (2) 
$ 2. ` 
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Fig. 2. Block-diagram of the apparatus 
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Fig. 3. Dependence of the ratios of the rates of coincidence ————— and ¬ оп the thickness 


of the lead absorber 


where x is the density of particles of extensive air showers falling on the absorber, 
W, = —3 + 3.27 —3*, Ма) = — + 3(1 + 2)” — 3(2 +2)” + 842)”, Кх 0 *? dx 
is the differential spectrum of densities of registered showers and r’ · r + S is the effective 
area of the telescope used in the measurements. From the ratio (TD)/T in measure- 
ments without absorber we evaluated y = 1.41 + 0.02. The immediate results of 
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measurements of different types of coincidences as compared with the rate of threefold 
coincidences (T) are shown in Tables 1 and 2. In Fig. 3 are given, as an example, two 
curves showing the dependence of the ratios of the rates of coincidences TDEFJT 


and TDEFg/T on the thickness of lead absorber. 


Results of Measurements of the Ratio р/е 


By means of the formula (2) we can find, for all types of coincidences, the values 
of the expression (P, + р/е : P, + k)r'r = R(t) which correspond to different thick- 
nesses of absorber. From these quantities, after subtraction of penetrating component 
k*r'*r and normalization of the obtained curves for the thickness of the absorber 
t = 0 we compute the functions R(t) = P,(t) + ple: P,. The functions obtained 
from different types of coincidences are shown in Fig. 4. 
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Fig. 4. Curves R (t) obtained from different types of coincidences 
] — when the absorber is placed over the telescope and one counter tray is operated beneath the 


absorber 

2 — when the absorber is placed over the telescope and a telescope composed of two trays is operated 
beneath the absorber 

3 — when the absorber is placed over the telescope and a telescope composed of three trays is 
operated beneath the absorber 

1 — when the absorber is placed inside the telescope and one counter tray beneath the absorber 
is operated in coincidence with the upper tray of the telescope 

2 — when the absorber is pla ced inside the telescope and simultaneously two trays of the telescope 
beneath the absorber are operated in coincidence with the upper tray of the telescope 
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To evaluate the ratio p/e from these curves we must carry out corrections for 
the counter wall thicknesses and normalize them anew. [1 this way we obtain from 
group A measurements the quantities R,(t), R(t) and R(t), corresponding to cases 
in which one, two or three trays of counters are operated beneath the absorber. The 
quantity R,(t) has to be corrected for penetrating low-energy photons; we obtain thus 
the quantity R(t). From the values of Rj, R, and В, we calculate, by means of the 
formula (1), the ratio p/e for different thicknesses of the absorber. They are shown 
in Table 3. The most convenient method we consider to be a modified anticoincidence 


Table 3 


Values of the Ratio p/e Evaluated from Measurements 


With modified anticoincidence method 
(from the differences of rates of coinci- 


With absorber placed 
р dences for two positions of the absorber: 


over the telescope ) 
Рр over the telescope and beneath its two 


upper trays) 


(TDEF) 


mean value: mean value: 


0.70 + 0.02 0.69 + 0.04 


method in which we make use of the differences of rates of corresponding coincidences 
registered with an absorber placed over the telescope (group A) and then between 
its two upper trays (group B). These differences (A—B) may be expressed by the 
formula: 


d De Be KR a e е-5*}з И. efle: Py r'r$2) . е-"5х . dx (3) 
0 


The rate of coincidences of group B may be expressed independently by formula (4) 
similar to that which gives the rate of coincidences of the group A: 


B= Като MUT е-5*)з XR e ы ЛАГЫ @х. (4) 
0 


The first term P, in the exponent in formula (4) determines the main contribution 
to coincidences arising from electrons which fall on the absorber and produce beneath 
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it the secondary electrons, which, in turn, operate the counter trays. The second 
term 1 p'le: P,} in the exponent causes only small contribution to registered rate of 
coincidences. It gives an account of the events when high energy photons fall on the 
telescope simultaneously with a low energy electrons; the latter operate the upper 
tray of the telescope and are stopped in the absorber, while photons actuate, through 
their secondary electrons, the counter trays beneath the absorber (obviously p' — p). 
ex a R(t), defined previously, and R(t) = P,(t) + p'le: Р, are shown in 
ig. 5. 


10 25 50 75 mmPb 


Fig. 5. Curves R (t) obtained from different types of coincidences normalized after extrapolation of the 
thickness of counter ‚walls to zero 
Curve К, (4) is obtained from curve ‘1’ 
Curve R,,; (t) is obtained from curve ‘2’. and 73” 
Curve В, (4) is identical with the curve qm - 


Ситуе R, (4) is identical with the curve 21 


We can immediately calculate (p/e* P,} from formula (3) using experimental 
data for R(t). To compute the ratio р/е we need to know only опе theoretical quantity P,. 
The main advantage of this method, as compared with the ordinary anticoincidence 
method, is that we do not need here any normalization of the experimental data at 
t = 0. The values of the ratio p/e obtained by this anticoincidence method are also 
given in Table 3. It is evident from the table that, within the limits of error, the values 
of pje obtained by both methods аге concordant. 

We cannot conclude from our measurements that there exists any dependence 
of the values p/e on the thickness of the absorber. As we have performed our measure- 


і 
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ments for all thicknesses (0—25 mm) of the lead absorber with the same statistical 
accuracy, we may evaluate, from our results, the mean values for both methods of 


measurements. These are shown in Table 3. 


Results of the Measurements of Penetrating Photons 


Zatsepin (1948), taking into consideration the dependence of the absorption 
coefficient on the photon energy, evaluated the cascade curves for electrons for lead 
absorbers. He found that the range of the electron-photon cascade in heavy elements, 
in particular in lead, is higher than that obtained from the simple cascade theory in 
which this effect is neglected. Zatsepin et al. (1948) have shown that the slope of the 
absorption curve of extensive shower particles at thicknesses greater than 12 cm. of 
lead is caused by radiative damping of the particles of the soft component and not 
by the absorption of penetrating particles. 

Greisen (1949) has also been concerned with explaining the slope of the absorption 
curve for particles of extensive showers in lead for thicknesses greater than 10 cm. 
The mechanism he assumed involves the prolongation of the cascade by photons of 
energies of 1.2—7 MeV for which the coefficient of absorption in lead has minimum 
values. These photons have a mean free path of about 2 cm, i.e. several times longer 
than that of electrons or photons of higher energies. Greisen studied the influence 
of this effect on the composition of the registered radiation beneath thick lead absor- 
bers. According to Greisen beneath absorbers more than 10 cm. thick the main remnants 
of high energy electron-photon cascades are photons of low energy (1.2—7 MeV), 
` the number of the latter notably exceeding the number of electrons. 

In the literature there are very little experimental data concerning the immediate 
observations of these penetrating low energy photons beneath thick lead absorbers. 
The first observation was made by Daudin (1949). Bassi et al. (1952) also observed 
these photons in their measurements with lead absorbers of 6 cm and 20 cm covering 
the counter telescope. Jurkiewicz (1954) carried out systematic measurements соп- 
nected with this problem. He showed that the penetrating low energy photons may 
influence the shape of the transition curve of extensive shower particles already for 
rather thin lead absorbers (for t > 2 cm). 

We have thoroughly studied the influence of these penetrating low energy photons 
on the shape of transition curves of extensive shower particles. These curves were 
obtained in absorption measurements for lead absorber thicknesses of 2 cm to 15 cm 
both with single trays of counters and with counter telescopes. 

Let us compare the curves R,(t) and К, g(t) shown in Fig. 5. The first of these 
two curves represents the function R(t) when the absorber is placed over the telescope 
and when one counter tray is operated beneath the absorber. The second curve shows 
the function R(t) when the absorber is. placed over the telescope and when two or 
three trays are operated in coincidences beneath the absorber. As the functions R,(¢) 
and R,(t) are indentical within the limits of error, we denote them by common symbol 
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К, 3(t). Of course all of these functions are normalized for the thickness of the absor- 
Бег == 0. 

We сап see from Fig. 5 that we have relations R,(t) = R,s(t) for t < 2 cm, 
and R,(t) > В, (1) for t > 2 cm. In our opinion, this phenomenon is caused by 
penetrating low energy photons. Beneath thick absorbers, photons of a few MeV are 
the main component of the cascade (Greisen). The number of these photons is so great 
as compared with the number of electrons which may pass through the counter walls 
that they give some contribution to the rate of registering such cascades in the extensive 
shower arrangement in which absorption is measured with a single tray, despite the 
low efficiency of G—M counters for low energy photon detection (е == 0.02). This 
contribution cannot be observed in practice when the absorption is measured by 
a counter telescope composed of two or three trays in coincidence because of the low 
probability of photons causing coincidences in a telescope. 


We can make a similar explanation for the curves R (t) and R,(t). The curve 


R,(t) corresponds to cases in which the absorber is placed between trays D and E, and 
one counter tray is operated beneath the absorber in coincidence with an upper tray 
of the telescope. This curve is connected with events registered primarily as follows: 
electrons falling on the telescope operate the tray D and then through their secondaries 
are detected in a single tray beneath the absorber. On increasing the thickness of the 
absorber we observe in this case some contribution to the registered rate of coincidences, 
this case some contribution to the registered rate of coincidences, this contribution 
due to penetrating photons. This contribution increases for greater thicknesses of 


lead absorber. The curve R,(t) corresponds to events in which a telescope of two trays 
is operated beneath the absorber by extensive shower particles simultaneously with 
the upper tray D (coincidences TDEF), the absorber being placed between trays 
D and E. As was stated above, low energy penetrating photons cannot make, in this 
case, an appreciable contribution to registered rate of coincidences, and therefore 


the ordinates of the curve Р.) аге lower than those of the curve R(t). 

An additional argument pointing to the presence of a rather great number of 
these photons beneath thick absorbers is obtained by comparing the rates of coinci- 
dences (TDF) and (TDEF). If, for these coincidences, we evaluate the respective values 
of R(t) = (P, Е ple: P, + k)r’r which we denote by G(TDF) and G(TDEF), when 
the absorber is put over the telescope, and by G(TDF) and R(TDEF) when the 
absorber is put inside the telescope — then we can write a relation holding for all 
absorber thicknesses (0—25 cm): 


@(ТРЕЕ) 


ОТОР = 09 


where R(TDEF) and R(TDF) are the unnormalized ordinates of the transition curves 
with the contribution of penetrating component of extensive showers. We may regard 
the above ratio as an ’efficiency’ of our arrangement for registering shower particles 


# 
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with the telescope. On the other hand we can write an inequality : 

&(TDEF) _ Ф(ТРЕР) 

G(IDF) ` G(TDF)' 
which for t > 5 cm. is very strong. For example for 2 = 7.5 cm. Pb the value of the left- 
-hand side of above inequality is equalto —0.74 whereas that one of the right-hand 
side is equal to ~0.92. 

We associate this phenomenon with the penetrating photons in the following way. 

If the apparatus is hit by a shower which operates the upper tray D of the telescope, 
and if it happened that from this shower only low energy penetrating photons remained 
beneath the absorber and triggered the discharge in the lower tray F of the telescope, 
then the probability of simultaneous operation of the middle tray Ё by these photons 
is very small. This event may be missed in measurements of the absorption of extensive 
shower particles with the telescope. It may be however observed only when a very 
large number of such photons are emitted from the bottom of the absorber. The con- 
tribution to coincidences when a single tray beneath the absorder has to be operated is 
sufficiently large to show the presence of the effect beyond the limit of the statistical error. 


From the excess of the ordinates of the curve R,(t) over those of the curve R,(2) 
we can estimate, for different thicknesses of absorber, the number of penetrating low 
energy photons which leave the absorber and whih are not accompanied by the elec- 
trons. This number corresponds to one electron passing through the absorber. It may 
be expressed by the formula: 


, R0 — Ryo) 
e+ R,(t) 
where в == 0.02 is the efficiency of the detection of gamma-ray photons with G—M 
counters. These data are collected in Table 4. The values shown in this table are in good 
agreement with those evaluated by Greisen. They have ben used in our calculation 


of p/e to correct the values of R(t) for which the effect of penetrating photons has to 
be taken into account. 


TABLE 4 


Number of Penetrating Low-energy Photons Corresponding to One Registered Electron Leaving the 
Absorber 
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Discussion of Results 


Аз concerns measurements of р/е, three different arrangements are given in the 
literature (Milone, 1954): a) all counter trays of extensive shower detector are covered 
with an absorber, b) only one counter tray which does not belong to the shower detector 
is covered with an absorber, c) anticoincidence method. The two latter methods have 
a main advantage that their results do not need any barometric or temperature cor- 
rections. But of these two methods, we consider as the most convenient an anticoinci- 
dence method in which we make use of only one theoretical value P,, besides its 
results do not need to be normalized. 

Since the results of measurements of transition curves obtained in arrangements 
with the registering trays of the telescope at different distances from the absorber 
have given concordant values of p/e, and since the latter practically do not depend on 
the thickness of the absorber, we may conclude that our measurements are not influen- 
ced noticeably by the scattering of particles in the absorber. From the fact that the 
results obtained are independent of the thickness of the absorber, we conclude that 
the result is not greatly influenced in our case by the change in the mean density of 
particles of the showers registered by our apparatus, connected with the change of 
thickness of the absorber. In our experiment the mean density of the particles is within 
the limits of 10 and 20 per equare meter. 

Аз has been mentioned above, it'is not quite correct to use the values of P, and Р, 
from the theory of Arley. In our case we are using only one of these functions, namely P,. 
But this function does not seem to depend strongly on the assumptions ‘made in the 
theory. The following consideration can give us a rough picture of this matter. From 
the position of the maximum of the transition curve we can roughly estimate that 
the main contribution to registering of coincidences is made by photons of the energy 
of about 108 eV, in addition to electrons. If on the basis of cascade theory (Jánossy 
and Messel, 1951) we estimate w,/w* and Qyax from the position of the maximum of 
the cascade curve for electrons generated by photons, and if we assume, for simplicity, 
the Poisson distribution for the number of electrons, taking namely Pt ax) == 
1 — e-Omex, then we can show that for t,,,,, in the limits of 1 to 2 cascade units, P, 
has a value between 0.45 for & ах = 1 cascade unit and 0.63 for tmar = 2 cascade units. 
We see that this rather rough evaluation of P, is in agreement with the value accepted 
for tnax from the theory of Arley: P,(tmax) = 0.5. Since the statistical error in our 
measurements is small, we assume that above limits define roughly the limits of error 
in the ratio pje. We then take pje = 0.70 + 0.15. 225 s i 

Let us compare the numerical values of the ratio p/e obtained by different authors. 
The result of the experiment of Bassi et al. (p/e = 0.30 and 0.12) differs greatly from 


the values obtained by other authors. It probably results from the poor statistics con- 
Б ---------------------------------- . 
1 Where w is the energy of primary photons, w* — the thereshold energy of registered electrons 

and Qax — the average number of registered electrons in the maximum of cascade. 
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nected with unsuitable geometry of the telescope in this experiment. Millar, performing 
his measurements by method ’a’, obtained for the ratio p/e the value 0.25. Milone (1954) 
called attention to the improper method used by Millar in his calculations. Emploing 
data from Millar’s experiment and using the curves P, and P, from Arley’s theory, 
Milone obtained the mean value of 0.66 for the ratio p/e. Massalski (1954) gives two 
values of pje : 1.0 + 0.3, obtained by anticoincidence method, and 0.8, obtained by 
modified anticoincidence method. The first result has been obtained from the maxi- 
mum of the transition curve. The correction applied for the beginning of the curve 
caused a big error in the value found for {pje * P,j.! If instead of extrapolation to zero 
absorber we use Massalski’s data for the values of P,, as given in Table 1 of his paper, 
we obtain for the ratio p/e the values given in Table 5. 


TABLE 5 


Thickness of the 
lead absorber (cm) 


We thus find an average value p[e — 0.69 on the basis of Massalski's measure- 
ments. 

In conclusion we may say that the results of Millar, Milone, of the present experi- 
ment and also of the experiment of Massalski, when we consider in the latter the cor- 
rection mentioned above, are in complete agreement. In all of these experiments the 
mean densities of the registered particles do not differ greatly. It seems therefore, 
that the value of p/e which we obtain from them applies to the extensive air showers 
with particle densities from a few to some tens of particles per square meter. 

The value of the ratio p|e obtained from experiment is smaller than that one derived 
from the one-dimensional theory of cascades. Probably, this is connected with different 
mechanism for the scattering of photons and electrons in extensive showers. 

The problem of determining the ratio р|е is connected with the first part of the 
transition curve. For lead absorbers of thicknesses greater than 2 cm the shape of the 
transition curve may be influenced by photons which remain after the degradation 
of the energy of the particles entering into the cascade. These photons having an energy 
of several MeV are in the range of the energy for which gamma-ray photons have 
a minimum absorption coefficient. Penetrating low-energy photons have been observed 
in this experiment in direct measurements with G—M counters. We consider the 
agreement of our results, shown in Table 4, with figures given by Greisen as an argu- 
ment for our explanation of differences of the obtained transition curves. 

The authors appreciate the cooperation of Mr. T. Niewodniczaüski during part 


of measurements. Our appreciation is also expressed to Mr. T. Owsiak for his valuable 
help in the course of the experiment. 


1 This has been pointed out by C. Milone in a private letter which is greatly appreciated by the 
authors. 
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KPATKOE СОДЕРЖАНИЕ 


M. Бабецки, И. M. Массальски, Л. Юркевич, М. Менсович. Кривая 
перехода электронно-фотонной компоненты широких атмосферных ливней в свинце тол- 
щиной между о и 25 см. 


Произведнены измерения кривой абсорбции частиц широких атмосферных 
ливней при помощи аппаратуры, состоящей из обыкновенно трхедискового детек- 
тора и телескопа, составленного из трёх счётчиковых дисков с большой поверх- 
ностью. Зарегестрировано одновременно 8 видов коициденций с очень хорошей 
статистикой при двух положениях свинцового абсорбера в отношении дисков 
телескопа. Толщина абсорбера менялась от 0 до 25 см. Анализ полученных 
результатов дал возможность вычислить отношение р/е количества фотонов: 
и электронов в широких атмосферных ливнях. Найденная величина этого отно- 
шения меньше чем 1,0 и составляет 0,70+0.15. Определена она с большей точно- 
стью, чем в измерениях других авторов. Проанализированы все результаты изме- 
рений, произведенных разными авторами, касающиеся вопроса отношения р/е. 
Обнаружено, что проникающие низкоэнергетические фотоны появляются в боль- 
шом количестве между частицами, выходящими из свинцового абсорбера толщи- 
ной в несколько сантиметров, и что они могут быть регестрированы последством 
единичных дисков телескопа, но не регистрируются телескопом, состоящим из 
двух или трех дисков. Они регистрируются в механизме, предложенном Гуей- 
зеном. 
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MULTI-PLATE IONIZATION CHAMBER FOR THE 
DETECTIONS OF SLOW NEUTRONS 


Ву ANDRZEJ BupzaNowskKi AND KAZIMIERZ GROTOWSKI 
Institute of Nuclear Research, Cracow 
(Received November 19, 1956) 


The authors describe the design of a multi-plate ionization chamber with electron 
pulses designed for the detection of slow neutrons. The electrodes of the ionization chamber, 
made in the shape of a multi-plate condenser, are coated with natural boron 3 mg/cm? 
thick. The chamber is filled with argon at a pressure of 1 atm. The chamber efficiency 
for thermal neutrons is about 2%. 


The problem of the detection of slow neutrons is, at the present time, most fre- 
quently handled by using the reaction °B(n,«)’Li. For thermal neutrons this reaction 
has a large cross section which, in this energy region, is inversely proportional to the 
speed of the neutrons. 

Boron is introduced into the ionization Бет of the proportional counter 
in two forms: either in the form of a gaseous compound BF,, which, in this case, 
is the working gas for the ionization detector, or in a metallic form by precipitating 
it on the high-tension electrode of the ionization detector. 

Neutron detectors with ВЕ» gas have a rather large efficiency — of the order 
of 20% — while their disadvantage is that it is quite difficult to purify ВЕ.. The ВЕ, 
gas has a rather large coefficient of attachment of electrons and a small amoünt of 
impurities worsens the situation even more. Аз a result, the work with this detector 
in the electron pulse region encountered considerable difficulty. This gave much 
poorer detector characteristics and decreased the stability of óperation. 

Neutron detectors using thin layers of solid boron may be filled with argon and 
operate in the electron pulse region, and as a result, afford very stable operation. The 
efficiency of these detectors is, unfortunately, very small — on the order of a pro mille. 
The efficiency of this type of detector may be considerably increased by using several 
layers of boron precipitated on the plates of a multi-plate ionization chamber. This 
method has been long used for the detection of fast neutrons with the use of the 
uranium fission reactor. The voltage pulse obtained from the multi-plate ‘chamber 
is rather low when one considers the relatively large capacity of the chamber. The 
fragments of the uranium fission have sufficiently large energy in order for this not 


(135) 
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to constitute a difficulty. In the literature published thus far, one paper describes 
the operation of such a chamber (Lowde 1950). It became possible to build such 
a chamber only after mastering electron pulse technique and as a result of the use 
of proportional amplifiers with a large signal to noise ratio. 
In the present work, the authors were concerned with the building of a slow neutron 
detector characterized by stable operation with the efficiency as high as possible. 
For this purpose, an ionization chamber was built having the shape of a multi- 
plate condensor with a capacity of 15 pF. The high-tension electrode and the collecting 
: electrode of this chamber consist of three 
f d _ plates (Fig. 1). The plates of the ionization 
chamber and the main plate of its vacu- 


-umseal shield were made ‚from alumi- 
nium, which absorbs a minimum of neu- 
trons. 

In building the chamber, particular at- 
tention was paid to decreasing the micro- 
phonics, precise location of the electrodes, 
which, for a relatively large capacity, can 
have a significant effect. 

The plates of the chamber were coated 
on both sides by a layer of powdered boron 
with a surface density of 3 mg/cm?. The 
boron was precipitated onto the plates from 
a suspension in alcohol to which a drop of 
Canadian balsam was added. The chamber 
was filled with argon with an addition of 
5% CO, to a pressure of about 700 mm. Hg. 
The voltage characteristic of the chamber 
is shown in Fig. 2. The chamber works at 
a voltage of 850 V. 

The chamber is connected to the elec- 
tronic circuit, the schematic diagram of 
which is shown in Fig. 3. The pulses from 
the chamber are amplified by the three 
stage preamplifier and a seven-stage amp- 


Fig. 1. Ionization chamber lifier, Los Alamos model 100 (Elmore and al. 
a — collecting electrode 1949). The doubly screened preamplifier 
b — high-tension electrode is located at the chamber and the signal 


с — through insulator (set on araldite) 


а — chamber shield with alumini Б 
кй. кте ЖИ ЖЫШ um top ed cable. The collecting electrode of the 


UGE for aiulditing chamber is connected to the grid of the first 
g — calite insulating rods valve of the preamplifier, the grid resistance 


is fed to the amplifier over a doubly screen- - 
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of which is 100 МО. The preamplifier is built with 6 AK 5 valves with small interelec- 
trode capacities. А delay line is located at the input to the amplifier which forms pulses 
5 psec. long. The overall amplification of the preamplifier and amplifier is 105. The 
amplifier has a rise time of 0.5 usec. The pulses from the amplifier are transferred 
to a Schmitt discriminator and then to a scaler. Before conducting the measurements, 


d 


2000 


7000 


100 500 V 


Fig. 2. Voltage characteristic of the chamber 


Fig. 3. Block diagram 


I — preamplifier 


II — amplifier 
Ш — discriminator 
IV — scaler 


V- — high-tension supply 
VI — square wave generator 


VII — oscillograph 


the equipment was checked with a square wave generator. The square pulses were 
applied to the high-tension electrode of the chamber. The amplifier and preamplifier 
valve filament currents are supplied by accumulators, the anode voltage is stabilized 
electronically. The high tension for the chamber is filtered by a 10 uF condensor moun- 
ted on the chamber. 

_ The efficiency of the multi-plate ionization chamber for a parallel beam of thermal 
neutrons is about 2%. This chamber was used for making cross-section measurements 
_ for thermal neutrons over a period of more than 2 months, and during this time no 

change of efficiency was observed. The efficiency of this chamber can be considerably 
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increased by coating the plates with boron enriched in isotope 10B. It is also possible 


to increase the number of plates in the chamber. 
We desire to thank Professor Н. Niewodniczanski for his valuable comments 


and suggestions during this work. 
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SURVEY OF THE HYPERFRAGMENT EXPERIMENTAL DATA* 


Ву A. Ғплркомузкі, J. GIERULA, P. ZIELIŃSKI 


Cosmic Ray Department, Institute of Nuclear Research, Warszawa 


(Received December 5, 1956) 


All the available cases (about 120) interpreted as hyperfragment decays underwent 
a selection in effect of which 72 cases were considered to be hyperfragment decays, whereas 
in the remaining ones there is a possibility of the existence of other interpretations. On 
the basis of such a “pure sample‘ of hyperfragments certain statistical conclusions in 
reagard to the production, decay, binding energy, structure, and life-time of light 
hyperfragments was reached. 

Tables facilitating the estimation of the binding energy of the A° hyperon in light 
hyperfragments have been drawn. 


Introduction 


During the last 4 years since the discovery of hyperfragments more than 120 events 
interpreted as possible examples of hyperfragment decay have been already published. 
On the basis of such material one may try to reach some conclusions of statistical 
nature, provided that one is dealing with a sample not seriously biased by an admixture 
of other processes. Phenomenologically the overwhelming majority of events interpre- 
ted as possible examples of hyperfragments corresponds to double stars recorded in 
nuclear emulsion. Other processes which yield a similar record in emulsion, and may 
bias the sample to some extent are: collisions leading to nuclear reactions and all 
kinds of o-stars (capture of negative mesons and hyperons). 

To reduce such a bias all the available experimental material was submitted to 


+ 


an appropriated analysis and selection. 


Selection criteria, 


АП the events in the material discussed here were divided into two groups: the 
group of hyperfragments and the group of possible hyperfragments. 
In order to be classified in the first group an event should satisfy at least one of 
the following conditions: 
В Ве ынты ктоо ОЕ аа 
ж Cosmic Ray Conference, Budapest, September 1956; communicated by J. Gierula. 
f (139) 
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(A) collision and o-star interpretation should be both excluded, 

(B) a well defined decay should support the hyperfragment interpretation. 

All the remaining events were classified into the second group of possible hyper- 
fragments. 

Discussing condition (A) we used data concerning the evaluation of the kinetic 
energy of the secondaries and also data concerning the interconnecting track: its 
scattering, thinning down, 6-rays, direct evaluation of the charge and mass, evaluation 
of the kinetic energy at the center of the secondary star. We eliminated collisions 
when we decided that there was good reason to consider the secondary event as happen- 
ing at rest or nearly at rest, or that the kinetic energy of the interconnecting particle 
was considerably insufficient to produce the secondary star. Similarly we eliminated 
o-stars when the direct evaluation of the charge of the interconnecting particle was 
considerably greater than 1, or when the direct evaluation of its mass was significantly 
greater than 2 m, or else, when additional information excluded its interpretation as 
being z^, K- or Y- particle. 

Discussing condition (B), we used also the information concerning the detailed 
analysis of the secondary star. We considered an event as fulfilling condition (B) 
when either: two or three tracks of the secondary star were colinear or coplanar respec- 
tively and the momentum balance of the charged secondaries consistent with all the 
observed data led to an interpretation of the decay process in good agreement with 
the hyperfragment interpretation of the event, or, when an admission of one neu- 
tral particle led to ап *uniquely defined* hyperfragment decay intepretation. 

Obviously this selection unavoidably was to some extent arbitrary. 


Results of the selection 


According to the selection criteria discussed above, 72 cases were classified in 
the hyperfragment group! including all the mesonic decays: 46 cases, and a part of 
the non-mesonic ones: 26 cases. The group of the remaining cases of possible hyper- 
fragments contained thus only non-mesonic decays?. One may expect that the group 
of possible hyperfragments may contain a large admixture of other processes. This 
supposition is to some extent supported by a comparison of the energy spectra of the 
fastest proton in a non-mesonic decay of hyperfragments, and of'possible hyper- 
fragments. Their spectra (Figs. la and 1b) differ considerably. 

Supposing that the reaction 


A’+N-n+N + 175 MeV 


1 References — see tables. I, II and Ш. 
2 Bacchella 55, Baldo 55, Brisbout 54. Castagnoli 55a, b, Ciok E b, 55 b, Crussard 56a, Eisen- 


berg 55, Fry 55 a, b, c, 56, тш. 56, Gottstein 55, 56, Grilli 54, ‘Hornbostel 56, Ladu io: 55, Lovera 53, 
Noon 54. 
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is responsible for the non-mesonic decay of hyperfragments, one would expect to 
observe an emission of fast protons with an energy spread around the value of some 
80 MeV — in accordance with the histogram for hyperfragments. The corresponding 
histogram for possible hyperfragments contains a large number of cases with low 
proton energy. Some of these events may be caused by non-mesonic decays with the 


emission of two fast neutrons. 


In our further discussion we restrict ourselves to the 72 cases classified as hyper- 
fragments — excluding all the remaining ones, the possible hyperfragments. We 
believe that the sample thus obtained is nearly free from bias introduced by any 


r 
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admixture of other processes. However, the relative numbers of different types of 
hyperfragments and different decay schemes present in the sample are, no doubt, 
strongly biased by selection, scanning, and observation conditions. Any conclusions 
concerning these numbers, if considered at all, must be taken with greatest reserve. 


Production 


There are two processes leading to hyperfragment emission: collisions of high 
energy particles with the nuclei, and capture of negative hyperons and heavy mesons 
by the nuclei. 

In the sample here discussed, 60 hyperfragments are produced in collisions, 
while 11 — in the process of capture. 

Collision. Among the 60 cases of collision, 52 hyperfragments are produced 
by cosmic ray particles, 2 — by 6.1 BeV protons, 4— by 1.5 BeV and 3 BeV л mesons, 
1 — by K- meson and 1 — by K^ meson. It is characteristic that our selection does 
not yield any example of hyperfragment emission produced by 3 BeV protons — 
although more than 20000 stars of this kind have been examined for hyperfragments. 

The 52 cosmic ray stars are characterized by a large evaporation (Fig. 2) and 
a rather small number of minimum tracks (Fig. 3), indicating that the more or less 
central collisions of not too energetic cosmic ray particles with heavy nuclei of the 
emulsion are responsible for the hyperfragment emission. Estimation of the primary 
energy using the formula: E, = 37N, + 4 №, „2 + 1000 п, yields the histogram 
(Fig. 4. The energy of hyperfragment emission? is rather low, as seen from 
the histogram (Fig. 5). 

There is scarce information concerning the angle of hyperfragment emission 
versus the direction of the primary, and their associated emission. Out of the 17 cases 
in which this angle was evaluated, 12 hyperfragments were emitted forward, while 
5 — backward. In 5 cases the associated emission of a heavy meson and а hyperfrag- 
ment has been found. à 

The frequency distribution in Z for hyperfragments (Fig. 6) is probably strongly 
biased Something about this distribution can be said considering the charge versus 
range distribution for hyperfragments (Fig. 7). This distribution indicates that hyper- 
fragments with larger Z have rather shorter ranges. The selection criterion (А) has the 
tendency to eliminate hyperfragments with short ranges, whereas, the selection crite- 
rion (B) — those with large Z decaying mostly non-mesonically (Fig. 6). 

Taking this into account it seems that in fact the ratio of hyperfragments with 
greater Z produced in collisions to those with small Z is greater than it would appear 
from the histogram (Fig. 6). 
qM SON dee MU BANTUR nS. Сомасы с d RE aa 


3 The correction for E. for 2 > 1 was calculated according to Baroni G., Castagnoli C., Cortini A.» 
Franzinetti C., Manfredini А., CERN BS9 and Reynolds Н. L., Zucker А., Phys. Rev. 96, 393 (1954)- 
The time of flight was based on the tables by Fay H., Gottstein K., Hain K., Suppl. Nuovo Cimento 11, 
234. (1954) with suitable corrections for higher Z. | 


Survey of the Hyperfragment Experimental Data 145 


15 75 
Q % 
5 & 
E © 

- 10 SD 
59 

= 5 > 5 

qup M 25-135 0 20 30 


Eneryy of the primary Energy per nucleon 


Fig. 4. Fig. 5. 


MESONIC 


С nor mesonic 


Number of events 


230 4 2 


Fig. 6. Fig. 7. 


Capture. Атоп the 11 cases of capture, 4 hyperfragments are produced by 
negative hyperon capture, while 7 — by negative heavy meson capture. 
A fast reaction as for instance 


SS TOCA O 
K+N>a+ Л°+0 


may be considered here as responsible for hyperfragment production. The capture of 
the A°-hyperon by the residual nucleus of the o-stars may here lead to unidentified 
$ 
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hyperfragment decays. It is of interest to stress that in favourable cases the capture 
process allows a good identification of hyperfragments yielding independent in- 
formation both from the analysis of the o-star and the hyperfragment decay star. 


Decay 


The following processes are considered as responsible for hyperfragment decays 
1) 4° >p + rn +37 MeV 
2) A’—n+n° +40 MeV 
3) A? +N >п- N + 175 MeV 


The material selected here can be wholly interpreted in terms of reaction 1 — the 
mesonic, and 3 the non-mesonic decays, excluding reaction 2. No experimental proof 
exists in favour of reaction 2. It is necessary to stress, however, that this mode of 
decay can hardly be expected to be noticed, as it seems reasonable to suppose that 
in the majority of such cases the decay process would be recorded in the form of 
a large-angle single scattering, thus escaping attention. 

Mesonic decays. Among the 46 mesonic decays, we decided, obviously to 
some extent in an arbitrary way, that 27 are “uniquely defined” (Table I). For the 
remaining 19, alternative possibilities of the decay process interpretation were con- 
sidered (Table II). 


The 27 “uniquely defined" decay schemes аге as follows: 


1. ЗН Не + x- — 4 cases 
3H — 10 cases 2. 3H >d +p +a —3 , 

3. 3H —2p +n +” —3 „ 

1. 4H — Не + 27 — 6 саѕеѕ 
4H — 9 cases 

2. Н —3He-d-n + л7--3 cases 

1. He —?He +p +x — 2 cases 
4He — 3 cases 

2. ÀHe—d +4 2р +x — 1 case 
5Не — 4 cases 1. 5Не — “Не +p + x^ — 4 cases 
711 — 1 case 1. Li ts Be dis» — 1 сазе 


Classifying all the 46 mesonic decays we can distinguish the following types: 


1) hfr > fr + p +m — 21 cases 
2) hfr > fr + a7 dan, 
3) hfr > fr +n + m- Бы 


4) other more complex or less definite — 7 , 
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Table I 


"Uniquely defined" mesonic decays of hyperfragments 


Disintegration 
scheme 


3H -?He+n” 

3H +3Не+ мо 

3H -3He+n” 

ЗН -ЗНе-+{ x^ 

3H -d+p+nr 

ЭН "арат 

SH >d >p zc 
?H^ pp tnic 
ЗН >p+p+n+n 
ЭН pd ping: 
4H -1He-4- n7 

ІН -#Не+т- 

4H -1He-4- x7 

aH -4Hr-4-- 

4H “ШОП 

4H +4He+ x7 

4H -?He4-n4-x- 
Н Hence 
4H —4He+n +r 
«Не-ЗНе--р--л” 
«Не -ЗНе-р--л” 
“Не-р-р--4--л” 
5Не >4He+p+7- 
5Не +*Не+-р-+т- 
5Не -4Не+р+т- 
5Hr -^4Hr4-p4-x- 
Li -’Be+7- 


B © ноя 
Binding E Е Pat Р 5 References 
ЗЕ 882 ЕЕ БЕ 

0.7 + 1:1 39.4 + 1 А, В | Bonetti 54a, b 

0.3 + 4 39.8 + 4 Te А,В | Debenedetti 54 
—02+1 40.3 Bn B Baldo 55 

144445 | 38.8 + 4.5 А, В | Castagnoli 55a, b 

0.6 + 0.6 31.8 FW B Fry 56 
—0.5 + 0.6 23.4 + 0.5 FW B Fry 56 

Т | Bus ов аас 

0.4 + 0.7 27 F,W | В Fry 56 

1.2 + 0.6 24.5 Bk B Schein 56, Haskin 56 
—0.9 + 2 26.8 + 2 Bk B Haskin 56 

2.2 +1 52.1 +1 Вп В Haugerud 55, Oslo & Copenhagen 55 

1.2 + 2.1 53.1 + 1 К Вп В Davies 55, Friedlander 55b, —— 55 
—0.4 + 2.4 54.5 + 2.4 Вп В Crussard 55 56b 

3.3 + 1 51 А, В | Schein 55, 56, Goldsack 56 

1.6 + 3 52.8 + 3 Вп В Hornbostel 56 

0.6 + 0.9 53.6 + 0.9 Bn A, B Crussard 56b 

12-09 | 30541 B Ladu 55, Bacchella 55, Goldsack 56 

1.2 + 0.7 30 Bn,S A,B Rosselet 56 

1.8 29.1 Kt Bn B Friedlander 55a, Goldsack 56 
то 305 В | Sorrels 55 

1.5 + 0.6 28.1 + 0.5 B Gilbert 56 

1.8 + 0.6 Bn А,В | Oslo & Copenhagen55, Fowler 56 

25 1 2540005 А,В | Hill 54a, b, 56 | 

1.9 + 1.0 21.6 Вп В Crussard 55, 56b 

2.9 + 0.5 22.4 + 0.4 B Oslo & Copenhagen55, Fowler 56 

2.0 + 0.6 29.5 FW | В Егу 56 

4.4 + 0.7 37.2 + 0.7 B Gilbert 56 


Bk — Barkas W. H., University of California Radiation Lab. Report UCRL-2579 Rev 
Bn — Baroni С, Castagnoli C, Cortini С, Franzinetti C, Manfredini A, CERN BS9 
F — Fay H, Gottstein K, Hain K, Supplemento al Nuovo Cimento 11, 234 (1954) 
S — Ѕіеівегі Phys. Rev. 83, 474 (1951) 

У — Vigneron L. J., Jour. de Physique et Radium 14, 145 (1953) 

W — Wilkins J. J, G/R664 (unpublished) | 
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_ The energy spectrum of z--mesons from three and more body decays is given 
in the histogram (Fig. 8). 


Non-mesonic decays. Out of the 26 cases of non-mesonic decays of hyper- 
fragments we considered only two cases as “uniquely defined”, i.e. | 


“Не >р-+4-+ п | 
HC — He + Li + p 
the others being less defined (Table Ш). 


All, but the two non-mesonic decays, are not inconsistent with the bound 4°- 
-hyperon interpretation. Аз stated above, the energy spectrum of the fastest proton 
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from these decays (Fig. la) suggests that the reaction 
A? +N >n +N + 175 MeV 


may be responsible for the non-mesonic deccays of 4° — hyperfragments. 

In the two exceptional cases, the observed energy release in the decay process 
considerably exceeds the admissible limit for a A°-hyperfragment. Both these cases 
are not inconsistent with а K-hyperfragment interpretation. 


Binding energy 


The binding energy В, of a A?-particle іп a hyperfragment with a mass number 
А can be defined as the difference between the sum of the masses of a fragment 
with the mass number 4-1 and that of the A°-particle, and the sum of the masses 
of the decay particles and their kinetic energies: 


Ва = (ma! + тл) (Dim; + Q4) 


or: 
Ва = (mg -14 Mg — Ут) — Qkin? 
B, 2-3 Qo a Qrin 


Tables of the О’ quantity for different hyperfragments and different decay schemes, 
(tables IV, V) were calculated using the Segré atomic mass values and Friedlander's 
value of Q, = 36.9 MeV. Using these figures the В, for different events may still be 
systematically biased due to differences in the range-energy relation and corrections 
as applied by different authors. 

Тһе weighted mean valuds of В 4 for different “uniquely defined" hyperfragments 
calculated according to the author's limits of error data are tabulated below: 


Hfr B, in MeV Number of cases 
SH 0.5 + 0.3 10 
4H 1.6 = 0.4 8 
“He 1.6 + 0.4 4 
5He қ 2.5 + 0.3 4 
711 Б 4.4 + 0.7 1 
AG 14.6 + 6 1 


These figures must be taken with due reserve, and considered as а rough evaluation 
only, until a recalculation of all these results using uniformly standardized methods · 
of analysis will be applied. Even before such a standardization the data discussed 
here seems not to be in variance with the assumption of unique B, values for hyper- 
fragments of defined А and Z values (Fig. 9). In accordance with the charge inde- 
pendence hypothesis of A°-nucleon forces В, for 4H and “Не do not differ within the 
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“Uniquely defined" decays 


Seeman 55 


Fry 55c 


114,7 + 2,7 
104 + 6 


4Не +p+d+n 
11C -3Не--714--р 


"Non uniquely defined" decays 


3 | Ww, СН od phe 2 He 106-121 97 +4 Ciok 54, 55b ? 
4 | Оз ТЕ! R 2 | 4 2 ~88 ~44 Solheim 54 4 
5 | Ww—Ph, 2.5 | < 0.5 | В 2 |5He 2 592-71 494-1,7 unpublished 5 
6 О и 113-123 | 10-4 Ceccarelli 55a, b 6 
7 | Loi, БЕТТЕ ЧЕ oos 54 48 + 3 Тїйтап 53 7 
8 | Ww—Ph; 129? 3.4 | «0.5 R 3 145 72 70-- 125 60 -- 10 unpublished 8 
9 |Pd, 639 | 51 ң<051 В |3—4 м 39 ~60 23 Baldo 55 9 
10 10 1 Е а-а ЕЛ 3 Brisbout 54 10 
и | Wwe R 21269236) 0.6 | К 3| 3—4 "Li 8 44 40 Ciok 54a, b, 55b 11 
12 | Wseg К. 11 + On Са Re ЗА 20 9,8 Ету 55с 12 
13 | Коз К. 8 + On 6.5 OAR Bean Belus ~12 6,5 Castagnoli 55a, b 13 
14 | Rog TU 5--9p | 41° 1.3 0.8 | В BED 66 33 (100”) Castagnoli 55a, b 14 
15 | Wsc, К. 17 + 5? 3.6 |<0.5 | R 4 |8Ве 2 ~105 58 — 75n Fry 54Ъ, 55c 15 
16.003 В. 27 + 4p BSR 4 |®Be 2 20 — 30 140 = 170" Gottstein 54, 55, 56 16 
17 Re, TU 18 4- 4p «0.5 5 R 4 Ве 3 ~40 6,6 Crussard 54 17 
18 | Pd, К. | Па Т 07| В (4—5 | %Ве 4 112 E Grilli 54, Baldo 55 18 
19 | Ww, К. оороо ка В 100 82 Danysz 53a, b, Ciok 54b, d, 55b| 19 
20 | Mnn, Ho 250-0105]0530:7 909 90 -0517Һ ив 136 4- 15 125 4- 10 Freier 54 20 
21 | Pd; .В. 12-ЕЭр| 51° 4.3 EOS | В 4—7 |"B 2 ~9 8,2 Baldo 55 21 
22 и; R: 8 + On 05 R клу рыу 753 —90 —50 Barret 54 29 
93 .R. 12 + бр 4,8 | <0.5 В БОВЕ 5 ~110 79 Hirschberg 55 23 
24 | Lu, .В. out = <2 E 9 id oe Waldeskog 54 24 
Anomalous decays 
25 | Вов CR; col | 26 + 9n 27 | 24 R 1B 3 139 +168 9,7 K A A Castagnoli 55 | 25 
26 | Wscy C.R. col | 22 + 9р 8 Da eer | в 4 | 340 300 + 125 | В б 


if non mesonic dec. 


~84 | E =45418 


if mesonic decay 


ie wt = 
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limits of error. The dependence of В, from the mass number А of the hyperfragments 
is shown on Fig. 10. It suggests in the range of masses considered here an increase 
of B, with the mass number corresponding to a slope of about 1 MeV per nucleon. 
For the 16 not uniquely defined cases, where an evaluation of B д With a reasonable 
accuracy was possible, there always exists a possibility to choose an interpretation 
yielding B , values consistent within the limits of error with the result discussed above. 


B, 


10 4 Mev 


f. О d Re eee ae ШИ А 
Fig. 10. 


In all the better defined hyperfragment cases: the 29 “uniquely defined", and the 
16 yielding some information concerning the B, and interpreted in accordance with 
the B, versus mass number relation discussed above (Fig. 10), the structure of the 
hyperfragment can be put down as follows | 


hfr = bound nuclear fragment + A? 


No convincing evidence exists for the production of hyperfragments with a structure 
differing from this scheme. 

There is no convincing evidence concerning the existence of the hyperdeuteron, 
moreover, the decay of such a particle, if it were not rare, should rather not escape 
observation. 

All the identified "He hyperfragments decay according to the scheme 

5Не > tHe + p+a, - 

so that it seems that in this case the presence of the /1?-particle has little influence 
on the structure of the residual fragment. ; 
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The lifetime of hyperfragments 


The time of flight of hyperfragments decaying at rest permits us to estimate 
roughly thé lower limit for the lifetime of hyperfragments: 


for singly charged > 10711 sec, 
for multiply charged >> 1071? sec. 
Six cases of hyperfragment decays in flight yield the following values for the time of 
their flight 
H—3 x 107° sec, (cloud chamber), 
3H — — 10-12 sec, 
4He—8 x 10-1 sec, 
4He —5 x 10-10 sec, (cloud chamber) 
О вес) ! 
К == 2x 102% sec. 
These figures seem to indicate that the lifetime of the A, hyperon is not strongly 
changed when bound in a light hyperfragment. 
We wish to express our deep gratitude to Professor M. Danysz for his stimulating 


discussions, help, and valuable advice. Thanks are also due to all the Authors who 
have kindly sent us the preprints of their papers prior to publication. 


KPATKOE СОДЕРЖАНИЕ 


A. Филипковский, Г. Геруля, П. Зелинский. Обзор зкспериментальных 
данных, касающихся гиперфрагментов . 


Все, доступные авторам случаи (около 120), истолкованные как распады гипер- 
фрагментов, подвергнуты соответствующей селекции, вследствие которой 72 слу- 
чая признаны распадами гиперфрагментов, тогда, как в остальных случаях, KOH- 
статирована возможность существования иных интерпретаций. Основываясь на 


так полученном ‚чистом образе” гиперфрагментов, сделаны некоторые статисти-. 


ческие выводы относительно выработки, распада, энергии связи, структуры 
и продолжительности жизни легких гиперфрагментов. 

Помещены таблицы, облегчающие определение энергии связи гиперона A? 
в легких гиперфрагментах. 
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LABORATORY EQUIPMENT AND TECHNIQUES 


MEASUREMENTS OF THE INTENSITY OF TRANSIENT MAGNETIC 
FIELDS BY THE FARADAY EFFECT 


Ву J. Mareckı, M. SURMA AND J. GIBALEWICZ 
Institute of Physics, Polish Academy of Sciences; Adam Mickiewicz University, Poznan. 
(Received July 23, 1956) 


Transient magnetic fields were produced by discharging a battery of capacitors 
through a coil. The magnetooptic Faraday effect was used for measuring the intensity 
of transient magnetic fields. The change in light intensity produced by a rotation of the 
plane of polarization was recorded using a photo-multiplier and cathode ray oscillograph. 
Magnetic fields of intensities of about 100 kOe were measured; CSa, ССІ, and H,O were ` 
used as optical media. 


Introduction 


Raoult (1949) applied an optical method for measuring the intensity of magnetic 
fields produced by discharging a battery of capacitors through a coil. This method 
is based on the magnetooptic Faraday effect. The rotation of the plane of polarization 
created by the magnetic field in the optical medium was measured in Raoult’s arrange- 
ment by a photocell and galvanometer. Magnetic field intensities of about 50. kOe 
were obtained and measured by Raoult. The value of the intensity of the magnetic 
field H can be determined from the magnetooptic Faraday effect if the rotation of the 
plane of polarization a, as well as the Verdet constant V and the optical length l of 
the medium are known. In the present paper a method of measuring the rotation @ 
caused by transient magnetic fields is given. 


Description of the Method 
The intensity of the magnetic field produced in a coil through which capacitors 
are discharged is given by the equation 
R 
= н с 1 
H= Hy sin exp | 2). (1) 
(151) 


Г 
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where w denotes the frequency of the discharge, t the time, К the resistance and L the 
inductance of the coil. 


Using the Verdet law ж = VIH we have: 
Ў R Р 
% = 40 sin wi exp Zar » 
where 
Xo = VIH . 


If the intensity.of the magnetic field is H = H max the rotation value reaches @ nox 
with 
| (2) 


max* 


The value of Ны, can be computed from (2) if the rotation @,,,, is measured. 
The value of & ах was measured as follows: a monochromatic light beam passes the 
polarizer, the optical medium, the analyzer and enters the photo-multiplier connected 
to the cathode ray oscillograph. If the polarizer and analyzer are parallel the intensity 
of the light entering the photo-multiplier is Jọ. The cell containing the liquid being 
examined is placed inside the coil between the crossed polarizer and analyzer. 


Fig. la ij Fig. 1b 


л 3 
Fig. 1. Oscillogram |a: = < Фах CT Di n= ar 2 


The dotted initial parts of the curves drawn in oscillograms in Fig. 1 and 2 were not recorded on the 
films because the brilliance intensifying device was started with a delay. 


4 


. At the time t = 0 the intensity of the magnetic field inside the coil is Н = 0 
(1) and the intensity of the light entering the photo-multiplier is J = 0. Let us assume 
that the maximum intensity Ha, is such that in a given optical medium the rotation 


л л a 
amounts 10 @ = > + 9 where 0 < 9 < 5. At the time of the maximum magnetic 


field the intensity of the light entering the photo-multiplier is Г = I, sin? & ax 
л | ) | 
= sin* | — 2 g |. The change in the photo-multiplier anode current resulting from 


the change = light intensity Г can be registered by a cathode ray oscillograph. The pic- 
ture on the screen of the cathode ray tube was photographed. Fig. 1 fium the change 


of intensity of the light entering the photo-multiplier when = — rax e ап! 
Лада < а : EE 


N او ی‎ a 
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From the oscillograms as shown in Fig. 1 it is possible to calculate æ nay by measur- 
ing h and Anax Since the photo-multiplier has a linear characteristic, the anode current 
is t = kI where Ё is the multiplier sensitivity. 


If the value of œ increases from 0 to @ nax = 


wla 


+ ф, the change of anode current 
is 4i=i,+1, 

F 
where 


A == у == дй 
2 


тах? 


i, = kl, sin? ф= ph, 

hax being the maximum deflection of the spot on the fluorescent screen corresponding 
л 

to an angle of rotation of 3 h is the deflection corresponding to 9 (maximum field) 


as explained on the diagram; p is the sensitivity of the device. 


Hence 
— arc sin В 
ф = arcs ERA 
ei s 
тах = 2 + arc sin hu E 


When „= п A + ф where n = 1,2... the value of ф can be determined in an 


and 


analogous manner; the number of maxima corresponding to Ai = i, determines the 
number n. From (2) we have 2 


л. 1 [oh 
"uer arem Т 
HE ——— M. (3) 


Vi 


л С . . . Б . 
Нах < — it is also possible to calculate the maximum field intensity. In this 


2 


case the analyzer is rotated by a certain angle Ф< а их in the direction of the rotation 
produced during the first quarter of the field period. 
With this bias, the change of the photo-multiplier anode current is given by: 


Ai = ig + iy, 

where | 
ig = kl, sin? Ф = ph, 

апа | 
iy = БІ sin? y = phy, 


Ф denoting the bias angle and ф = Amas — Ф. 
f 
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sin? Ф 

From this equation we obtain = а 

| hy sin? yp 

V = arc sin £c Фф. 

ho 
Hence 
hy . 
@ + arcsin We sin Ф | 
Ho MNT : (4) 


The bias angle Ф is read from the analyzer scale, and the value of hy and hg are 
taken from the oscillogram. 

hy corresponds to an angle of rotation of y and hg corresponds to Ф. The oscil- 
logram obtained with a bias angle is shown in Fig. 2. 


i л 
Fig. 2, Oscillogram (в... < 3 — with a biasing angle 


Apparatus 


The schematic diagram of the apparatus used by the authors for measurements 
of the intensity of the transient magnetic fields is shown in Fig. 3. 

The light source S is a d.c. point lamp operated at 12 volts and about 50 watts. 
The light beam passes through condenser C, microscope lenses L4, L, and the inter- 
ference filter F transmitting light of wave length A = 572 mu. The monochromatic 
light passes polarizer P, cell T, analyzer A with scale, and is focused by a lens L, 
on the photocathode of a multiplier. Glazebrook prisms were used as the polarizer 
and analyzer. 


Fig. 3. Schematic diagram of apparatus. 


The photo-multiplier VpG-69 is connected through an amplifier to the cathode- 
-тау oscillograph. A trigger sweep generator and a device for intensifying the bright- 
ness of the spot were used. The photo-multiplier output is regulated by a potentiometer. 
А glass tube 0.6 cm in external diameter and 10 cm in length was used as the cell. 
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The tube is closed at both ends by glass windows about 0.01 cm thick. The effective 
length of the optical medium was calculated from the topography of the magnetic 
field by graphical integration. 

The field was mapped by using a search coil connected to a cathode-ray oscillo- 
graph. This method has been described by A. Piekara and J. Malecki (1956). Transient 
magnetic fields of high intensity were produced inside a coil through which was 
discharged a battery of condensers having a capacity of 2uF and charged up to 20 kV. 
Before the measurements were made the photo-multiplier was set at the required 
sensitivity, after which the polarizer and analyzer are crossed, or the analyzer is rotated 
by a known angle Ф. At time t = 0 (beginning of the discharge) the trigger sweep 
generator and the spot intensifier energized. The screen of the cathode ray tube is 
photographed. The values of h and Ray, or hg and h,, are determined from the 
oscillograms using the Abbé comparator; hence, by applying (3) or (4), the value 
of H max is obtained. 


х 


0 
rer ОЗА BO: М2, а 16: 718: 20 2 x4. 2069/8: 710255122915 5116 ee 18 20 
Voltage (kV) Voltage (kV) 


Fig. 4 жем Fig. 5 
Fig. 4. The rotation & лу of the plane of polarization versus condenser voltage U. 
Fig. 5. The maximum field intensity Нуу versus condenser voltage U. 


Results 


The points on the CS,-diagram in Fig. 4 represent the angle of rotation Са 
in CS, placed within а coil consisting of 220 turns in 5 layers and having an internal 
diameter of 0.8 cm. 

Тһе effective length of the medium 1, as obtained by graphic integration based 
on the magnetic field intensity distribution within the coil, was [= 4.7 cm.; the 


^ min 
Verdet constant V = 4.4 HE as 


"E 
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The rotation in CCl, was measured inside a coil of 138 turns and of 0.6 cm. dia- 


meter. The effective length 1 for the CCl, medium was 2.3 cm. and the Verdet constant 


V = 1.64 x 10? . The results of measurements in CCl, are represented on 


i 
Oe cm 
the CC],-curve in Fig. 4. From the curves in Fig. 4 the value of H max can be computed 
using (2). 

The curve in Fig. 5 shows the results of measurements of transient magnetic 
fields intensity obtained with both the optical and search coil methods. The method 
of measuring transient magnetic field intensities by using the search coil is described 
by А. Piekara and J. Małecki .(1956). The coil within which the transient magnetic 
field was produced had 220 turns and an internal diameter of 0.8 cm. As an optical 
medium, distilled water was used. The effective length of the medium was T= 4.1 cm. 

min 
Ое ст ' 

Verdet’s constant for CS,, ССІ, and H,O was taken from Landolt-Bórnstein tables 
(1950). The curve in Fig. 5 shows that the results obtained both by the induction and 
optical methods are in very good agreement, the error being within the limit of 1%. 


and the Verdet constant V = 1.4 x 10-2 
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KPATKOE СОДЕРЖАНИЕ 


E. Малэцкий, M. Сурма, Я. Гибалевич. Измерение напряженности импуль- 
сного магнитного поля зффектом Фарадея. 


Разряжением батареи конденсаторов через катушку получались импульсные 
магнитные поля, напряженность которых измерялась с использованием магнито- 
оптического эффекта Фарадея. 

Изменения интенсивности света вызванные вращением плоскости поляриза- 
ции подвергались записи посредством фотоумножителя и электронного осцилло- 
графа. В качестве оптической среды применялись С$2, CCl, и H20. Измерялись 
магнитные поля до 100 000 Ое. 
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LETTERS TO THE EDITOR 


ON THE 3D ELECTRON BAND IN BODY CENTERED STRUCTURE 


By MACIEJ SuFFCZYNSK1 
Institute of Theoretical Physics, University of Warsaw 


(Received November 19, 1956) 


The tight binding method can serve in energy band calculations as an interpolating 
scheme in the manner indicated by Slater and Koster (1954). In the very detailed 
discussion given by these authors for the 3d electron band in body centered structure 
the nearest neighbour interactions are taken into account. When one assumes the 
two-centre approximation and takes nearest and next nearest neighbour interactions 
into account one has more parameters, while the secular equation can still be solved 
analytically in several directions in the Brillouin zone. These directions of high sym- 
metry will be enumerated below. The notation of Slater and Koster will be followed. 
From Tables II and III of their paper we have the symmetric energy matrix Te 


Ни = (xy|xy) = — 4r cosé cosy cos È + B, cos 22 + В, cos2 т] + B, cos2 6, 
H, = (xy|yz) = Ag біп & cos 7 sin б, 
Нуз = (xy|zx) = А, cos € sin 7 sin б, 
Hn (90757) =9 | 
Hyg = (ху|322 — r?) = 24; sin € sin y соз €, | 
Н» = (yz|yz) = — A, cos & cos 7 cos ЕВ, cos 26 + B, cos 2 n + В, cos 2 C, 
Hy = (yz|zx) = Аз біп E sin 7] cos 0, 
Ни = (yz|x? — у?) = үз Аз cos Ё sin 7 sin Ё, 
Н» = (у2|3=? — r?) = — Аз cos 6 sin 7 sin 0, 
Hyg = (zx|zx) = — A, cos & cos y cos С + В» cos 2 € + B, cos 2 + В; cos 2 $, 
Ни = (zx|x? — у?) = — ҮЗ A, sin 4 cos 7 sin C, қ 
Н» = (2х2? — r?) = — Аз sin $ cos y sin % 

ч (157) 
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Hy, = (x? — y?|x? — y?) = А, cos 4 cos y cos £ + В, cos 2 & + В, cos 2 n + B, cos 2C, 

Ну = (х2 — у? 32 — r?) = В, cos 2 £ — Bg соз 2 y, 

Нь = (322 — r?|3z2 — r?) = A,cos&cosncos& + В, соз 2 £ + B5 cos 2 y + Въсоз26. 
Неге mee) = а (kx, ky, kz) the lattice constant being а. 


The 47 and B's are the following combinations of two-centre integrals 


8 


[3(ddo), + 2(ddn), + 4(ddé)s], 
Ay = Ê [-3(ddo), + (ddr), + 2(ddé),), 
Ay = м (адл), — (446), 


А = = [2(4ал)з + (dd9),], 


B, = 445), В, = 5 [3(ddo), + (445, 
B, = 2айл), В, 22. (440), + 3448, 
В, = 2(ddo),, = у V3 (dd), + (da 


When the nearest and next nearest neighbours are taken into account there is no 
possibility of solving immediately the secular equation over all the ¢ = 0 plane but 
still the secular equation in that plane splits into two factors: a quadratic equation and 
an equation of the third degree: 


(Hog — E) (Hy, — E) — Н, = 0, 
(Ни — E) (Нц-- E) (Ags — E) — His] — НН — Е) = 0, 


where in the matrix elements ¢ = 0. 
Likewise in the plane = 7 the secular equation can be split into 


(Hy, — Hyg — E) (Ни — E) — 2H}, = 0, 


and an equation of the third degree 


Hy — E, Hy, Hs 
2 Hiz Н», + В. == E, 2Н = 0, 
Ho Hen Haske 


where now in the matrix elements £ = 7 is understood. As the plane & = 7 is per- 
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pendicular to the ¢ = 0 plane, the possibility of finding swiftly the roots in such 
a plane facilitates the calculation of the density of states. 

Along the directions of highest symmetry the secular equation splits into 
equations of the second degree at most. The roots are 


for £— n = 0, ЕО, 0 <х<л, that is for the Д line 
E = — A, cosa + Ву cos 2 a + 2B,, 
Еҙҙ = — 4, соза + B, cos 2 æ + B, + B, 
E, = .A,cosa% + В, cos 2 « + 2B, 
E; = Aycos@ + В, cos 2 « + 2B,, 
FOR pes == @® =, 0 <а<-, that is for the X line 


1 
йыз = zu — А) соз + (B, + В;) cos2« + 5 (В, + Bs) F 


al 
= |} [A + 40 costa + 2 0 — B) cos 2a + В, — Bilt + darsi 


Ез = — Ау cos?« + A, sin?a + (В, + В.) cos 2а + By, 
Е, = A, cosa + 2B, cos2a + B,, 


ог = == = ДД] <a < that is for the A line 
E, = [245 — (A, + 24,) cosa] cose + (B, + 2B,) cos 2а, 
1 
Eon => [(A, + A, — Ay) cos x — Ag] cosa + Е (В, + Bj) + В, + 5, | cos 2% F 
+ 2 [((A, + A, — А») cosa + A) cosa + (2 (Bs — Ba) + 
2. % 
+ В, — B4) cos2a]? + 6.Аҙ sina cos?a( › 


for £ = = 2 5 =а, 0<a<5, that is for the D line ` 


Б» = 5 (By + В) cos 2a — (В, + B) F 


2 
+ р [(B, — В.) costa + 2(В; — Bj]? + 443 sea] ; 


Ey, = В, cos 20 — В, — В; + А» cosa, 
E, = B, cos 20 — 2B,, 


160 _ М. Suffczyäski 

for 6 -- 0, А (4-5 (); «ac. that is for the G line | 
1 ; 

Еш- 2 (41 — 4) соѕёх + (В, + B,)cos2a-+ 5 (B, + Bj) F 


F i [(А, + 44) costa + 2 (By — В,) cos 2a + B, — B3]? + 443 sintal d | 


Е, = А, cos?a F А, sin?a + (В, + Во) cos 2% + Bs; 
Е, = — А, соѕ?а + 2B, cos 2a + By, 


и =a, =л—а@а, 0<a< E that is for the F line 
E, = (А, + 243) cos?a — 24;| cosa + (B, + 2B,) cos 2a, 


Iac $4 — A, — А,) cos*a + Az] cos& + Б (B, + B) + B, + в, соз2@ | 
+ |: O + А, — Ay) соға + Ay) соға + (2(By — В+ 


+ В, — B) cos 2a]? + 64} sina esta] a; 


É- 


Here the upper signs uses] to first subscripts. The n are chosen - 
in such a manner that they agree with the numbers of d-wave functions in the I’ point. = 
r3 S Ne d E 


R Slater, 7: С. and Kal Gr Ы ый 94, 1498 (954). 
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ON THE WIDTH OF THE 3D ELECTRON BAND IN IRON 


Ву МАСІЕЈ SUFFCZYNSKI 
Institute of Theoretical Physics, University of Warsaw 
(Received November 19, 1956) 


The simplified tight binding or LCAO method with two-centre approximation 
has been used to discuss the energy band of 3d electrons in body-centered iron. The 
discussion is closely analogous to the one given by Slater and Koster (1954) for body- 
centered structure. Here the nearest and next nearest neighbour interactions are taken 
into account. The values of the energy integrals have been computed (1956) in Fletcher's 


Fig. 1. The Brillouin zone for body-centered cubic structure. 


(1952) two-centre approximation from atomic nonorthogonalized wave functions cal- 
culated for iron by J. H. Wood (1955). The signs of these integrals are therefore the 
same for the nearest and next nearest neighbours. The width of the 3d band turns 
out to be 2,47 eV with nearest neighbour interactions and 2,88 eV with nearest and 
next nearest neighbours. The width here is the difference between top and bottom 
of the energy bands calculated along directions of highest symmetry in the Brillouin 
zone (Fig. 1). | 

(161) 
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The variation of energy along these directions is shown in Figs. 2—7, А and B. 
Figs. А refer to the nearest neighbour interactions only and are therefore essentially 


the same as those given by Slater and Koster. Figs. B show the differences which 


arises by introducing second neighbour interactions. One notices that the sym- 
metry of the energy band about the midpoint of energy is suppressed by introducing 
these interactions. In the figures the energies are in electronvolts, and as usual &, 7, 6 


denote 4 a(k,, Ry, k,). 
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Fig. 2--7. Variation of energy of the 3d band along the directions of highest symmetry in the Brillouin 


zone: 
A with nearest neighbour interactions, 
B with nearest and next nearest neighbour interactions. 
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EINE NEUE GESETZMASSIGKEIT AUF DEM GEBIETE DER PLASTISCHEN 
DEFORMATION 


Ву К. WorRUBA 
Ustav pro technickou fysiku, Praha, CSR 


( Eingegangen am 10. Dezember 1956) 


Beschreibung der Versuche 


Diese wurden an flachen Nickelbändern kommerzieller Reinheit durchgeführt, 
die ein einheitliches Ausmass von 100 x 2 x 0,1 mm hatten. Die Bander wurden 
zuvor in technischem Wasserstoff bei etwa 780 Grad C geglüht und sodann jeweils 
einer definierten plastischen Deformatión im Bereiche von 0 bis etwa 30% unter- 
worfen. Letzterer Wert war die Bruchdehnung. Nach der Dehnung wurde an jedem 
Muster die Anfangspermeabilität in einem Felde von 0,045 Oe gemäss Bild 1 und 
sodann die Koerzitätskraft gemäss Bild 2 gemessen. Leztgenanntes Messgerät stellt 
ein von uns entworfenes neuartiges Magnetometer mit magnetischer Sonde dar. 

Wenn wir den Kehrwert der Anfangspermeabilität in Abhängigkeit von der 
Koerzitätskraft gemäss Bild 3 graphisch auftragen, ersieht man, dass der Zusammenhang 
innerhalb der unvermeidlichen Messfehler ein linearer ist. 

Wir glauben uns berechtigt, die geschilderte Tatsache auf folgende Weise zu 
deuten: 

Wie bekannt, hat die Koerzitätskraft (Kondorskij 1937, 1949, Kersten 1945, 
Néel 1946) und die Anfangspermeabilität (Kersten 1955) zwei unterschiedliche Ursachen, 
und zwar einerseits in der Gegenwart unvermeidlicher, nicht-ferromagnetischer 
Inklusionen (Н, resp. и), andererseits in der Gegenwart von unvermeidlichen inneren 
Spannungen (H, resp. и). Es ist naheliegend, die gesamte Koerzitätskraft H, resp. 
die gesamte Anfangspermeabilität u, auf folgende Weise zu definieren: 


H, = Н, + Н, (1) 

Hp, uli echa, ue | 2 (2) 

Die plastische Dehnung, die uns wenigstens angenähert die Möglichkeit gibt, 

die Spannungsanteile von Koerzitätskraft und inverser Anfangspermeabilität zu 
ermöglicht uns zu entscheiden, ob für die Koerzitütskraft die gleichen Span- 


į (165) 


steigern, 
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nungsanteile in Betracht kommen, wie für die Koerzitatskraft, oder nicht. Bild 3 scheint 
für Nickel zu beweisen, dass es sich in diesem Falle um die gleichen Ursachen handelt. 
Weitere, jedoch nur vorläufige, Messungen an Mumetall haben uns gezeigt, dass auch 
hier der Zusammenhang linear bleibt (Fig. 4). Für kommerzielles Eisen und Hipernik 
gilt jedoch nach unserem bisherigen Ergebnis der lineare Zusammenhang nicht. 


L-Lora; Ro +0 


(rot AT)? (Ko -(о + X) 
A= RO -lor lro FATT "Co 
P TQ + AT г 
р w lko -Co* X) rota 75257 Ro 


Fig. 1 
1) Messmuster, 2) M fur ein Marea von etwa 40 Millioersted, 3) Messkondensator, 
4) Stromquelle, 5) induktionsfreier Messwiderstand, 6) konstante Kapazitüt, 7) Präzisionsübertrager 
(2 + 2) :1, 8) Resonanzverstärker, 9) Kathodenstrahlróhre als Nullinstrument, 
Fig. 2. 
1) Messmuster, 2) Feldspule max. 75 Oersted, 3) Kompensationsspule gleicher Dimensionen, 
4), 5), 6) Magnetsonde aus Permalloydraht, (siehe Literaturangabe 6), 7) ihre Tertiärwicklung, 
8) Verstärker, 9) Kathodenstrahlröhre als Nullinstrument, 10), 11) Stromquellen, 13) hier nicht 
benutzer magnetometrischer Teil des Messgerütes, bestehend aus Stromspule der Länge des Mess- 
"~~~ musters, mittels derer sich dessen Magnetisierung kompensieren lässt. 
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DOMAIN FORMATION IN FERROELECTRICS 


Ву ARKADIUSZ JASKIEWICZ 


Physical Laboratory B. Bierut University, Wroclaw 


( Received Semptember 22, 1956) 


The dielectric properties of ferroelectries depend on the domain structure and 
domain dynamics. The formation and movement of domains have been studied by 
various authors: Forsbergh (1949), Merz (1954), Little (1955). Our main aim is to find 
how the formation of the c domain of BaTiO depends on the electric field. 

The nucleation theory supposes that local fluctuations in free energy produce 
small, ordered regions called embryos. The free energy has to overcome various energy 
expenditures. Owing to the increase in the surface energy and other expenditures 
associated with the formation of an embryo, most embryos are unstable and vanish. 
Embryos greater than a critical size are stable and grow with a decrease in the free 
energy of the substance. This implies that a critical size can be reached only by the 
kinetic process of statistical fluctuations, see e.g. R. Smoluchowski (1951). 

The change in free energy AF associated with the formation of a nucleus 
in ferroelectrics can be represented by the expression 


AF = AF, = AF, + AF, + AF, (1) 


The first term represents the surface energy (wall energy); the second, the 
volume energy (dependent on the parameter of order); the third, the energy of inter- 
action between the polarized dielectric and the field; the last term, the depolarisation 
energy. 

If we take the spontaneous polarization as the parameter of order, then for the 
order energy we obtain the expression 4 


AF, = — «Р? (2) 
where а depends on the shape of the embryo, V is the volume, and P, is the spontaneous 
polarization. 

Formula (1) giving the free energy change in the absence of a field, can be written 
in the form: 
AF = oS — «РУ + ВР? l (3) 
where ø is the wall energy/cm?, S the wall area, В the depolarizing factor. 
/ (169) 
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If a phase transformations is to occur, @ must be greater than В, because the free 
energy of the new phase is lower than that of the metastable phase. 
We assume that the embryo has а cone-like shape with length J and radius r in 


accordance with optical observations Hippel (1950), Merz (1954). Using (3) we have 


1 
AF = ло (r? + п) — (a — В) x a LP (4) 
if we write 
х— В = у (5) 
we get 
AF = пс (r? + rl) — 5 yarl P. (6) 
We obtain the critical dimensions from the maximum of (6): 
DELETE T 
[^ УР? , с yP? ( ) 


We see that there will be only large domains about the transition point. The 
nucleation rate of domains we can get from well known expression for the rate of 
nucleation. 

Now, we take under consideration a ferroelectric crystal in an external field Е. 
The free energy change of nucleus formation is given by (1), thus, from (6) we can 
write: 


АЕ? = по (® +r) — 5 yari P?— ЕР, (8а) 
ДЕ = ло (r? + rl) a yarlP + EP,V (8b) 


where — EP V is the electrostatic energy of a nucleus with parallel polarization, 
EP,V the energy of a nucleus with antiparallel polarization. 

If we repeat our procedure we obtain the critical dimensions: 

a) domains with parallel polarization 


30 60 


FP, WP, E) P, (yP, TE) Өз) 
b) domains with antiparallel polarization 
a 30 a 60 
в =; L= 9b 
P, (yP, — E) P, (yP, — E) s 


We see that equations (9a) and (9b) explain the behaviour of domains in the electric 
field. In the case of a nucleus with antiparallel polarization, the critical dimensions 
are greater than in the case of parallel polarization. АП the antiparallel nuclei which 
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have dimensions less than the critical one are unstable and vanish. The destroyed 
antiparallel domains have a greater probability to nucleate parallelly than antiparallely. 
The probability may be obtain from the nucleation rate, The last equations are well 
suited for a study of domain formation in connection with experiments such as those 
of Piekara and Pająk (1953) and of Badian and Jaśkiewicz. 

The author wishes to express his thanks to Dr. Z. Galasiewicz and Professor 
B. Makiej for helpful discussion 
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